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The fiber bundle model describes a collection of elastic fibers under load. The fibers fail sucessively 
and for each failure, the load distribution among the surviving fibers changes. Even though very 
simple, this model captures the essentials of failure processes in a large number of materials and 
settings. We present here a review of the fiber bundle model with different load redistribution 
mechanisms from the point of view of statistics and statistical physics rather than materials 
science, with a focus on concepts such as criticality, universality and fluctuations. We discuss 
the fiber bundle model as a tool for understanding phenomena such as creep, and fatigue, how it 
is used to describe the behavior of fiber reinforced composites as well as modelling e.g. network 
failure, traffic jams and earthquake dynamics. 
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I. INTRODUCTION 

In materials science and engineering, a class of sim- 
ple models, known as fiber bundle models (FBM), has 
proven to be very effective in practical applications such 
as fiber reinforced composites. In this context, su ch mod- 
els ha ve a history that goes back to the twenties ( Peircd . 
Il926l ). and they constitute today an elaborate toolbox 
for studying such materials, rendering computer stud- 
ies orders of magnitudes more eff icient than brut e force 
methods. Since the late eighties (ISornettel . [l989l) . these 
models have received increasing attention in the physics 
community due to their deceivingly simple appearance 
coupled with an extraordinary richness of behaviors. As 
these models are just at the edge of what is possible an- 
alytically and typically not being very challenging from 
a numerical point of view so that extremely good statis- 
tics on large systems are available, they are perfect as 
model systems for studying failure phenomena as a part 
of theoretical physics. 

Fracture and material stability has for practical rea- 
sons interested humanity ever since we started using 
tools: our pottery should be able to withstand handling, 
our huts should be able to withstand normal weather. As 
science took on the form we know today during the Re- 
naissance, Leonardo da Vinci studied five hundred years 
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ago experimentally the strength of w ires — fiber bun- 
dles — as a function of their length ( Lund and Bvrnd . 
l200l[ ). Systematic strength studies, but on beams, were 
also pursued systematic by Galileo Galilei one hundred 
years later, as was done by Edme Mariotte (of gas law 
fame) who pressurized vessels until they burst in connec- 
tion with the construction of a fountain at Versailles. For 
some reason, mainstream physics moved away from frac- 
ture and breakdown problems in the nineteenth century, 
and it is only during the last twenty years that fracture 
problems have been studied within physics proper. The 
reason for this is most probably the advent of the com- 
puter as a research tool, rendering problems that were 
beyond the reach of systematic theoretical study now ac- 
cessible. 

If we were to single out the most important contribu- 
tion from the physics community with respect to fracture 
phenomena, it must be the focus on fluctuations rather 
than averages. What good is the knowledge of the aver- 
age behavior a system when faced with a single sample 
and this being liable to breakdown given the right fluctu- 
ation? This review, being written by physicists, reflects 
this point of view, and hence, fluctuations play an im- 
portant role throughout it. 

Even though we may trace the study of fiber bundles 
to Leonardo da Vinci, their m odern s tory s tarts with 
the already mentioned work bv iPeirci ( 1926f ). In 1945, 
Daniels published a seminal review cum research arti- 
cle on fiber bundles which still t oday must be regarded 
as essential reading in the field (lDanielsl . [l945h . In this 
paper, the fiber bundle model is treated as a problem 
of statistics and the analysis is performed within this 
framework, rather than treating it within materials sci- 
ence. The fiber bundle is viewed as a collection of elastic 
objects connected in parallel and clamped to a medium 
that transmits forces between the fibers. The elongation 
of a fiber is linearly related to the force it carries up to 
a maximum value. When this value is reached, the fiber 
fails by no longer being able to carry any force. The 
threshold value is assigned from some initially chosen 
probability distribution, and do not change thereafter. 
When the fiber fails, the force it carried is redistributed. 
If the clamps deform under loading, fibers closer to the 
just-failed fiber will absorb more of the force compared 
to those further away. If the clamps, on the other hand, 
are rigid, the force is equally distributed to all the surviv- 
ing fibers. Daniels discussed this latter case. A typical 
question posed and answered in this paper would be the 
average strength of a bundle of N fibers, but also the vari- 
ance of the average strength of the entire bundle. The 
present review takes the same point of view, discussing 
the fiber bundle model as a statistical model. Only in 
Section V we discuss the fiber bundle model in the con- 
text of materials science with all the realism of real ma- 
terials considered. However, we have not attempted to 
include any discussions of the many experimental studies 
that have been performed on systems where fiber bun- 
dles constitute the appropriate tool. This is beyond the 



scope of this statistical-physics based review. 

After introducing (Section II) the fiber bundle model 
to our readers, in Section III, we present the Equal Load 
Sharing Model, which was sketched just a few lines back. 
This seemingly simple model is in fact extremely rich. 
For example, the load at which catastrophic failure oc- 
curs is a second order critical point with essentially all 
the features usually seen in systems displaying such be- 
havior. However in this case, the system is analytically 
tractable. In fact, we believe that the equal load shar- 
ing fiber bundle may be an excellent system for teaching 
second order phase transitions at the college level. Un- 
der the heading of Fluctuations, we discuss the burst 
distribution, i.e., the statistics of simultaneously failing 
fibers during loading: When a fiber fails and the force it 
was carrying is redistributed, one or more other fibers 
may be driven above their failing thresholds. In this 
equal load sharing model, the absolute rigidity of the 
bar (transmitting forces among the fibers) suppresses the 
stress fluctuations among the fibers. As such, there is no 
apparent growth of the (fluctuation correlation) length 
scale. Hence, although there are precise recursion rela- 
tions and their linearized solutions are available near the 
fixed point (see Sec. HI), no straight forwa rd appl i cation 
of the Renormalization group techniques ( Fisherj . Il974l ) 
has been made to extract the exponents through length 
scaling. 

In Section IV Local Load Sharing is discussed. This 
bit of added realism comes at the added cost that an- 
alytical treatment becomes much more difficult. There 
are, still, a number of analytical results in the literature. 
One may see intuitively how local load sharing compli- 
cates the problem, since the relative positions of the fibers 
now become important. Under global load sharing, ev- 
ery surviving fibers gets the same excess force and, hence, 
where they are do not matter. There are essentially three 
local load sharing models in the literature. The first one 
dictates that the nearest surviving neighbors of the fail- 
ing fiber absorb its load. Then there are "softer models" 
where the redistribution follows a power law in the dis- 
tance to the failing fiber. Lastly, there is the model where 
the clamps holding the fibers are elastic themselves, and 
this leads to non-equal redistribution of the forces. 

Section V contains a review of the use of fiber bun- 
dle models in applications such as materials science. We 
discuss fatigue, thermal failure, viscoelastic effects and 
precursors of global failure. We then go on to review 
the large field of modeling fiber reinforced composites. 
Here fiber bundle models constitute the starting point of 
the analysis, which, by its very nature, is rather complex 
seen from the viewpoint of statistical physics. Lastly, we 
review some applications of fiber bundle models in con- 
nection with systems that initially would seem quite far 
from the concept of a fiber bundle, such as traffic jams. 

We end this review by a summary with few concluding 
remarks in Section VI. 
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II. FIBER BUNDLE MODELS 



Imagine a heavy load hanging from a rigid anchor point 
(say, at the roof) by a rope or a bundle of fibers. If the 
load exceeds a threshold value, the bundle fails. How 
does the failure proceeds in the bundle? Unless all the 
fibers in the bundle have got identical breaking thresholds 
(and that never happens in a real rope), the failure dy- 
namics proceeds in a typical collective load transfer way. 
One can assume that in this kind of situation the load is 
equally shared by all the intact fibers in the bundle. How- 
ever, the breaking threshold for each of the fibers being 
different, some fibers fail before others and, consequently, 
the load per surviving fiber increases as it gets redis- 
tributed and shared equally by the rest. This increased 
load per fiber may induce further breaking of some fibers 
and the avalanche continues, or it stops if all the surviv- 
ing fibers can withstand the redistributed load per fiber. 
This collective or cooperative failure dynamics and the 
consequent avalanches or bursts are typical for the failure 
in any many-body system. It captures the essential fea- 
tures of failure during fracture propagation (recorded by 
acoustic emissions) , earthquake avalanches (main and af- 
tershocks), traffic jams (due to dynamic clustering), etc. 

The model was first introduced in 1926 by iPeircd 
(I1926D m the context of textile engineering. Since then 
it was modified a little and investigated, mainly numeri- 
cally, with various realistic fiber threshold distribut ions 
by th e engineering communi t v ICoiemanl.ll957al:fDaniel£ , 
1945 : iHarlow and Phoeni j . Il978l: iPhoenix and Smith . 
19831 ). Starting from late eighties, physicists took 
interest in the a valanche distribut ion in the model and in 
its d yn amics ([Andersen et al. I . 



1993; 



1997 



1991 



Kun et 



Hemmer and Hansen 
20001: 



■ 19971 : 
19921: 



Gomez et 



Kloster et al 



Pradhan et all 



iNewman and Gabrielov, 
20021 : iSornette , T989 , 1992 ; 



Zapperi et all Il997l ). A recursive dynamical equa- 



tion was set up fo r the equal-load-sharing version 
recen tly ( de Silveiral . Il999l : iPradhan and Chakrabartil 
l200lh an d the dynamic critica l behav io r is now solved 

exactly (iBhattacharvva et al. j 20031 : iPradhan et all 

l2002t IPradhan and Hemmei 7120071) . In addition to the 



extensive numeri c al res ults ( Hansen and Hemmed . Il994l : 
IZhang and Di^ . I1994D on the effect of short-range 
fluctuations (local load sh aring cases), some p r ogress 



1994 : 



-P 

with analyt i cal studies (Duxbury a nd Leathl. 
Gomez et al ]. [1993.: Harlo w., .1985. : .Harlow and Phoenixl 
19911 : iKloster et aLl Il997l)" have also been made 



There are a large number of experimental studies of 
various materials and phenomena that have successfully 
been analyzed within t he framework of th e flber bun- 
dle model. For example. iLavton and Sastrvl ( 2004J ) have 
used the fiber bundle model to propose explanations for 
changes in fibrious collagen a nd its relation to neuropathy 
in connection with diabetes. iToffoli and LehmanI (I2001D 
propose a method to monitor the structural integrity of 
fiber-reinforced ceramic-matrix composites using electri- 
cal resistivity measurements. The basic idea here is that 



FIG. 1 A fiber bundle model having parallel fibers. The 
original position of the rigid platform on which force has been 
applied, is indicated. All the fibers are assumed to have the 
same elastic constant (normalized to unity here) until break- 
ing, while the breaking strength of these fibers are assumed to 
be randomly dispersed. As the bundle gets strained (a; > 0), 
some weaker fibers (having lower stress or strain capacity) 
fail, and the number of intact fibers decrease from its starting 
value N. Consequently, the stress on them (or the strain) in- 
creases due to the redistribution or transfer of loads from the 
failing fibers. Some of these fibers may fail further as they can 
not support this extra redistributed load. The process stops 
if there is no further failure, and the bundle will show non- 
linear elastic response (although each fiber has linear elastic 
behavior until breaking). Otherwise, the bundle fails when 
the stress concentration on the fibers (due to the dynamic 
stress redistribution) becomes so high that none of the fibers 
can withstand that. 



when the fibers in the composite themselves fail rather 
than just the matrix in which they are embedded, the 
structure is about to fail. The individual fiber failures is 
recorded through changes in the electrical conductivity 
of the material. Acoustic emission, the crackling sounds 
emitted by materials as they are loaded, provide yet an- 
other example where fiber bundle models play and im- 
portant role, see e.g. jNechad et al. I (120051 . 



III. EQUAL LOAD SHARING MODEL 

The simplest and the oldest version of the model is 
the equal load sharing (ELS) model, in which the load 
previously carried by a failed fiber is shared equally 
by all the remaining intact fibers in the system. As 
the applied load is shared globally, this model is also 
known as global load sharing (GLS) model or demo- 
cratic fiber bundle model. Due to the consequent mean- 
field nature, some exact results could b e extrac t ed for 
this model and this was demonstrated bv lDanielsl (|1945D 
in a classic work some sixty years ago. The typical 
relaxation dynamics of this model has been solved re- 
cently which has clearly established a robust critical be- 
havior (IBhattacharvva aLl. [20031: [P radhan et "al[. l2002l: 



IPradhan and Hemmer . 2007). It may be mentioned at 
the outset that the ELS or GLS models do not allow 
for spatial fluctuations (due to the absolute rigidity of 
the platform in Fig. 1) and hence such models be- 
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long to the the mean field category of critical dynam- 
ics, see e.g., IStanlevI (I1987D . Fluctuations in breaking 
time or in avalanche statistics (due to randomness in fiber 
strengths) are of course possible in such models and are 
discussed in details in this section. 

A bundle can be loaded in two different ways: Strain 
controlled and force controlled. In the strain controlled 
method, at each step the whole bundle is stretched till 
the weakest fiber fails. Clearly, when number of fibers N 
is very large, strain is increased by infinitesimal amount 
at each step until complete breakdown and therefore the 
process is considered as a quasi-static way of loading. 
On the other hand, in the force controlled method, the 
external force (load) on the bundle is increased by same 
amount at each step until the breakdown. The basic 
difference between these two methods is that the first 
method ensures the failure of single fiber (weakest one 
among the intact fibers) at each loading step, while in 
the second method sometimes none of the fibers fail and 
sometimes more than one fail in one loading step. 

Let X denote the strain of the fibers in the bundle. 
Assuming the fibers to be linearly elastic up to their re- 
spective failure point (with unit elastic constant), we can 
represent the stresses on each of the surviving fibers by 
the same quantity x. The strength (or threshold) of a 
fiber is usually determined by the stress value x it can 
bear, and beyond which it fails. We therefore denote 
the strength (threshold) distribution of the fibers in the 
bundle by p{x) and the corresponding cumulative distri- 
bution by P{x) = p{y)dy. Two popular examples of 
threshold distributions are the uniform distribution 



P(x) = 



1 



for < X < Xr 
for X > Xr, 



and the WeibuU distribution 



P{x) = 1 - exp(-(a;/a;^)''). 



(1) 



(2) 



Here reference threshold, and the dimensionless 

number p is the Weibull index (Fig. [2]). 

In the strain controlled loading, at a strain x, the total 
force on the bundle is x times the number of intact fibers. 
The expected or average force at this stage is the refore 



ine expectea or average torce at tnis stage is tne i 
(jHemmer and HanseDl . ll992l: ISornetti . ll989l[T993 ) 

F{x) = Nx{l- P{x)). (3) 

The maximum Fc of F(x) corresponds to the value Xc for 
which dF/dx vanishes: 



1 - P{Xc) ~ Xcp{Xc) = 0. 



(4) 



Here the failure process is basically driven by fluctua- 
tions an d can be analyzed using extreme order statis- 
tics (iHemmer an ^ Hanseni . Il992 : iKloster et al. I, Il997t 
ISornettel . ll989lll992D . 

In the force controlled method, if force F is applied on 
a bundle having N fibers, when the system reaches an 
equilibrium, the strain or effective stress x is (see Fig. [T]) 



x{F) 



F 



N[l~P{x)]' 



(5) 



Q. 



(A) 



Q. 




FIG. 2 The uniform distribution (A) and the Weibull distri- 
bution (B) with p = 5 (solid line) and p = 10 (dotted line). 



Therefore, at the equilibrium state, Eq. ^ and Eq. ^ 
are identical. It is possible to construct recursive dyna m- 



ics (|de Silveiral . ll999l : lPradhan and Chakrabartill200lh of 

the failure process for a given load and the fixed-point so- 
lutions explore the average behavior of the system at the 
equilibrium state. 



A. Average behavior 

Fig. [T] shows a static fiber bundle model in the ELS 
mode where N fibers are connected in parallel to each 
other (clamped at both ends) and a force is applied at 
one end. At the first step all fibers that cannot withstand 
the applied stress break. Then the stress is redistributed 
on the surviving fibers, which compels further fibers to 
break. This starts an iterative process that continues 
until an equilibrium is reached, or all fibers fail. The 
average behavior is manifested when the initial load is 
macroscopic (very large N). 



1. Recursive breaking dynamics 

The breaking dynamics can be repre- 
sented by recursion r el ations (|de Silveiral . Il999l: 
iPradhan and Chakrabartil l200l[ ) in discrete steps. 
Let Nf be the number of fibers that survive after step 
t, where t indicates the number of stress redistribution 
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steps. Then one can write ( de Silveiral . Il999l ) 



Nt+i = N 



Nt 



(6) 



Now we introduce a = F/N, the appUed stress and 
Ut = Nt/N, the surviving fraction of total fibers. Then 
the effective stress after t step becomes Xt = y- and 
after t + 1 steps the surviving fraction of total fibers 
is Ut+i = 1 — P{xt)- Theref ore we can construct 
the following recursion re l ation s ( Pradhan et all l2002t 
iPradhan and ChakrabartiL l200lh : 



2^4+1 = 



l-Pixt) 



and 



Ut+i = l-Pia/Ut);Uo = l. 



(7) 



(8) 



At equilibrium Ut+i = Ut ^ U* and Xt+i = Xt = x* . 
These equations (Eq. [7] and Eq. [S]) can be solved at 
and around the fixed points for the particular strength 
distribution p{x). 



2. Solution of the dynamics: Critical behavior 

Let us choose the uniform density of fiber strength dis- 
tribution (Eq. [T]) up to the cutoff Xr ~ 1. Then the cu- 
mulative distribution becomes P{a/Ut) = cr/Ut- There- 
fore from Eq. (O and Eq. ([8]) we can construct a pair of 
recursion relations 



Xt+l 



I - Xt 



and 



Ut 



t+i 



1 



a 

TTt' 



(9) 



(10) 



This nonlinear recursion equations are somewhat charac- 
teristic of the dynamics of fiber bundle models and such 
dynamics can be obtained in many different ways. For 
example, the failed fraction 1 — Ut+i at step 1 is given 
by the fraction F/NUt = cr/Ut of the load shared by the 
intact fibers at step t and for the uniform distribution of 
thresholds (Fig. [HA), one readily gets Eq. PH)) . 

At the fixed point the above relations take the 
quadratic forms 



and 







U* -U* + (j = Q, 



with the solutions 



x*{a)^-±{a,-af'^ 



(11) 



(12) 



(13) 



u 



1 

0.8 
0.6 
0.4 
0.2 


0.05 0.1 0.15 0.2 0.25 0.3 

a 

FIG. 3 Variation of effective stress (x*), fraction of unbroken 
fibers {U*) and the order parameter (O) with the applied 
stress cr for a bundle with uniform distribution (Eq. [T| of 
fiber strengths. 



and 



U*{a) = -±{a, 



Nl/2 



(14) 



Here ac~\\s the critical value of applied stress beyond 
which the bundle fails completely. Clearly, for the effec- 
tive stress (Eq. [T5|) solution with (— ) sign is the stable 
fixed point and with (+) sign is the unstable fixed point 
whereas for fraction of unbroken fibers (Eq. [T4|) . it is 
just the opposite. Now the difference U*{(j) — U*{ac) 
behaves like an order parameter signaling partial failure 
of the bundle when it is non-zero (positive), although 
unlike conventional phase transitions it does not have a 
real- valued existence for a > Or- 



0=U*{a)-U*{ac) = (a^^aY 



(15) 



Fig. [3] shows the variation of x* , U* and O with the 
externally applied stress value a. One can also obtain 
the breakdown susceptibility Xi defined as the change of 
U*{a) due to an infinitesimal increment of the applied 
stress cr: 



dU*ia) 



da 



(16) 



Such a divergen ce in \ ha d alread y been reported in 
several studies (Ide Silveiral . 11999, : .Moreno et all \200(k 
IZapperi ei"aZI . Il997lll999al) . 

To study the dynamics away from criticality (cr —^ ac 
from below), the recursion relation (Eq. [T0|) can be re- 
placed by a differential equation 



dU_ 
'dt 



U'-U 



U 



Close to the fixed point, Ut{<j) = U*{(j) 
AJ7 0) and this gives 



(17) 

AC/, (where 



A(7 = Ut{cT) - U*{a) « exp(-t/r). 



(18) 
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where t = ^ [^('''c ~ + l] • Therefore, near the 

critical point: 

T^{a,-ar'-0 = ^. (19) 

At the critical point (cr = Uc), a dynamic critical behavior 
has been observed in the relaxation of the failure process 
to the fixed point. From the recursion relation (Eq. 
it can be easily verified that the fraction Ut{ac) follows a 
simple power-law decay: 



+ (20) 

starting from Uq ^ 1. For large t (t — > cxd), this reduces 
to [/( — 1/2 oc t^''; rj = I; indicating critical slowing down 
which is a robust characterization of the critical state. 



3. Universality class of the model 



(a) (b) 




FIG. 4 (a) The density function p(x) and (b) the cumulative 
distribution P(x) of random fiber strengths x distributed with 
linearly increasing density in the interval [Cl,Cr]. In the 
particular instance shown in the figure Cl = 0.15 and Cr = 
0.75. 



The critical properties, obtained above, are for the uni- 
form threshold distribution, and the natural question is 
how general the results are. To check the universality of 
the ELS model, two other types o f fiber strength distri- 
butio ns can be easily considered (iBhattach arvva et all . 
|2Q03[ ): linearly increasing density distribution and lin- 
early decreasing density distribution. 

For linearly increasing density of fiber strengths in the 
interval [Cl^Ch], the normalized density function and 
the cumulative distribution are given by (illustrated in 
Fig.©: 



and 

Ut+i = l-(J^^~ri}j , Uo^l. (25) 

The fixed point equations, (Eq. [11] and Eq. (T^]) , now 
assume cubic form: 

(F*)^ - 2Fl (F*)^ + (fI - l) F* + Fo = (26) 
where F* = x* /{gr — ctl), and 



p{x) 




(21) 



and 



P(x) 



0, < X < Cl 

1, Cr < X. 



(22) 



Now we introduce the transformed quantities: 



Cr — Cl ' 



Fl = 



Cl 



Cr — Cl ' 



Ft = 



Cr — Cl 
(23) 

For an initial stress Cl < cr < Cr, (or, F^ < Fq < 
Fl + 1) along with the cumulative distribution given by 
Eq. ([22| . the recursion relations (Eq. [7|and Eq. [8]) appear 
as: 



F 



t+i 



To 

l-(Ft-FL)2 



(24) 



iu*y + (Fi - 1) {u*y - (2FlFo) u* + tI^ o. (27) 

Consequently each of the recursions (Eq. [24] and Eq. [25]) 
have three fixed points - only one in each case is found 
to be stable. For the redistributed stress the fixed points 



where 



and 



2 $ 

-Tl + 2K0COS-, 

2 $ $ 

-Fl - Ko cos — + VSKo sm — , 

2 $ $ 

-Fl ~ ifocos— - VSKosm—, 



^0 = gVs + ri 



COS* = 



(9 - rj) - 27Fo/2 
(3 + Fi)^/^ 



(28) 
(29) 
(30) 

(31) 
(32) 
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(a) (b) 




FIG. 5 The fixed points of (a) the redistributed stress, and 
(b) the surviving fraction of fibers for the distribution of fiber 
strengths shown in Fig. [J] In each part of the figure the curve 
for the stable fixed points is shown by a bold solid line and 
those for the unstable fixed points are shown by bold broken 
lines. We have Cl = 0.15 and Cr — 0.75, so that Uc = 0.3375; 
the position of the critical point is marked by an arrowhead. 
For a < Cl the fixed points are trivial: since there are no 
broken fibers x* = a and U* = Uo = 1. 

Similarly, for the surviving fraction of fibers the fixed 
points are: 



and then each of the quantities T and U have one stable 
and one unstable fixed point. The critical point has the 
trivial lower bound: CTc > Cl- The expression of Tc in 
Eq. ([55)1 shows that it approaches the lower bound as 
Ti — > cxD which happens for finite values of C'l and 
when {Cii — Cl) — > 0. It follows that the upper bound 
for the critical point is also trivial: Uc < Cr. Also, at 
the critical point we get from Eq. ([5^ and Eq. ([57)1 : 

COS$crit = COSecrit = -1 (39) 

or, 

$crit = Ocrit = TT. (40) 

The stable fixed points r2 and Ui are positive real- 
valued when Fq < Tc', thus the fiber bundle always 
reaches a state of mechanical equilibrium after partial 
failure under an initial applied stress (Tq < Cc- For a > (Tc 
(or, Fo > Fc), F2 and C/j* are no longer real- valued and 
the entire fiber bundle eventually breaks down. The tran- 
sition from the phase of partial failure to the phase of 
total failure takes place when cr just exceeds ctc and the 
order parameter for this phase transition is defined as in 
Eq. ini): 

O ^ [/* - C/r_„it- (41) 

Close to the critical point but below it, we can write, 
from Eq. §7^ and Eq. (001), that: 



c/r = 

m = 



e 



^-^ + 2Jocos-, 

1 - F^ 6 
— ^— ^ - Jo cos — + \/3 Jo sii. ^ , 

— Jo cos — - V3 Jo sm — , 



where 



and 



Jo = ^^{Tl-ir + 6TLTo 



(33) 
(34) 
(35) 

(36) 



cos 9 



(ri - 1) 



9FlFo 



27rg/2 



(F2 - if + 6FlFo 



3/2 



(37) 

Of these fixed points F2 and Ui are stable whereas F^, 
and U2, C/g are unstable (Fig.El). 
The discriminants of the cubic equations (Eq. and 
Eq. [27]) become zero at a critical value Cc (or, Fc) of the 
initial applied stress: 



TT — ~ sin O 

3V3Fc(3 + Fi)3/4(F,-Fo)i/^ 
[(Fi-l)2+6FiFd3/2 



(42) 



and the expressions for the fixed points in Eq. ([55]) and 
Eq. ([M[) reduce to the forms: 



rc(3 + ri)3/4 ^^^_^^^i/2 (43) 



(Fi - 1)2 + 6FiFc 



and 



2--2-crit (r2_l)2+6FLF 

where 



TJ* — TJ* 
•^l-crit — '^2-crit 



1 - F2 1 



-^ + ^^(Fi-l)2 + 6FiF, 



(45) 

is the stable fixed point value of the surviving fraction 
of fibers under the critical initial stress ac- Therefore, 
following the definition of the order parameter in Eq. (1411) 
we get from the above equation: 



2 



27 



Ft 9 - Ff 



(3 + F 



2\3/2 
L) 



(38) 



Fc(3 + Fi)3/4 _ 
^ " (Fi - 1)2 + 6FiF/^ ■ 



To 



(46) 
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On replacing the transformed variable Tq by the original 
a, Eq. shows that the order parameter goes to zero 
continuously following the same power-law as in Eq. (fT5|) 
for the previous case when a approaches its critical value 
from below. 

Similarly the susceptibility diverges by the same 
power-law as in Eq. p6p on approaching the critical point 
from below: 



X 



d[/i* 



da 



cx (Tc - To) 



-1/2 



(47) 



The critical dynamics of the fiber bundle is given by 
the asymptotic closed form solution of the recursion (Eq. 
[25]) for To = r^: 



Ut-Ul 



crit 



(C/r-crit)' 



3(r,)^-2rir,c/*_„.i, 



(48) 

where Tc and U^^^j.^^ are given in Eq. ([55]) and Eq. (liS)) 
respectively. This shows that the asymptotic relaxation 
of the surviving fraction of fibers to its stable fixed point 
under the critical initial stress has the same (inverse of 
step number) form as found in the case of uniform density 
of fiber strengths fEq. l^ . 

We now consider a fiber bundle with a linearly decreas- 
ing density of fiber strengths in the interval [C l , C'^j] . 
The normalized density function and cumulative distri- 
bution (illustrated in Fig. (5]) are: 



0, 



p{x) 



< X < Cl 

Cr<X 



2{Cr-x) Ct < rr < 
0, 



(49) 



and 



P(x) 



Cr-x 
Cr — Cl 



0<x<Cl 
Cl<x<Cr 

Cr < X 



(50) 



With the transformed quantities defined in Eq. 
the recurrences (Eq. [7] and Eq. [5]) for < cr < Cr 
appear as: 



(a) 



(b) 





FIG. 6 (a) The density function p(x) and (b) the cumulative 
distribution P{x) of fiber strengths x distributed with hnearly 
decreasing density in the interval [Cl,Cr]. Similar to the 
cases shown in Fig. |4]we have Cl = 0.15 and Cr = 0.75 in 
this example also. 



(C/*)^-(l -I- Fl)^ {U*f + 2 (1 + Tl) FoC/*-r2 = (54) 
and they have the following solutions: 



^(l + ri) + 2i^^cosy, 



(55) 



^(H-ri)-X^cosy + \/3i^'siny, (56) 
^{1 + Tl) - K[,cos^ - VSK' sin^, (57) 



where 



^0 



cos$' 



l + T, 



27Fo 



2(i + ri)' 



(58) 



(59) 



and 



(I + Fl - Ft) 



and 



Ut+i 



(51) 



Uo = 1. (52) 



The fixed point equations are again cubic: 
(r*)3 -2(1 + Fl) (F*)' + (1 + Fl)' F* - Fo = 0, (53) 



u: = 



where 



ii±I^ + 24cosf 



(60) 



V cosy + V3J^ sin y, (61) 

^ -J^osy-^/34siny, (62) 



•/o^ 3^(1 + Tl)*- 6(1 + ri)ro, (63) 
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(a) 




(b) 













uv 



0.05 0.1 0.16 t 0.2 
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JT* — JJ* 

'-'l-crit — '^2-crit 



4 9 



and 



cos e;^,it 



cos $' .it = 1 



or, 



e' 



(67) 
(68) 
(69) 



Comparing Eq. 



and Eq. with Eq. (gD]) we see 



that the critical values of 8 and 8' 



the same whereas 



those of $ and differ by tt radians. 

Near the critical point, but below it, we get from 
Eq. ^ and Eq. 



FIG. 7 The fixed points of (a) the redistributed stress, and 
(b) the surviving fraction of fibers for the distribution of 
fiber strengths shown in Fig. |6] The curve for the stable 
fixed points is shown by a bold solid line and those for the 
unstable fixed points are shown by bold broken lines. In 
this example too we have Cl = 0.15 and Cr = 0.75; here 
a"c = 0.173611, marked by an arrowhead. The critical point is 
located lower than that in Fig. [5] due to abundance of fibers 
of lower strengths compared to the previous case. 



(i + Fl)^ (i + rL)^-9ro +27rV2 

cos 8' = ^ i — . (64) 



(i + rL)^-6(i + ri)ro 



3/2 



Here Fg and Ul are stable fixed points while the rest are 
unstable (Fig.El). 

The discriminants of Eq. ((53)) and Eq. ([54]) show that 
the critical applied stress in this case, a'^ (or, F^), is given 
by: 



F' 



or, 



21{Ch-Cl) 



2 ■ 



(65) 



(66) 



In order to satisfy the condition a'^ > C'l , it requires from 
Eq. ([66]) that Cr > 3Cl, which imposes an upper bound: 

Uf. ^ 3 • 

Like before, for Fq < F^ the stable fixed points are 
real-valued, which indicates that only partial failure of 
the fiber bundle takes place before a state of mechanical 
equilibrium is reached; for Fq > F^ the fixed points are 
not real and a phase of total failure exists. The order 
parameter O of the transition is given by the definition 
in Eq. gH). 

For Fq = F^ we get the following properties from 
Eq. dSlD, Eq. ® , Eq. (EU and Eq. 



,1/2 



and 



m ~ u. 



2-crit 



Fi)^/' (F^-Fc 



,1/2 



(70) 



(71) 



Therefore, by the definition of the order parameter in 
Eq. (|¥T|l and that of the susceptibility in (|T7)l . we get in 



this case O cx (F;^ - Tq)^^^ and x oc (F^^ - Fq) 



-1/2 



ro 



FJ,— . These power laws have the same exponents as the 
corresponding ones in the previous cases and differ from 
those only in the critical point and the critical amplitude. 

At the critical point the asymptotic relaxation of the 
surviving fraction of fibers to its stable fixed point [ob- 
tained as an asymptotic solution to Eq. ([5^ ] is again 
found to be a power law decay similar to Eq. (j20p and 
Eq. 



4 ^r-crit 

3 t ' 



(72) 



Ut - C/r-crit 



The two density functions, Eq. and Eq. (|49|) . 

can be transformed from one to the other by a reflec- 
tion on the hne x = {Cl + Cr)/2 (compare Fig. (D^a) 
and Fig. [Bl^a)). But the fixed point equations and their 
solutions do not have this symmetry. This is because 
the density function p{x) does not appear directly in 
the recursion relations for the dynamics. It is the cu- 
mulative distribution P{x) which appears in the recur- 
sion relations. Eq. ([22)) and Eq. ([50)) show that the cu- 
mulative distributions of these two cases are not mutu- 
ally symmetric about any value of the threshold stress 
X (compare Fig. IHb) and Fig. IHlIb)). However a cer- 
tain relation exists between the critical values of the 
applied stress for a special case of these two models: 
if Cl = 0, we get from Eq. ^ and Eq. ^ that 
Oc/Cr = a/4/27 and o'JCr = 4/27 respectively; there- 
fore we have (j'^/Cr = [(Jc/Cr)^ . 
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The critical behavior of the models discussed in this 
section show that the power laws found here are in- 
dependent of the form of the cumulative distribution 
P. The three threshold distributions studied here have 
a common feature: the function x* [1 — P {x*)] has a 
maximum which corresponds to the critical value of the 
initial applied stress. All threshold distributions hav- 
ing this property are therefore expected to lead to the 
same universality class as the three studied here. If the 
threshold distribution does not have this property we 
may not observe a phase transition at all. For exam- 
ple, consider a fiber bundle model with P{x) — 1 — 1/x, 
X > 1. Here x* [1 — P {x*)] = 1 and the evolution 
of the fiber bundle is given by the recursion relation 
Ut+i = Ut/(T which implies that there is no dynamics 
at all for a — I and an exponential decay to complete 
failure, Ut = (c) S for cr > 1. There are no critical 
phenomena and therefore no phase transition. However 
this gen eral conclusion may not be tr ue for finite-sized 
bundles (jMcCartnev and Smithl . Il983l ). 

Thus the ELS fiber bundles (for different fiber thresh- 
old distributions) show phase transition with a well de- 
fined order parameter which shows similar power law 
variation on the way the critical point is approached. 
For all the cases, discussed here, the susceptibility and 
relaxation time diverge following similar power laws and 
the failure processes show similar critical slowing at the 
critical point. This suggests strongly that the critical be- 
havior is u niversal, which we now prove through general 
arguments ([Hemmer et al ]. l2006fl . 

When an iteration is close to the fixed point, we have 
for the deviation 



V+AUt 



to lowest order in AUt- This guarantees an exponential 
relaxation to the fixed point, AUt oc e"*/"^, with param- 
eter 



1 / In 



U 



*2 



ap{a/U*) 



(74) 



Criticality is determined by the extremum condition (Eq. 
IH), which by the relation (Eq. [5]) takes the form 

= apia/Uc). 

Thus T = oo at criticality. To study the relaxation at 
criticality we must expand Eq. (|73p to second order in 
AUt since to first order we simply get the useless equation 
AUt+i — AUt- To second order we obtain 

AUt+i = AUt ~ CA[/2, 

with a positive constant C. This is satisfied by 

1 



AUt = 



Ct 




FIG. 8 Branching ratio as a functio n of applied stress f or 
three different p (WeibuU index) [from (jMoreno et aZ.Ll2000l )]. 



1. The values a — (3 = 9 = ^ can be shown to be 
consequences of the parabolic maximum of the load curve 
at criticality. Thus all threshold distributions for which 
the macroscopic strength function has a single parabolic 
maximum, is in this universality class. 

It is clear that at the critical stress value CTc, ELS fiber 
bundles show phase transition from partially broken state 
to completely bro ken state. Wh a t is the orde r of this 
phase transition? IZapperi et al\ i 19971 Il999at ) consid- 
ered the fraction of unbroken fibers as the order param- 
eter and as it has a discontinuity at the critical stress 
value, they suggested, after a mean-field analysis, that 
it can be seen as a first-order phase transition similar 
to spinodal instability (Mgncttc, 1994). The additional 



Hence in general the dominating critical behavior for the 
approach to the fixed point is a power law with rj = 



reason for identifying the trans ition at a = ar. a s a first- 
order spinodal point had been (jKun et all |2000( ) that in 
the presence of short-range interactions (as in LLS, see 
Section IV), the transition becomes discontinuous and 
first-order like. It is indeed hard to identify continu- 
ously changing order parameter there. We, however, be- 
lieve the transition in ELS to be second-order. Chrono- 
l ogically, a little late r, a new parameter was identified 
( Moreno et al\ . [2OOOI ): the branching ratio {(), which is 
defined as the probability of triggering further breaking 
given an individual failure. The branching ratio con- 
tinuously approaches (Fig. [8]) the value 1 at the criti- 
cal stress ((Tc) starting from value (for very small a). 
Also it shows a power law variation: 1 — C oc {a^ — cr)^, 
with (3 — 1/2. Therefore 1 — C acts as the order pa- 
rameter showing a continuous transition at the criti- 
cal point, signalin g a second-order phase transition . As 
mentioned earlier Pradhan and Chakrabartil (|200lh and 
iPradhan et d\ (|2002l ) considered the difference between 
the fraction of unbroken fibers at any a and at CTc, as the 
order parameter (O) : it shows a similar continuous varia- 
tion with the applied stress: O oc {ac — fj}^, with (3 — 1/2. 
Apart from this, the susceptibility and relaxation time di- 
verge at the critical point following power laws having 
universal exponent v alues ( Bhattacharvva et al\ . 120031 : 
iPradhan et al. I. l2002fl . One may therefore conclude that 
at the critical point the ELS fiber bundles show a second- 
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FIG. 9 Number of relaxation steps tf{a) for a fiber bundle 
with a uniform threshold distribution (Eq. [!}. Here tJc = 
0.25. The figure is based on 1000 samples, each with A'^ — 10^ 
fibers. 



order phase transition with robust critical behavior as 
discussed here. 

Finally we compare the ELS fiber bundle model stud- 
ied here with the mean-field Ising model. Though the 
order parameter exponent (equal to ^) of this model is 
identical to that of the mean-field Ising model the two 
models are not in the same universality class. The sus- 
ceptibility in these models diverge with critical exponents 
i and 1 respectively on approaching the critical point. 
The dynamical critical exponents are not the same ei- 
ther: in this fiber bundle model the surviving fraction 
of fibers under the critical applied stress decays toward 
its stable fixed point as t~^, whereas the magnetization 
of the mean-field Ising model at the critical temperature 
decays to zero as 



4. Relaxation behavior and critical amplitude ratio 

When an external load F is applied to a fiber bun- 
dle, the iterative failure process continues until all fibers 
fail, or an equilibrium situation with a nonzero bundle 
strength is reached. Since the number of fibers is finite, 
the number of sl eeps, tf, in this sequen ti al pro cess is fi- 
nite. Following iPradhan and Hemmeil ( 20071 ). we now 
determine how tf depends upon the applied stress a. 

The state of the bundle can be characterized as pre- 
critical or post-critical depending upon the stress value 
relative to the critical stress CTc ~ Fc/N, above which the 
bundle collapses completely. The function tf{<T) that we 



now focus on, exhibits critical divergence when the criti- 
cal point is approached from either side. As an example, 
we show in Fig. [5] the tf(a) obtained by simulation for a 
uniform threshold distribution. 

We study the stepwise failure process in the bundle, 
when a fixed external load F = Na is applied. Let Nt 
be the number of intact fibers at step no. t, with A^o = 
A''. We want to determine how Nt decreases until the 
degradation process stops. When iV is a large number, 



we recall the basic recursion (Eq. ^ to formulate the 
breaking dynamics: 

Ut+i = 1 - Picj/Ut), (75) 

where Ut = Nt/N is considered as a continuous variable. 

(a) Post-critical relaxation 

We study first the post-critical situation, a > ac, with 
positive values of e = ct — CTc, and start with the sim- 
plest one, uniform threshold distribution (Eq. [T]) with 
the critical point at Xc = 1/2, CTc = 1/4. Then the basic 
recursion relation (Eq. I75p takes the form 



(76) 



This nonlinear iteration can be transformed into a linear 
relation. We introduce first Ut — \ — Vt\f^-, into Eq. (|76p . 
with a result 



Vt+\ - Vt 



1 + VtVt+i 
Then we put yt — tan Vt , which gives 



(77) 



tanwt+i — tanuf 



= tan(?;t+i - vt). (78) 



1 -|- tanwf+i tanuf 
Hence we get Vt+i — vt — tan^^(2\/e), with solution 

vt^VQ + t isai-^{2^fe). (79) 
In the original variable the solution reads 



Ut 



e tan tan { 



1 

1,2 



-)+ t tan-i(2V^) (80) 



= I - tan (-tan~^(l/2\/J) + i tan"^(2x/J)) , (81) 

where Uq = 1 has been used. 

Eq. (fZS)) shows that when Ut obtains a value in the 
interval (0,(t), the next iteration gives complete bundle 
failure. Taking Ut — <J as the penultimate value gives a 
lower bound, i'^, for the number of iterations, while using 
C/t = in Eq. ([5T|) gives an upper bound t"^. Adding 
unity for the final iteration, Eq. (|8ip gives the bounds 



t}{a) = 1 + 



2tan-i(l/2Ve) 
tan-i(2Vi) ' 



and 



1 



tan-i((i - e)/Vg) + tan-\l/2^) 
tan"i(2Vi) 



(82) 



(83) 



Fig. \TO\k shows that these bounds nicely embrace the 
simulation results. 

Note that both the upper and the lower bound behave 
as e~2 for small e. A rough approximation near the crit- 
ical point is 



tf{a) « K+((T - acY 



(84) 
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FIG. 10 Simulation results with post-critical stress for (A) the 
uniform threshold distribution (Eq. [TJ, and (B) the WeibuU 
distribution (Eq. |2]) with index 5. The graphs are based 
on 10000 samples with N = 10^ fibers in each bundle. Open 
circles represent simulation data and dashed lines are the the- 
oretical estimates Eq. ((HH), Eq. in (A) and Eq. (|92} in 
(B). 



with = 7r/2. 

Due to the inherent simplicity, uniform distribution is 
somewhat easy to analyze. Therefore we now discuss 
how to handle other distributions. Let us start with a 
WeibuU distribution (Eq. ^ with index 5. The critical 
parameters for this case are Xc = 5~^/^ = 0.72478 and 

= (5e)-i/^ = 0.5933994. 

The interesting values of the external stress are close 
to (Tc, because for large super-critical stresses the bundle 
breaks down almost immediately. For a slightly above 
(Tr the iteration function 



f{Ut) = 1 - Pi<j/Ut) - e-('^/^')' 



(85) 



takes the form sketched in Fig. [TT] 

The iteration function is almost tangent to the reflec- 
tion line Ut+i = Ut and a long channel of width propor- 
tional to e appears. The dominating number of iterations 
occur within this channel (see Fig. [TT]). The channel wall 




0.7 0.75 0.8 0.85 0.9 0.95 1 

FIG. 11 The iteration function f(U) for the WeibuU distri- 
bution (Eq. [2} with index 5. Here a — 0.6, slightly greater 
than the critical value — 0.5933994. 



formed by the iteration function is almost parabolic and 
is well approximated by a second-order expression 

Ut+i ^Ue + (Ut - Uc) + a{Ut - Ucf + b{a, - a). (86) 

Here Uc = e^^l^ is the fixed point, Ut+\ = Ut, of the 
iteration ai a = ac- With u = (U — Uc) /b and e = a — (Tc 
Eq. ([86]) takes the form 



H+l 



Ut 



-Aui 



(87) 



with A = ah. In the channel u changes very slowly, so 
we may treat the difference equation differential 
equation: 



du , 9 
— = -Au^ 
dt 



with solution 

tVAe = - tan"^ {u^jA/ej + constant . (89) 

Thus 

te-ts = {Aey^ ^i&Yr^{us\/ A/f) - tan"^(uev^A/e)| 

(90) 

is the number of iterations in the channel, starting with 
Us, ending with Ug. This treatment is general and can 
be applied to any threshold distribution near criticality. 
Although the vast majority of the iterations occur in the 
channel, there are a few iterations at the entrance and 
at the exit of the channel that may require attention in 
special cases. The situation is similar to type I intermit- 
tency in dynamical systems, but in our case the channel 
is traversed merely once. 

For the WeibuU distribution the expansion (Eq. 
has the precise form 



Ut = e~("/^)' ~e-i/5 + (?7-C/c) 

-|ei/5((7-(7,)2_5i/5(a-ae) 



(91) 
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where Uc = e'^l^, a = |el/^ b = S^/^ and A = |(5e)i/^ 
For completeness we must also consider the number of 
iteration to reach the entrance to the channel. It is not 
meaningful to use the quadratic approximation (Eq. I9ip 
where it is not monotonously increasing, i.e. for U > 
U,n = Uc + l/(2a) = f e~^/^ ~ 0.98. Thus we take Us = 
Um as the entrance to the channel, and add one extra 
iteration to arrive from [/q = 1 to the channel entrance. 
(Numerical evidence for this extra step: For cr = CTc the 
iteration (Eq. [55]) starts as follows: Uq — 1.00, Ui — 
0.93, U2 — 0.90, while using the quadratic function with 
Uq = Um = 0.98 as the initial value, we get after one 
step Ui = 0.90, approximately the same value that the 
exact iteration reaches after two steps.) With C/g = we 
obtain from Eq. (1^1) . in the WeibuU case, the estimate 



if 



1 + (Ae)-i/2 |tan~i(e-i/5y]47^/55) 

+ tan-i(e-i/5V^/&)}, (92) 



with A = |(5e)i/5 and h = S^/^. 

Near the critical point Eq. ([5^ has the asymptotic 
form 



tf^n{Ae)-^'^ = ^i+{a-a,)-^'\ 



(93) 



with K+ = 7r(2/5)^/^(5e)"i/i". The critical index is the 
same as for the uniform threshold distribution. The the- 
oretical estimates give an excellent representation of the 
simulation data (see Fig. [TOB). 

(h) Pre-critical relaxation 

We now assume the external stress to be pre-critical, 
(7 < (Tc, and introduce the positive parameter e — ac — <J 
to characterize the deviation from the critical point. 
Starting with uniform threshold distribution and intro- 
ducing Ut = \ + -Jej zt and a ~ j — e into Eq. ((76l) , one 
gets 



2y/I : 



Zt+1 ~ Zt 



1 - Zt+1 Zt 

In this case we put zt — tanh wt , which gives 



(94) 



_ tanh wt+i - tanh Wt 

2Ve = 7 — -, 7 — ; = tanh(u)f+i - wt). (95) 

1 — tanhwt+i tanhwt 

Thus Wt+i — Wt = tanh^"^ (2y^), i.e. 

wt^wa+t tanh"^(2Ve). (96) 
Starting with t/g = 1, we obtain zq = 2y^ and hence 

Wt = (1 +t) tanh"^(2Ve). (97) 
This corresponds to 



tanh {(1 +t) tanh ^(2^^)} 



(98) 
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FIG. 12 Simulation results with pre-critical stress for (A) the 
uniform threshold distribution (Eq. [!}, and (B) the WeibuU 
distribution ([2]) with index 5. The graphs are based on 10000 
samples with A'^ = 10^ fibers in each bundle. Open circles rep- 
resent simulation data and the dotted lines are the theoretical 
estimates, Eq. (fTOTj) in (A) and Eq. (|111)I112|I in (B). 



in the original variable. 

Apparently Ut reaches a fixed point U* = \ + ^ after 
an infinite number of iterations. However, for a bun- 
dle with finite number of fibers, only a finite number 
of steps is needed for the iteratio n to arrive at a fixed 
point N* which is approxiin ately ( Pradhan et al\ . [20021 : 
iPradhan and Hemmerl . |2007^ 



^* = f (1 - ^/T^) + ^ (i + (1 - 4^^) 



-1/2 



(99) 



AT* _ 1 
- ~ 2 



TT7 2 + 



-1/2 



(100) 



2 ^ ' 2 

As a consequence, we can use 

1 

4/V 

as the final value in Eq. ((98)l . Consequently we obtain 
the following estimate for the number of iterations to 
reach this value: 

coth"^ {! + (! + 2Ve)/4Af£} 



t/(^) = -1 



tanh"^(2Ve) 



(101) 



14 



Fig. [I2K shows that the simulation data are well approx- 
imated by the analytic formula (Eq. llOip . 

For very large N Eq. (jlOip is approximated by 



ln(AO ,/2 



-1/2 



(102) 



with K_ = ln(7V)/4. The critical behavior is again char- 
acterized by a square root divergence. 

Again we use the WeibuU distribution (Eq. ^ as an 
example threshold distribution. In principle, the itera- 
tion, 



Ut+i = 1 - Picr/Ut), 



(103) 



will reach a fixed point U* after infinite many steps. The 
deviation from the fixed point, Ut — U*, will decrease 
exponentially near the fixed point: 

Ut-U* (X e~^l\ (104) 

with 

T = \l\Yv{U*'^a-^ lv{olU*)] . (105) 

For the WeibuU threshold distribution with index — 5 

p(a/t/*) = 5(a/f/*)4 exp (-(a/f/*)^) = SctVC/*^, 

(106) 

and thus 



T = l/ln(C/*V5^T^)- 



(107) 



If we allow ourselves to use the exponential formula (Eq. 
I104p all the way from C/q = 1, we obtain 



(108) 



For a finite number N of fibers the iteration will stop 
after a finite number of steps. It is a reasonable supposi- 
tion to assume that the iteration stops when Nt — N* is of 
the order 1. This corresponds to take the left-hand side 
of Eq. (jlOSp equal to 1/A^. The corresponding number 
of iterations is then given by 



tf^TlniN{l-U*)) 



in general, and 



In (A^(l - [/*)) 
ln(C/*V5o-^) 



(109) 



(110) 



in the WeibuU case. Solving the WeibuU iteration U* — 
exp(— (tj/?/*)^) with respect to a and inserting into Eq. 
(IllOp . we obtain 



tf = - 



ln{iV(l- [/*)} 



ln{5(-lnC/*)} 

(7 = C/*(-lnC/*)i/^ 



(111) 
(112) 



These two equations represent the function t(a) on pa- 
rameter form, with U* running from Uc — to 
U* = 1. 



For U* = Uc = e-i/^ Eq. (fTTTI) shows that tj is infi- 
nite, as it should be. To investigate the critical neighbor- 
hood we put U* = Uc{l + AC/) with AC/ small, to obtain 
to lowest order 



tf = 



In(iV) 
5AC/ 



<Jc-a = ^a(AC/)2 



(113) 
(114) 



The combination of Eq. (|113p and Eq. pi4p gives, once 
more, the square root divergence 



now with the magnitude 

= 10-i/2(5e)-i/i"ln(iV). 



(115) 



(116) 



Simulation results for the pre-critical WeibuU distribu- 
tion are shown in Fig. I12B . which shows good agreement 
with the analytic solution (Eqs. 11111 1112"|l . 

For a general threshold distribution the divergence and 
its amplitude are most easily deduced by expanding both 
the load curve a — x[l — P{x)] and the characteristic time 
T around the critical threshold x^- To lowest contributing 
order in a;c — a; we find 



1 



[2p{Xc) + Xcp'iXc)] + (Xc - xY 



and 



Xcp{Xc) 



2p{xc) + xIp'{xc) 



-{Xc 



(117) 



(118) 



Inserting for (x^. — x) from the equation above, and using 
(fTTO]) . we find 



</ = K-(o-c - 0-) 



-^^-1/2 



(119) 



with 



K_ = Xcp{x,)[Ap{x,) + 2xcp'(xc)]-^'^ In(iV). (120) 

To show how the magnitude of the amplitude k_ de- 
pends on the form of the threshold distribution, we con- 
sider a WeibuU distribution 



P{x) = 1 - e^"/"^)" 



(121) 



with varying coefficient p, and constant average strength. 
With a — r(l -f- 1/p) the average strength {x) equals 
unity, and the width takes the value 

w = {{x^) - {x)^Y = (r(l + 2/p)/T\l + l/p) -ly. 

(122) 

Here F is the Gamma function. Using the power series 
expansion T{1 + z) ^ I - 0.577z + 0.9892^ + ... we see 
how the width decreases with increasing p: 



1.52 



w ' 



(123) 
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For the Weibull distribution (Eq. I121[) we use Eq. 
(|120p to calculate the amplitude with the result 

K_ = {T{l + l/p)/2p)h{pe)-^'^nn{N) ^ (2p)-^ ln(7V), 

(124) 

the last expression for large p. Comparison between Eq. 
(|123p and Eq. (|124p shows that for narrow distributions 



(125) 



That narrow distributions give small amplitudes could 
be expected: Many fibers with strengths of almost the 
same magnitude will tend to break simultaneously, hence 
the relaxation process goes quicker. 




Cl 1/2 
FIG. 13 The stress-strain curve for an ELS bundle having 
uniform fiber strength distribution with a low cutoff Cl = 0.2. 



(c) Universality of critical amplitude ratio 



As function of the initial stress a the number of relax- 



-1/2 



at the 



ation steps, tf{a), shows a divergence \a — ac 
critical point, both on the pre-critical and post-critical 
side. This is a generic result, valid for a general proba- 
bility distribution of the individual fiber strength thresh- 
olds. On the post-critical side tf{a) is independent of the 
system size N for large N . On the pre-critical side there 
is, however, a weak (logarithmic) iV-dependence, as wit- 
nessed by Eqs. (46), (47) and (55). Note that the critical 
amplitude ratio takes the same value = ln(A'^)/27r 

for the uniform and the Weibull distributions. This shows 
the universal nature of the critical amplitude ratio, inde- 
pendent of the threshold distributio n. N ote the difference 
with normal critical phenomena |Aharonvl . 119761 ) due to 
the appearance of the In(A^) in this amplitude ratio here. 



5. Non-linear stress-strain behavior 

Fiber bundle model captures c orrectly the non - linear 
elastic behavio r in ELS mode ( Pradhan et I2OO2I: 
ISornettd . [l989l) . In case of strain controlled loading, us- 
ing the theory of extreme order statistics, it has been 
shown (Sornettej, 1989 ) that ELS bundles shows non- 
linear stress-strain behavior after an initial linear part up 
to which no fiber fails. Similar non-linear behavior is seen 
in the force controlled loading case as well. Moreover, 
from the recursive failure dynamics, the amount of stress 
drop at the breaking point can be calculated exactly 
( Pradhan et al\ . |2002| ) . To demonstrate the scenario we 
consider an ELS bundle with uniform fiber strength dis- 
tribution, having a low cutoff Cl, such that for stresses 
below the low cutoff, none of the fibers fail. Hence, un- 
til failure of any of the fibers, the bundle shows linear 
elastic behavior. As soon as the fibers start to fail, the 
stress-strain relationship becomes non-linear. This non- 
linearity can be easily calculated in the ELS model, using 
Eq. ([5]) for the failure dynamics of the model. 

Fibers are here assumed to be elastic, each having unit 
force constant, with their breaking strengths (thresholds) 



distributed uniformly within the interval [Cl , 1] : 
I \ / 0, 0<x<Cl 



i-Cl ■ 



For an applied stress a < Cl none of the fibers break, 
though they are elongated by an amount e = x = a. The 
dynamics of breaking starts when applied stress a be- 
comes greater than Cl- For a > Cl, the basic recursion 
relation (Eq. [8]) takes the form: 



Ut+i = 



1 



1-Cl 

which has stable fixed points: 
1 



2(1 - Cl) 



1 



a 

1 



1/2- 



(127) 



(128) 



The model now has a critical point (Tc — 1/[4(1 — Cl)] 
beyond which the bundle fails completely. At each fixed 
point, there will be an equilibrium elongation e{cr) and a 
corresponding stress S ~ U*e{cr) develops in the system 
(bundle). From Eq. p27)) . one gets (for a > Cl) 



1 



1-Cl^ 



(7 



(129) 



Also, from the force balance condition, at each fixed- 
point e{a) — X* . Therefore, the stress-strain relation for 
the ELS model finally becomes: 



S^l eil 



0, 



Cl 



< cr < Cl 
Cl < cr < ac 
a > ffc 



(130) 



The stress-strain relation in a ELS bundle is shown in 
Fig. [131 where the initial linear region has unit slope (the 
force constant of each fiber) . This Hooke's region for the 
stress S continues up to the strain value e = Cl, until 
which no one of the fibers breaks. After this, nonlinearity 
appears due to the failure of a few of the fibers and the 
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consequent decrease of U*{a). It finally drops to zero 
discontinuously by an amount x*U*{ac) = 1/[4(1 — Cl)] 
at the breaking point cr = CTc or e = = 1/2 for the 
bundle. It may be noted that in this model the internal 
stress X* is universally equal to 1/2, independent of Cl 
at the failure point a = ac- 



6. Effect of a low cutoff: Instant failure situation 
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A low cutoff in the fiber threshold distribution ex- 
cludes the presence of very weak fibers in a bundle. The 
weaker fibers mainly reduces the strength of a bundle. 
But in practice we always try to build stronger and 
stronger materials (ropes, cables etc.) from the fibrous 
elements. Therefore this situation (exclusion of weaker 
fiber s) is very reahstic. In t his se ction we discuss the ef- 
fect ( Pradhan and Hansenl . |2005| ) of a low cutoff on the 
failure properties of ELS bundles. 

We follow the weakest fiber break ing approach 
(|Danielsl . [l94llHemmer and Hansenl . Il992f ): The applied 
load is tuned in such a way that only the weakest fiber 
(among the intact fibers) will fail after each step of load- 
ing. We first find the extreme condition when the whole 
bundle fails instantly after the first fiber ruptures. As the 
strength thresholds of N fibers arc uniformly distributed 
between Cl and 1 (Eq. I126p . the weakest fiber fails at a 
stress Cl (for large N). After this single fiber failure, the 
load will be redistributed within intact fibers resulting a 
global stress Xf = NCl/{N — 1). Now, the number of 
intact fibers having strength threshold below Xf is 



NP{xf) = N / p{y)dy 



/■Xf 

Jcl 



Njxf - Cl) 

(1-Cl) ■ 



(131) 



Stress redistribution can break at least another fiber 
if NP{xf) > 1 and this 'second' failure will trigger an- 
other failure and so on. Thus the successive breaking of 
fibers cannot be stopped till the complete collapse of the 
bundle. Clearly, there cannot be any fixed point (criti- 
cal point) for such 'instant failure' situation. Putting the 
value of Xf we get 



Ni 



NCl 
N-1 



Cl) 



(1-Cl) 



> 1; 



which gives 



Cl > 



{2N- 1)' 



(132) 



(133) 



In the large N limit the above condition can be written 
as Cl > 1/2. Therefore, the condition to get a fixed 
point in the failure process is Cl < 1/2. 

We can also calculate how many steps are required to 
attain the final catastrophic failure for Cl < 1/2. Let 
us assume that we have to increase the external load n 
times before the final failure. At each step of such load 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

c, 



FIG. 14 The number of steps of load increase (till final fail- 
ure) is plotted against Cl for a ELS model having 50000 
fibers. The dotted line represents the analytic form (eqn. 
I136|) . triangles are the simulated data for a strictly uniform 
strength distribution, and the circles represent the data (av- 
erages are taken for 5000 samples) for a uniform on average 
distribution. 



increment only one fiber fails. Then after n step the fol- 
lowing condition should be fulfilled to have a catastrophic 
failure: 



rx,ll + l/{N~n)] 

N I p{y)dy > 1. (134) 



where 

Xi = Cl 
The solution gives 



^(1 - Cl) 
N 



Cr 



2 1-Cl 



(135) 



(136) 



The above equation suggests that at Cl = 1/2, n = 0. 
But in reality we have to put the external load once to 
break the weakest fiber of the bundle. Therefore, n = 1 
for Cl > 1/2 (Fig. \T^. To check the validity of the above 
calculation we take 'strictly uniform' and uniform on av- 
erage distributions of fiber strength. In our 'strictly uni- 
form' distribution the strength of the k-th fiber (among 
N fibers) is Cl + (1 - CL)k/N. We can see in Fig. M 
that the 'strictly uniform distribution' exactly obeys the 
analytic formula (Eq. I136P but the uniform on average 
distribution shows slight disagreement which comes from 
the fiuctuation in the distribution function for a finite 
system size. This fiuctuation will disappear in the limit 
iV — !■ oo where we expect perfect agreement. 



B. Fluctuations 

If the contribution to breakdown phenomena in mate- 
rials science by statistical physics were to be expressed 
in one word, that word would have to be "fiuctuations" . 
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In the context of fiber bundles, this concept refers to the 
effects of the fibers each having properties that are statis- 
tically distributed around some mean, which cannot be 
reproduced by substituting the fiber bundle by an equiv- 
alent one where each fiber is identical to all the others. 

Intuitively, it is not difficult to accept that fluctuations 
must play an important role in the breakdown properties 
of fiber bundles — or in fracture and breakdown phe- 
nomena in general. A plane ride in turbulent weather 
compared to one in smooth weather is a reminder of this. 

Closely connected to the question of fluctuations i s 
that of phase transitions and criticality ( Stanley! . Il987( ). 
Leaving the fiber bundles for a moment, consider a fiuid 
whose temperature is slowly raised. At a well-defined 
temperature determined by the surrounding pressure, the 
fluid starts to boil. Each gas bubble that rise to the sur- 
face is due a fluctuation being larger than a well-deflned 
size for which the bubble grows rather than shrinks away. 
At a particular pressure the character of the boiling 
changes character. There is no longer any size that deter- 
mined whether a nascent bubble grows or shrinks. There 
are bubbles of all sizes. At this particular point, the sys- 
tem is critical and undergoes a second order phase tran- 
sitions. The boiling process at other pressures signals a 
flrst order transition. 

A brittle material under stress develops microcracks. 
These appear where the material is weak or where the 
local stress field is high. As the stress increases, more 
and more microcracks accumulate until either one or a 
few microcracks go unstable and grow to macroscopic di- 
mensions causing failure. The spatial fluctuations of the 
local material properties cause the appearance of micro- 
cracks. Their subsequent growth accentuate these initial 
fluctuations, but in a highly complex manner due to in- 
teractions between the growing cracks. There are simi- 
larities between this scenario and a first order transition 
( Zapperi et aZ.I . [l997l ). On the other hand, stable mode 
I crack growt h as studied experimentally by Mal 0y and 
Schmittbuhl (|Mal0v and Schmittbuhil . ll997l . l200ll ) shows 
all the signs of the advancing crack front showing a dy- 
namics compatible with being at a critical point. 

We now turn to the global load sharing fiber bundle 
model in light of the preceeding remarks. 



1. Burst distribution for continuous load increase 

When a fiber ruptures somewhere, the stress on the 
intact fibers increases. This may in turn trigger further 
fiber failures, which can produce bursts (avalanches) that 
either lead to a stable situation or to breakdown of the 
whole bundle. A burst is usually defined as the amount or 
number (A) of simultaneous fiber failure during loading. 
One may study the distribution D(A) of the bursts ap- 
pear during the entire failure process until the complete 
breakdown of the bundle. 

The property of the fiber bundle model of interest in 
the present context, is the fluctuation driven hurst dis- 
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FIG. 15 The solid curve indicates the total force F{x) as 
a function of x. However, when our control parameter is F 
rather than x, the system will follow the dotted line, Eq. ([138} . 
The bursts are the horizontal parts of Fph{x). Here = 100. 



tribution. In order to define this property, we again con- 
sider a finite bundle containing N elastic fibers whose 
strength thresholds are picked randomly from a proba- 
bility density p{x). Let Xk be the ordered sequence of 
failure thresholds: xi < X2 < ■■■ ^ xn- Then the the 
external load or force F on the bundle (Eq. [3]) at the 
point where fcth fiber is about to fail can be written as: 



Fk = iN + l-k)xk, 



(137) 



where elastic constant of the fibers is set equal to unity as 
before. Note that the sequence of external loads Fk is not 
monotonously increasing. This may be readily seen from 
Eq. (|137p ; the total load is the product of a monotonously 
increasing fluctuating quantity Xk and a monotonously 
decreasing quantity {N + 1 — k). Suppose now that our 
control parameter is the total load F, and that /c— 1 fibers 
have broken. In order to be in this situation, F > Fk > 
Fj for all j < k. The latter inequality ensures that the 
situation we are studying is not unstable. We increase 
F until it reaches Fk, at which fiber k breaks. If now 
Fk+i < Fk, then fiber k + 1 will also break without the 
external load F being further increased. The same may 
be true for Fk+2 and so on until the (A: -I- A — l)th bond 
breaks. Thus, Fk+j < Fk for j < A. If now Fk+A > Fk, 
the burst of breaking bonds then stops at this point, and 
we have experienced a burst event of size A. 

The total force F expressed as a function of elongation 
X, is shown in Fig. 1151 When the control parameter is 
elongation x, the solid curve is followed. However, when 
the force F is the control parameter, the broken line, 
given by 

Fph = LMF Fix) , (138) 
where LFM designates the least monotonic function, 
(a) Generic case 
It was shown in lHemmer and HansenI (|l992t) that the 



average number of burst events of size A per fiber. 
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D{A)/N, follows a power law of the form 

D{A)/N = CA-« 
in the limit N ^ oo. Here, 



(139) 



2 (140) 

is the universal burst exponent. The value (Eq. I140p is, 
under very mild assumptions, independent of the thresh- 
old distribution P{x): the probability density needs to 
have a quadratic maximum somewhere in the interval 
a^min < a; < a;max- We demonstrate this in Fig. [TBI The 
prefactor C in Eq. p39p is given by 

C = ^ ^ ^^^^^ 

where Xc is the solution of the equation 

XcP{Xc) = 1 - P{xc) , (142) 

and is the value of x for which the characteristics 
ha s a maximum. Eqs . ([13911 to (|142p were derived 
in [Hemmer and Hanserj (|1992( ) using combinatorial ar- 
guments. We will, however, in the following, take 
an alternative route based on a mapping between the 
global load sharing model and a Brown ian process 
( Hansen and Hemmed . 119941: ISornettd . I1992D . 3efore we 
explain this mapping we quote, for later comparison, 
the pertinent results of the Hemmer-Hansen analysis 
([Hemmer and Hansenl . [1992[ ): The probability <i>(A,a::) 
that a burst event at elongation x will have the size A is 



$(A,a;) = 



where 



A 



A-l 



m(a;) 
A! 1 - m{x) 



m{x) = 1 — 



[1 — m(x)]e 



xp{x) 



va{x) — 1 



1 - P{x) 



(143) 



(144) 



Note in particular that by Eq. (|142p m{xc) — 0. Let us 
now assume that we do not load the fiber bundle until 
complete collapse, i.e., until x = x^, but stop at a value 
Xs < Xc- We may then ask for D{A, Xs)/N, the expected 
number of burst events of size A during the breakdown 



process that occurs between x = and x = 
given by the integral 

^^^/;>(a:)dx$(A,x) = 

^ dx p{x)j^, [[1 - m(x)]e'"(^)l 



/2ir 



This is 



,(145) 



where on the right-hand side the Stirling approximation 
A! « v^A^+i/^e"^ for large A has been used. The 
integrand in Eq. (|145p is strongly peaked near x — Xc- 
We therefore expand it to second order in y = Xc — x to 
find 



D{A,x,) _ A-^^^ 



N 



/2ir 



p{xc)m'{xc) 



r dy ye-^'i^^)"y'^n 



where we have extended the upper integration limit to 
CO. We may do this integral to get 



P>{A,x,) 
N 



where C is defined by Eq. (|14ip . 
We may write (jl47l) in scaling form, 

Ri^^ ^ A-«G(A, X,) = A-«G {A^ixc - x.)) , 

(148) 

where 



G{y) = C e-"'(^=)'2''/2 



(149) 



In particular G{y) tends to the constant C for y — ^ 
0. Two universal critical exponents appear, ^ — 5/2, 
Eq. pID[) . and 



(150) 



It is, thus, in the above sense, that the fracture pro- 
cess of the fiber bundle approaches a critical point at 
total breakdown: The distribution of burst events fol- 
lows a power law with an upper cutoff that diverges as 
the bundle approaches total failure. 

In [S ornettc (1993) and later on in full detail in 
[Hansen and Hemmen ( 1994[ ). the burst distribution (Eq. 
I139p was derived from the assumption that F^. may be 
directly interpreted as a biased random walk. The pre- 
cise nature of this random walk is elucidated below. It is 
a peculiar asymmetric walk with variable step length. In 
the limit iV — > oo and continuous time variable k /N — > t 
and Ak/N —^ St > 0, this random walk may be mapped 
onto a continuous Brownian process. Such Brown i an pro - 
cesses have been stu d ied by [Phoenix and Tavlod ( 1973[ ). 
[Daniels and Skvrmd ([1985() and [Panielsl ([l989t) m con- 
nection with the distribution of the strength S of fiber 
bundles. We will in the following derive Eq. (|145p by 
means of a biased random-walk model with variable step 
length. We find this an interesting example of univer- 
sality in statistical physics: The asymptotic behavior of 
one model is found by using a different model with the 
same asymptotic behavior as the first one, but which is 
simpler to solve before the continuum limit is taken. 

Under increasing load the variation of the force per 
fiber, / — F/N, will consist of a systematic nonfluctuat- 
ing part, given by the average load-elongation character- 
istics, with a small fluctuation of order superim- 
posed. 

The precise value of the force fluctuation depends upon 
whether one studies the force f{x) at given elongation, 
or the force fk at which flber number k breaks. Let us for 
both quantities calculate the variance of /, aj, starting 
with the constant-fc ensemble. 

The force per fiber when the fcth fiber is about to break 
is, hy Fk = iN+l- k)xk, 



(146) 



fk = [l- P{xk)]Xk 



(151) 
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FIG. 16 The burst distribution D{A)/N for the uniform 
distribution (A) and the Weibull distribution with index 5 
(B). The dotted hues represent the power law with exponent 
^ = 5/2. Both figures are based on 20000 samples of bundles 
each with iV = lO" fibers. 



where Xk is the elongation when the fcth fiber breaks, and 
we have defined Xk by 



P{xk) 



N + 1 



(152) 



For large N Xk is essentially the average value of Xk- For 
a fixed k the variance of fk is by Eq. (|151[) given by the 
variance of Xk- 



a}ik) = [l-Pixk)rcTl , 



(153) 



and we seek therefore the probability ip{x)dx that the 
fcth threshold in the ordered threshold sequence lies in 
the interval {x,x + dx). This probability equals 



fix) dx = 

(k-mN-ky. n^f~'\^ p{x)r-'pix)dx 



(154) 



For large k and N, and using Eq. (|152p . this is close to 
the Gaussian distribution 



(p{x) dx - (^2^p^^P^j^^l^p(_^^j^ 



1/2 



^-Npixty(x~Xkr/2P(xk)[l-P{xk)] 

This gives the variance of Xk and thus of fk, 
P(a;fe)[l-P(xfc)]^ 



(155) 



aUk) 



Npixkf 



(156) 



Let us now compare with the force fluctuation at con- 
stant elongation. The force per fiber is the following func- 
tion of elongation x, 



N 



f{x)^N-^Y.^Q{U-x) , 



(157) 



where ti is the breakdown threshold for the iih fiber, and 
Q{t) is the Heavyside function. This gives immediately 
the average force 

(/), -x[l-P(x)] , (158) 
i.e., the characteristics, as well as the variance 

x^P{x)[l ~ P{x)] 



aj{x) = 



N 



(159) 



Although the two types of force fluctuations have dif- 
ferent variances, in both cases a oc 1/^/N. 

The nonmonoticities of the force / within the fluctu- 
ation zone produce bursts. Since the fluctuations are so 
small for large N, one can treat the burst events locally. 

We now consider the force sequence Fk as a stochastic 
process. Since we seek the asymptotic burst distribution, 
we are interested in the behavior after many steps of the 
process. It is convenient, however, to start with the one- 
step process. 

Let us determine the probability distribution of the 
force increase AF = Pfc+i — Fk between two consecutive 
bursts, the first one taking place at elongation Xk with 
Fk = {N-k + l)xk. 

Since AF = (N — k){xk+i — Xk) — Xk, we have 



AF > -Xk . 



(160) 



The probability to find the fc + I'th threshold in 
(xfc+i, Xfe+i + dxk+i) for given Xk, 



N-k-2 



-p{xk+i)dxk+i , (161) 



gives directly by using the connection Xk+i = Xk + 
{AF{xk)+Xk)/{N-k), the probability density p{AF; Xk) 
of AF 

piAF;xk) = 



N-k-l 



N-k 



Pi 



[l-PCxfc)]"-"- 
, AF+x, . 



(162) 
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For large N — k this simplifies to 
piAF;xk) 

^ jl-P(..)exp PV";tS?'' J 



for AF > -Xk, 
for AF < -Xk- 



(163) 

This one-dimensional random walk is asymmetric in more 
than one way. First of all, it has nonzero bias 



{AF){xk) = 



1 - Pjxk) - XkpjXk) 

P{xk) 



(164) 



In addition the probability distribution around this av- 
erage is very asymmetric. 

The variance is easily determined, 



l~P{xk) 



(165) 



The Brownian-motion limit of a one-dimensional ran- 
dom walk is completely determined by the first and sec- 
ond moments of the single-step probability distribution. 
The results just obtained enables us therefore to select 
an "ordinary" biased random walk with constant step 
length a which has the same Brownian motion limit as 
the burst process. 

We imagine having a one-dimensional random walk 
along the z-axis with a constant bias. Each step is of 
length a. Let the probability to take a step in the nega- 
tive z direction be q, and let p be the probability to take 
a step in the positive z direction. The walk is biased 
when p is different from q. The probability distribution 
of the position zi after one step has the average 



(zi) = a{p-q) , 



and variance 



af = Apqa? 



(166) 



(167) 



Elimination of a yields 



p-q ^ (fi) 
2VM cri 



Since p + g = 1, the bias parameters are determined. 
After k steps a Gaussian distribution. 



[zk—Z(:)~ka{p — q)]'^/^pqa'^k 



^Snpqa^k 



(169) 



is approached when k increases. 

The two processes will have the same asymptotic be- 
havior when we make the identification 



p-q ^ (AF) 



XkP(xk) 
l-P{Xk) 



mixk) , (170) 



where m{x) is defined by Eq. (|144p . When the bias is 
small, both p and q are close to 1/2, and we have to 
lowest order 



p= i[l + m(a;)] , 
g = i[l - mix)] 



(171) 



We have now made the promised mapping between the 
fiber bundle problem and a random walk with a constant 
bias. A constant bias may be used since bursts can be 
treated locally. 

The next step is to calculate the burst distribution for 
such a biased random walk. Since this biased random 
walk by construction has the same asymptotic behavior 
in the limit A'^ — > cxd as the original fiber bundle problem, 
the two burst distributions will asymptotically be the 
same. 

In terms of the biased random walk, a burst event of 
size A at "time" k may be defined as follows: (i) Zk+i < 
Zk for < i < A and Zk+A > Zk- (ii) Furthermore, to 
ensure that we are not counting burst event inside other 
burst events, the condition Zk > Zj for fc > j is necessary. 

The first condition i s in fac t a spe cial case of the "Gam- 
bler's ruin" problem ( Felleij . Il966l ). A gambler plays a 
series of independent games against a bank with infinite 
resources. In each game, the gambler either looses or 
wins one Euro, and the probability that the bank wins 
is p — {1 + B)/2, while the probability that the gambler 
wins is g = (1 — B)/2. If the gambler starts out with a 
capital of z Rs, the probability that she is ruined after 
precisely A games is 



/AN Z f A 

^(z,A) = -U_^ 

^ \ 2 2 



p(A-.)/2^(A+.)/2 (^72) 



The probability that condition i is fulfilled for a biased 
random walk burst of size A is then 



H^ = l,A) = ^7II|(l-i?r/^ (173) 



where we have assumed that A >> 1. 

The probability that a biased r andom walker returns at 
least once to the origin is (iFellej . [19661 l-\p-q\ = l-B. 



(168) The probability that condition (ii), namely that Zj > Zk 
for all j < fc is fulfilled, is then simply 1 — {1 — B) — B, 
and we have that the probability for having a burst of 
size A happening at "time" k is 



RW 



(A,B) = ii3^(^ = l,A) 



A-3/2 



2tt 



B e 



-B^A/2 



(174) 



where we in addition have assumed that i? << 1. 

Returning to the fiber bundle model, the bias B — 
m{x). When x is close to Xc, we have B — m'{xc)y where 
y = {xc — x). Thus, the probability to have a burst of 
size A between y and y + dy is 

$w(A, m'{xc)y)p{xc)dy 
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FIG. 17 The macroscopic bundle strength {F}{x) for the dis- 
tribution (Eq. [177}, with xo — 2xr, and for ao = ^ (upper 
curve), I (middle curve), and | (lower curve). The broken 
part of the ao = 2/3-curve is unstable and the macroscopic 
bundle strength will follow the solid line. 



^-3/2 



(175) 



Thus, the cumulative burst distribution up the elonga- 
tion a;, is 



$w(A, m'{xc)y)p{xc)dy 



(176) 



Comparing this expression to Eq. (I146p . we see that 
they are identical. This completes the derivation of the 
asymptotic burst distribution via the mapping between 
the random-walk problem and the burst process. 

(b) Special cases 



The burst distribution given in Eq. (|145p is valid when 
the threshold distribution has a parabolic maximum in- 
side the interval of the thresholds. We now consider 
threshold distributions that do not reach their maximum 
at the boundaries of the interval ( Kloster et al 
Model examples of such threshold distributions are 



P{x) 



for X < xq 

1 — [1 -f (a; — xo)/xr]~°'° for x > xq 



(177) 



Here ao and xq are positive parameters, and refer- 
ence quantity which we for simplicity put equal to unity 
in the following. These distributions are all characterized 
by diverging moments. When ao ^ Ij even the first mo- 
ment — the mean — as well as all other moments diverge. 
This class of threshold distributions are rich enough to 
exhibit several qualitatively different burst distributions. 

The corresponding macroscopic bundle strength per 
fiber is 



{F){a 
N 



X for a; < ^0 

{i+x-xo^o for x>Xo 



(178) 



In Fig. [17] the corresponding macroscopic force curves 
{F){x) are sketched. We note that when ao — > 1, the 
plateau in Eq. (|178p becomes infinitely wide. 

The distribution of burst sizes is given by Eq. (|145p . 
In the present case the function m(x) takes the form 



m{x) 



xp(x) 



aQX 



1 - P{x) 1 + x-xq' 
A simple special case is xo = 1, corresponding to 

p{x) = ao x^"""^ for a; > 1, 

since then the function (|179p is independent of /: 
m{x) = ao. 

This gives at once 

^(A) 1-ao A^-i 



(179) 



N 



ao A! 

1-ao 



ao V 27r 



aoe 

A-t \aoe^-^°] 



1 A 



A 



(180) 



In other cases it is advantageous to change integration 
variable in Eq. (jl45p from x to m: 



J (ao — m)"''^^(l — m)TO^^ (me^^™) dm. 

(181) 

The asymptotics for large A, beyond the A~i depen- 
dence of the prefactor, is determined by the A-dependent 
factor in the integrand. The maximum of me^~™ is unity, 
obtained for m = 1, and the asymptotics depends cru- 
cially on whether m — I falls outside the range of integra- 
tion, or inside (including the border). If the maximum 
falls inside the range of integration the D{A) cx A~^ 
dependence remains. A special case of this is ao = 1, 
for which the maximum of the integrand is located at 
the integration limit and the macroscopic force has a 
"quadratic" maximum at infinity. Another special case 
is aoXo — 1 (and ao < 1), for which again the standard 
asymptotics A~2 is valid. In this instance the macro- 
scopic force has a quadratic minimum at x = xo (see 
Fig. [TTjfor ao = 1/2), and critical behavior arises just as 
well from a minimum as from a maximum. 

In the remaining cases, in which m = 1 is not within 
the range of integration in Eq. (|18ip . the burst distribu- 
tion is always a power law with an exponential cut-off, 



^(A) 
N 



A-« . 



(182) 



Here ^ and A depend on the parameter values xo and ao, 
however. This is easy to understand. Since 



dm{x) ao(l — Xo) 



dx 



(1 



Xo) 



(183) 
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Parameters 



Asymptotics 



0<xo<l,ao<l A-2-"o{ae 
< 2:0 < l,ao = 1 A"3 
a;o = 1, ao < 1 



1 < xo < a^^ 
1 < xo = a^^ 



A-i(Qoe'-"«)'^ 



A-^ 



1 < Qq ^ < a:o A 



aoxoe 
A-i 

-A/Ao 



TABLE I Asymptotic behavior of the burst distribution for 
strong threshold distributions in the ELS model. 



we see that m{x) is a monotonically decreasing function 
for Xo > 1, so that the maximum of me^~™ is obtained 
at the lower limit x — xq, where m = aoXQ. The asymp- 
totics 



D{A) (X A 2 (aoXoe 



l-aoxa 



y 



(184) 



follows. 

This is true merely for aoxo < 1, however. For ao^o > 
1 the macroscopic force {F)(x) decreases near x = Xq so 
that a macroscopic burst takes place at a force xq per 
fiber, and stabilization is obtained at a larger elongation 
xi (Fig. [T7|) . The subsequent bursts have an asymptotics 



Z?(A)cx A-i (a(/i)ei-"(-i))'^, 



(185) 



determined by the neighborhood of x = xi. For to < 1, 
the maximum of me^~'^ is obtained at x = c«, leading 
to the asymptotics 



D{A) cx A" 



(aoe^--)^. 



(186) 



reflecting the power-law behavior of the integrand at in- 
finity. 

The results are summarized in Table H] Note that the 
xq = 1 result (I180p cannot be obtained by putting xq = 1 
in Eq. p84p since in ()181|) the order of the limits A ^ 00 
and Xo — > 1 is crucial. 

(c) Crossover behavior 

When all the bursts are recorded for the entire failure 
process, we have seen that the burst distribution D(A) 
follows the asymptotic power law D oc A"^/^. If we 
just sample bursts that occur near the breakdown point, 
a different behavior is seen. As an illustration we con- 
sider the uniform threshold distribution, and compare the 
complete burst distribution with what one gets when one 
samples merely burst from breaking fibers in the thresh- 
old interval (0.9xc,Xc). Fig. [T8l shows clearly that in the 
latter case a different power law is seen. 

This observation may be of practical importance, as it 
gives a criterion for th e imminence of catastrophic failure 
(jPradhan efoll . I2005D . This proposal has so far not been 
tested experimentally. However, it is enticing to note the 
recent observation by Kawamura of a crossover behavior 



D(A) 
D(l) 




FIG. 18 The distribution of bursts for thresholds uniformly 
distributed in an interval (a;o, Xc), with a;o = and with a;o = 
0.9xc- The figure is based on 50 000 samples, each containing 
= lO'' fibers. The exponents of the distributions are ^ = 
5/2 and ^ = 3/2 respectively. 



in the magnitude distribution of earthquakes before large 
earthquakes appear (Kawamura, 200 6). We return to this 
result in Section V.C and Fig. [HH 

We introduce the following notation in Eq. (|145p . 

D{A) _ A^-^e-^ 
N 



A! 



Jp ''p(x)r(x)[l — r(x)] ^ exp [Ar(x)] dx , 



where 



r(x) 



1 



xp(x) 

Q{x) Q{x) dx 



^ :^[a;Q(x) 



(187) 



(188) 



and Q{x) — p{x) dx. We note that r[x) vanishes 
at the point Xc- If we have a situation in which the 
weakest fiber has its threshold xo just a little below the 
critical value Xc, the contribution to the integral in the 
expression (Eq. I187P for the burst distribution will come 
from a small neighborhood of Xc- Since r(x) vanishes at 
Xc, it is small here, and we may in this narrow interval 
approximate the A-dependent factors in Eq. p87p as 
follows 



(1 - r)^ e^'' = exp [A(ln(l - r) + r)] 
exp[-A(rV2 + ©(r^))] w exp [-Ar(x)V2] 



We also have 



r(x) « r'(xc)(x - Xc) . 



(189) 



(190) 



Inserting everything into Eq. (|187p , we obtain to dom- 
inating order 



D(A) 

N 



A! 



p(xc) r'(xc)(x - Xc) 
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FIG. 19 The distribution of bursts for the uniform threshold 
distribution for a single fiber bundle with A'^ = 10^ fibers. 
Results with xo = 0, i.e., when all bursts are recorded, are 
shown as squares and data for bursts near the critical point 
{xo = 0.9xc) are shown by circles. 



A 



A-2 „-A 



|r'(xe)| A! 



A 



A-2 „-A 



A! 



k'(Xc)| 



with 



A, 



(191) 



(192) 



By use of the Stirling approximation A! ~ 
A^e"'^V27rA, the burst distribution (Eq. [TgT|) may be 
written as 



CA-5/2 fi_g-A/A, 



^(A) 

N 

with a nonzero constant 

C=(2^)-l/2p(a;,)/|r'(Xe)| 



(193) 



(194) 



We can see from Eq. (|193p that there is a crossover at a 
burst length around Ac, 



g(A) ^ I A-3/2 for A < Ac 
N 1 A-5/2 for A > Ac 



(195) 



We have thus shown the existence of a crossover from 
the generic asymptotic behavior D cx A~^/^ to the power 
law D (X A"'^/^ near criticality, i.e., near global break- 
down. The crossover is a universal phenomenon, inde- 
pendent of the threshold distribution p{x). 

The simulation results we have shown so far are based 
on averaging over a large number of samples. For appli- 
cations it is important that crossover signal can be seen 




FIG. 20 A fuse model of size 100 x 100. Each bond is a 
fuse with a burn-out threshold t drawn from a probability 
distribution p{t). 



also in a single sample. We show in Fig. [19] that equally 
clear crossover behavior is seen in a single fiber bundle 
when is large enough. Also, as a practical tool one 
must sample finite intervals {xi, Xf) during the fracture 
process. The crossover will be observed when the interva l 
is close to the failure point ( Pradhan et aLl . [20051 . l200(t l. 

The ELS fiber bundle model is a simple model in that 
it is analytically tractable. A step up in complexity 
from t he ELS fiber bundle mode l, is the random fuse 
model (jHerrmann and Ro"iju3 . ll990l ). While resisting most 



analytical treatments, this model retains computational 
tractability. The fuse model consists of a lattice in which 
each bond is a fuse, i.e., an ohmic resistor as long as the 
electric current it carries is below a threshold value. If 
the threshold is exceeded, the fuse burns out irreversibly. 
The threshold t of each bond is drawn from an uncor- 
related distribution p{t). The lattice is placed between 
electrical bus bars and an increasing current is passed 
through it. The lattice is a two-dimensional square one 
placed at 45° with regards to the bus bars, and the Kirch- 
hoff equations are solved numerically at each node assum- 
ing that all fuses have the same resistance. We show the 
model in Fig. [5D1 The ELS fiber bundle model may be 
interpr eted as a mean-field version of the random fuse 
model ( Zapperi et all . Il997t ). Hence, the random fuse 
model may be used as a testing ground for results (see 
Table ini found with the ELS fiber bundle model to ex- 
plore their robustness when other effect not present in 
the fiber bundle model enter. 

To test the crossover phenomenon in a more complex 
situation than for ELS fiber bundle model, we consider 
the random fuse model ( Pradhan et aLl , [20061 ). When one 
records all the bursts in the random fuse model, the dis- 
tribution follows a power law D{A) (x A"'' with ^ ~ 3, 
which is consistent with the value reported in recent stud- 
ies. We show the histogram in Fig. With a system 
size of 100 X 100, 2097 fuses blow on the average before 
catastrophic failure sets in. When measuring the burst 
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FIG. 21 The burst distribution based on 300 samples random 
fuse lattices of size 100 x 100. The threshold t are uniformly 
distributed on the unit interval. On the average, catastrophic 
failure sets in after 2097 fuses have blown. The circles de- 
note the burst distribution measured throughout the entire 
breakdown process. The squares denote the burst distribu- 
tion based on bursts appearing after the first 1000 fuses have 
blown. The triangles denote the burst distribution after 2090 
fuses have blown. The two straight lines indicate power laws 
with exponents ^ = 3 and ^ ~ 2, respectively. 



distribution only after the first 2090 fuses have blown, a 
different power law is found, this time with S. — 2. After 
1000 blown fuses, on the other hand, ^ remains the same 
as for the histogram recording the entire failure process 

(Fig. mi). 

In Fig. [22l we show the power dissipation E in the net- 
work as a function of the number of blown fuses and as 
a function of the total current. The dissipation is given 
as the product of the voltage drop across the network V 
times the total current that flows through it. The break- 
down process starts by following the lower curve, and 
follows the upper curve returning to the origin. It is in- 
teresting to note the linearity of the unstable branch of 
this curve. In Fig[23l we record the avalanche distribu- 
tion for power dissipation, Dd{A). 

Recording, as before, the avalanche distribution 
throughout the entire process as well as recording only 
close to the point at which the system catastrophically 
fails, result in two power laws, with exponents ^ = 2.7 
and ^ = 1.9, respectively. It is interesting to note that 
in this case there is not a difference of unity between the 
two exponents. The power dissipation in the fuse model 
corresponds to the stored elastic energy in a network of 
elastic elements. Hence, the power dissipation avalanche 
histogram would in the mechanical system correspond 
to the released energy. Such a mechanical system could 
se rve as a simple model for eart hquakes. 

iDivakaran and Duttal ( 2007al) studied the critical be- 
havior of a bundle of fibers under global load sharing 
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FIG. 22 Power dissipation i? as a function of the number of 
broken bonds (upper) and as a function of the total current 
/ flowing in the fuse model (lower). 
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FIG. 23 The power dissipation avalanche histogram Dd{A) 
for the fuse model. The slopes of the two straight lines are 
—2.7 and —1.9, respectively. The circles show the histogram 
of avalanches recorded after 1000 fuses has blown, whereas 
the squares show the histogram recorded after 2090 fuses have 
blown. This is close to catastrophic failure. 
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FIG. 24 Mixed uniform distribution. 



scheme with threshold strength chosen randomly from a 
distribution which is uniform, but discontinuous. The 
form of the distribution is 



p{x) 



1 



I - {X2- Xi) 







1- {X2- Xi) 



< X < xi 

Xi < X < X2 

X2<x<l (196) 



where X2 — xi is the gap in the threshold distribution as 
shown in Fig. [M] 

Here, a fraction / of the fibers belong to the weaker 
section (0 < x < xi) and the remaining to the stronger 
section {x2 < a; < 1). The condition of imiformity of the 
distribution demands 



Xi = 



f 



1-/ 



(1-2^2) 



(197) 



so that fixing xi and / immediately settles the value of 

X2. 

To study the dynamics of this model, once again the 
recursive equation approach was used. Importantly, the 
redistributed stress must cross X2 for the complete failure 
of the bundle to take place. When the external load is 
such that the redistributed stress at an instant t, i.e., x(t) 
is greater than X2, the fixed point solution has the form 



U* 



1 



2(1 - {X2-Xl)) 



1 



1 



(J 



so that the critical stress Ur is 



4[l-(2;2-xi)] 



(198) 



(199) 



and the redistributed stress at the critic al point is found 
to be 1/2 as in the uniform distribution ( Pradhan et al\ . 
[2005 ). This immediately restricts the value of X2 to be 
less than 1/2 and therefore, X2 = 0.5 is defined as the 
critical distribution in this model. On the other hand, 
uniformity condition, Equation (|197p . puts another re- 
striction, namely xi < f. Interestingly, the critical stress 
is a function of the gap X2 — xi and reduces to one fourth 



when the gap goes to zero. The exponents related to the 
order parameter and susceptibility stick to their mean 
field values. However, the existence of a forbidden region 
shows a prominent signature in the avalanche size distri- 
bution of the mixed model. The expression for the total 
avalanche size distribution D{A) in this model includes 
two terms, one due to the contribution from thresholds 
between to xi and the other from the stronger sec- 
tion of fibers. Hence, the total avalanche size -D(A) is 
D{A)/N = Di{A) + Da (A) where 



A! 



1-xl+xi lo '^^( 



1—3:2+2:1 —2x > 



f X eXp(-r; f ) 

1 — X2+XI—X l — X2-\-Xl—X' 



(200) 



and 



D2{A) = 



A 



A-l 



A! I-X2 + X1 



0.5 



d.(i^) 

X 



--^exp(-— ^) 

1 — X 1 — X 



(201) 



The leading behav ior of -Di(A) is given by 

(jPivakaran and Duttal . [2007al) 



i?i(A) = A- 



where 



-5/2g(l- 



l~ X2 



(202) 



(203) 



which clearly indicates a rapid fall of the contribution 
of weaker fibers. On the other hand, D ^jA) resembles 
the im minent failure behavior studied by Pradhan et al\ 
(I2OO5D . where the avalanche size exponent shows a 
crossover from 5/2 to 3/2 as X2 — > 0.5. For the mixed 
model, the total avalanche size distribution D{A) shows 
a nonuniversal behavior for small A values, though even- 
tually there is a crossover to the universal mean field 
value. The most fascinating observation is the follow- 
ing: though the gap in the distribution is always present, 
nonuniversality is only prominent in the limit X2 — > 0.5. 
Divakaran and Dutta showed that this nonuniversal be- 
havior stems from the avalanche of fibers in the weak 
section and only in the vicinity of the critical distri- 
bution, the contribution of Di{A) overcomes I?2(A). 
Otherwise, the faster fall of Di{A) and large value of 
D2{A) together force the avalanche size exponent to be 
5/2. It is to be noted that the nonuniversal behavior 
is most prominent at a critical distribution where the 
avalanche size exponent crosses over to 3/2 in the asymp- 
totic limit. The typical behavior of D{A) is shown in 
Fig- ESI for two different distributions highlighting the in- 
crease in nonuniversal region as X2 approaches 0.5. How- 
ever, for many discontinuities in the threshold distribu- 
tions, avalanche size dis tribution shows a nonuniver sal, 
non power-law behavior (jPivakaran and D^,[2003) for 
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FIG. 25 Total avalanche size distribution D{A), Di{A) and 
D2(A) obtained by numerical integration of Eqs. (I200|l and 
(|201[) . Figure (a) corresponds to xi = 0.08, a::2 = 0.28 and 
/ = 0.1 and Figure (b) to xi = 0.25, X2 = 0.42 and / = 0.3. 
As X2 0.5, the nonuniversal region increases in the small A 
region, whereas A oc A~^^^ for large A. The dotte d line has 
a slope of -5/2. From IPivakaran and Duttal (|2007al ). 



small-size avalanches, although the large avalanches still 
exhibit similar crossover behavior as we discuss here. 

Divakaran and Dutta also looked at the model where 
fibers fro m two difi'erent Weibu U distributions are mixed 
fPivaka ran and Duttal . l2007bl ). Though an interesting 
variation of the critical stress with the mixing parame- 
ter was obtained using a probabilistic method introduced 
by iMoreno et al. I (|2000D . there is no deviation in the 
avalan che size exponent. In a recent paper. [Hidalgo et al\ 
(I2008D studied the infinite gap limit of the discontinuity 
model. Here, they considered a fraction ainf of the fibers 
having infinite threshold strength mixed with fibers hav- 
ing threshold chosen from a distribution p{x). They ob- 
served a critical fraction Ur such that for > an, the 



inf 

avalanche size exponent switches from the well known 
mean field exponent ^ = 5/2 to a lower value ^ = 9/4. 
It was also showed that such a behavior is observed for 
those distributions where the macroscopic constitutive 
behavior has a maxima and a point of inflexion. It is 
also claimed that below a critical gap, the Hidalgo et al. 
model redu ces to the discontinu ity model of Divakaran 
and Dutta. iKun and Nagvl (|2008t) studied the global load 
sharing fiber bundle model in a wedge-shaped geome- 
try. That is, the fibers are connected to two rigid blocks 
placed at an angle with respect to each other. The fibers 
are loaded by rotating the blocks with respect to each 
other, resulting in a linear loading gradient on them. In 
the limit of a threshold distribution tending towards zero 
width, the fibers break in an orderly fashion according to 
the load, and hence, position in the wedge. As the width 
is increased, a process zone — i.e., a zone where some 
fibers fail whereas others stay intact — develops. When 



the width of the threshold distribution is wide enough, 
the process zone spans the entire bundle. In this limit 
a burst size exponent ^ = 5/2 is recovered. However, 
with a narrower distribution so that a well-defined pro- 
cess zone smaller than the size of the bundle develops, 
the burst exponent ^ = 2.0 is found. 



2. Burst distribution for discrete load increase 

When the bundle is stretched continuously from zero, 
fluctuation plays crucial role and the ge neric result is a 
power law ijHemmer and Hansenl . 1 19921) D(A) oc A"?, 
for large A, with ^ = 5/2. However, experiments may 
be performed in a different manner, where the load is in- 
creased in finite steps of size 5 . The valu e of the exponent 
incre ases ( Hemmer and Pradhanl l2007t iPradhan et al\ . 
l2002f ) then to 3: D(A) oc A'^. The basic reason for 
the difference in the power laws is that increasing the ex- 
ternal load in steps reduces the fluctuations in the force. 
The derivation (see Section III.B.l) of the asymptotic size 
distribution r'(A) oc A~^/^ of avalanches, corresponding 
to stretching by infinitesimal step s, shows the imp ortance 
of force fluctuations ( Hemmer a nd Hansenl . [1993 ) . An ef- 
fective reduction of the fluctuations requires that the size 
i5 of the load increase is large enough so that a consider- 
able number of fibers break in each step. 

Here is an a n alytic derivation, following 
iHemmer and Pradhanl (I2007D . how to calculate the 
burst distribution in such situation. For the uniform 
distribution of thresholds (Eq. [T|), the load curve is 
parabolic. 



{F)=N x{l- x). 



(204) 



so that the expected critical load equals Fc = N/A. With 
a sufficiently large 5 we may use the macroscopic load 
equation (Eq. I204p to determine the number of fibers 
broken in each step. The load values are m5, with m 
taking the values m = 0, 1, 2, ... , N/A5 for the uniform 
threshold distribution. By Eq. (|204p the threshold value 
corresponding to the load m5 is 



i (l - v/1 - 4mJ/Ar) . 



(205) 



The expected number of fibers broken when the load is 
increased from m5 to (m -I- 1)5 is close to 



A = Ndxr,Jdm = 6/y/l-AmS/N. 



(206) 



Here the minimum number of A is 6, obtained in the first 
load increase. The integral over all m from to N/AS 
yields a total number N/2 of broken fibers, as expected, 
since the remaining one-half of the fibers burst in one 
final avalanche. 

The number of avalanches of size between A and A -t- 
c?A, -D(A) dA, is given by the corresponding interval of 
the counting variable m: D{A) dA — dm. Since 



dA 

dm 



2P 
N 



{l~Am6/N)-^/'^ 



NS 



(207) 
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FIG. 26 Avalanche size distribution for the uniform threshold 
distribution (Eq. [T| when the load is increased in steps of 
S = 10 and ^ = 50 (upper curve). The dotted lines show the 
theoretical asymptotics (Eq. I208|) for 5 = 10 and S = 50. The 
figure is based on 10000 samples with A'^ — 10^ fibers in the 
bundle. 




FIG. 27 Avalanche size distribution for the WeibuU distribu- 
tion ((2)) with index 5. Open circles represent simulation data, 
dashed lines are analytic expressions (Eq. I213|) and dotted 
line is the asymptotic power law with exponent —3. The load 
is increased in steps of (5 = 20. The figure is based on 10000 
samples of bundles with N — 10^ fibers. 



we obtain the 

( Hemmer and Pradhanl . 



fo ll owing dist ri bution 

120071 : iPradhan et al. I, I2002D 



of avalanche sizes: 



D{A) 



dA ^ 



(A > 5). 



(208) 



For consistency, one may estimate the total number of 
bursts by integrating ^(A) from A = S to oo, with the 
result N/AS, as expected. 

Fig. shows that the theoretical power law (Eq. I208P 
fits the simulation results perfectly for sufficiently large 
A. The simulation records also a few bursts of magnitude 
less than 6 because there is a nonzero probability to have 
bundles with considerably fewer fibers than the average 
in a threshold interval. However, these events will be of 
no importance for the asymptotic power law in the size 
distribution. 

In order to see whether the asymptotic exponent value 



^ = 3 is general, simulations for another threshold dis- 
tribution has been performed, the WeibuU distribution 
(Eq. [2) with index 5, which confirms similar asymptotic 
behavior (Fig. \T7\ . 

For a general threshold distribution P{x) a load in- 
terval S and a threshold interval are connected via the 
load equation (F) = N x {1 - P{x)). Since d{F)/dx = 
N[l — P{x) — xp{x)\, an increase 5 in the load corresponds 
to an interval 



dx — 



5 



N[l - P{x) - xp{x)] 



(209) 



of fiber thresholds. The expected number of fibers broken 
by this load increase is therefore 



A = N p(x) dx = 



pix) 



1 — P{x) — xp{x) 



(210) 



Note that this number diverges at the critical point, i.e. 
at the maximum of the load curve, as expected. 

Following the similar method, as in case of uniform 
distri bution, we can determine ( Hemmer and Pradhanl 
l2007f ) the asymptotic distribution for large A: 



^(A) ~ C A" 



with a nonzero constant 
C = N6 



p{xc 



2p{Xc) + XcP'{Xc 



(211) 



(212) 



where we have used that at criticality 1 — P{xc) = 
XcP{xc)- Thus the asymptotic exponent value ^ = 3 is 
universal. 

For the WeibuU distribution considered in Fig. [27] we 
obtain 



D{A) = NS A- 



25x^ 



5x° 



and A = 



5^ ; 



1 



5x^ 
(213) 

This burst distribution must be given on parameter form, 
the elimination of x cannot be done explicitly. The criti- 
cal point is at a; = 5~^^^ and the asymptotics is given by 
Eq. [2ni with C = NS{625e)-^/^. 

If we let the load increase S shrink to zero, we must 
recover the asymptotic D{A) oc A^^/^ power law valid 
for continuous load increase. Thus, as function of i5, there 
must be a crossover from one behavior to the other. It 
is to be expected that for 5 < 1 the D(A) cx A"^/^ 
asymptotics is seen, and when S ^ 1 the D{A) cx A~'^ 
asymptotics is seen. 



3. Energy bursts in fiber bundle model 

So far we have discussed in detail the statistical 
distri bution of the size of ava l anches in fiber btm- 
dles ([Hemmer and Hansen . 1992|: Hemmer et al. 



iPradhan et all 2005 : Raischel et all 2006 ). Sometimes 
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the avalanches cause a sudden internal stress redistribu- 
tion in the material, and are accompanied by a rapid re- 
lease of mechanical energy. A useful experimental tech- 
nique to monitor the energy release is to measure the 
acoustic emissions (AE), the elastically radiat e d waves 



produced in the bursts (iDiodati et aL ^ 1997t 



1991 



- -- iFazzini 

iGarcimartm et all Il997l : iPetri et all 11994 IScott 



19911 ). Experimental observations suggest that AE sig- 
nals follow power law distributions. What is the origin 
of such power laws? Can we explain it through a general 
scheme of fluctuation guided breaking dynamics that has 
been demonstrated well in ELS fiber bundle model? 

We now determine the stati stics of the energies released 
( Pradhan and Hemmed . [20081 ) in fiber bundle avalanches. 
As the fibers obey Hooke's law, the energy stored in a 
single fiber at elongation x equals ^x'^, where we for sim- 
plicity have set the elasticity constant equal to unity. The 
individual thresholds assumed to be independent 

random variables with the same cumulative distribution 
fmiction P{x) and a corresponding density function p{x). 

(a) Energy statistics 

Let us characterize a burst by the number A of fibers 
that fail, and by the lowest threshold value x among the 
A failed fibers. The threshold value Xmax of the strongest 
fiber in the burst can be estimated to be 



A 



Np{x) ■ 



(214) 



since the expected number of fibers with thresholds in an 
interval Sx is given by the threshold distribution function 
as N p{x) Sx. The last term in (Eq. I214p is of the 
order 1/A^, so for a very large bundle the differences in 
threshold values among the failed fibers in one burst are 
negligible. Hence the energy released in a burst of size A 
that starts with a fiber with threshold x is given as 



(215) 



Following iHemmer and HansenI ( 19921 ) the expected 
number of bursts of size A, starting at a fiber with a 
threshold value in the interval {x,x + dx), is 



f{A,x) dx = N- 



1 — P{x) — xp{x) 



ni X 
X Xix)"^ e-'^^'-'^^ dx, 



where 



X{x) = 



x p(x) 

1 - p{xy 



(216) 



(217) 



The expected number of bursts with energies less than E 
is therefore 



\/2S/A 



dx. 



(218) 



with a corresponding energy density 
dG _ 

A 



dE 



^Y.i^EA)-'/' f{A,^/2E/A). (219) 



Explicitly, 



g{E)^NY,9A{E), 



(220) 



A 



with 



9AiE) 



2E A! 

sp{s) 



(1 - Pis) ~ spis)) 
sp{s) 



l-P(s) 



exp 



I -Pis) 



(221) 



Here s = ^/2E/A. With a critical threshold value Xc, it 
follows from (Eq. I215p that a burst energy E can only 
be obtained if A is sufficiently large, A > 2E/x'^. Thus 
the sum over n starts with A = 1 -I- [2E/x1], where [a] 
denotes the integer part of a. 

(b) High energy asymptotics 

Bursts with high energies correspond to bursts in which 
many fibers rupture. In this range we use Stirling's ap- 
proximation for the factorial A!, replace 1 -I- [2E/x'^] by 
2E/x'^, and replace the summation over A by an integra- 
tion. Thus 



N 



2^3/2^1/2 J /^3/2 



(1 - Pis) - spis)) 



spis) 



■ exp 



spis) 



dA. (222) 



l-P(s) ' V l-^(s). 
By changing integration variable from A to s we obtain 



9iE) 



N 



2P3/27ri/2 



spjs) 
1-Pis) 

N 

2P3/27ri/2 



(1 _ p(s) _ spi^s)) 



spis) 



exp 1 



ds 



1-Pis), 
(1 - P(s) - sp(s))e-^''(^) ds, 

(223) 



with 
his) 



1 - Pis) - spis) ^ 1 - Pis) 



1 - Pis) 



spis) 



(224) 



For large E the integral (Eq. I223P is dominated by 
the integration range near the minimum of his). At the 
upper limit s = Xc we have hixc) = 0, since 1 — Pixc) = 
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Xcp{xc). This is also a minimum of h(s), having quadratic 
form, 



his) 



xIp{xc 



{Xc - sY 



(225) 



Inserting these expressions into (Eq. I223P and integrat- 
ing, we obtain the following asymptotic expression. 



g{E) c N 



C 



E- 



where 



C 



X^p{Xcf 



4^1/2 [2p{Xc)+XcP'{Xc)]' 



(226) 



(227) 



In Fig. [28] we compare the theoretical formula with 
simulations for the uniform distribution (Eq. [T]), which 
corresponds to Xc = ^, and C — 2^^7r^^/^, and for the 
WeibuU distribution (Eq. ^ with index p — 2, which 
corresponds to Xc — 2^^/^ and C — 2~''^(27re)^^/^. 

The corresponding asymptotics (Eq. I226P are also ex- 
hibited in Fig. [551 For both threshold distributions 
the agreement between the theoretical asymptotics and 
the simulation results is very satisfactory. The expo- 
nent —5/2 in the energy burst distribution is clearly 
universal. Note that the asymptotic distribution of the 
burst magnitudes A i s gove rned by the same exponent 
([Hemmer and Hansenl . [l992f) . 

(c) Low- energy behavior 

The low-energy behavior of the burst distribution is 
by no means universal: g{E) may diverge, vanish or 
stay constant as E' ^ 0, depending on the nature of the 
threshold distribution. In Fig. [511 we exhibit simulation 
results for the low-energy part of g{E) for the uniform 
distribution and the WeibuU distributions of index 2 and 
index 5. 

We see that g{E) approaches a finite limit in the 
WeibuU p = 2 case, approaches zero for WeibuU p — 5 
and apparently diverges in the uniform case. All this 
is easily understood, since bursts with low energy pre- 
dominantly correspond to single fiber bursts (A = 1, i.e. 
E = x^ 12) and to fibers with low threshold values. The 
number of bursts with energy less than E therefore cor- 
responds to the number of bursts with x < V2E, which 
is close to N P{^/2E). This gives 
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FIG. 28 Simulation results for g{E) characterizing energy 
bursts in fiber bundles with (A) the uniform threshold dis- 
tribution (Eq. d} and (B) the WeibuU distribution (Eq. [2]) of 
index 2. The graphs are based on 1000 samples with A'^ — 10^ 
fibers in each bundle. Open circles represent simulation data, 
and dashed lines are the theoretical results (Eqs. 12261 - [2271) 
for the asymptotics. 



borderline case between divergence and vanishing of the 
low-energy density. The same lowest-order results can be 
obtained from the general expression (Eq. I220p . which 
also can provide more detailed low-energy expansions. 



g{E) ~ N 



when i? — > 0. 



(228) 



For the uniform distribution g{E) should therefore di- 
verge as {2E)-'^/'^ for E ^ Q. The simulation resuUs 
in Fig. [2S1 are consistent with this divergence. For the 
WeibuU Distribution of index 2, on the other hand, (Eq. 
27V when E 



12281) gives g{E) 2N when E ^ 0, a, value in agree- 
ment with simulation results in the figure. Note that for 
a WeibuU distribution of index p, the low-energy behav- 
ior is giE) cx £'(''-2)/2. Thus the WeibuU with p = 2 is a 



For high energies the energy density obeys a power 
law with exponent —5/2. This asymptotic behavior 
is universal, independent of the threshold distribution. 
A similar power law dependence is found in some ex- 
perimental observa t ions on acoustic emiss ion studies 
i|Garcimartfn et all Il997l: iPetri all Il994l) of loaded 
composite materials. In contrast, the low-energy behav- 
ior of g{E) depends crucially on the distribution of the 
breakdown thresholds in the bundle. g{E) may diverge, 
vanish or stay constant for E 0. 
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FIG. 29 Simulation results for the burst distribution g{E), in 
the low-energy regime, for the uniform threshold distribution 
(circles), the Weibull distribution with p = 2 (triangles) and 
WeibuU distribution with p = 5 (squares). The graphs are 
based on 1000 samples with A*" — 10® fibers in each bundle. 



Exponent for Value Comment 

Order parameter (a) 1/2 — 

Breakdown susceptibility (/3) 1/2 — 




FIG. 30 A fiber bundle with periodic boundary conditions. 
The externally applied force F is the control parameter. 

loading of the fibers — as can be seen in pulling on the 
hairs on one's arm. 



Relaxation time (f?) 



1/2 amplitude ratio — lnA''/27r 



Avalanche size distribution (^) 3 

5/2 



discrete load increase 
continuous load increase 



Energy burst distribution (^e) 5/2 in the asymptotic limit. 

For the low energy limit, 
distribution is non universal. 



TABLE II Exponents for order parameter (O), breakdown 
susceptibility (x), relaxation time (r), avalanche size distri- 
bution -D(A) and energy burst distribution g{E) in the ELS 
model. 



IV. LOCAL LOAD SHARING MODEL 

So far we have studied fiber bundles where the force 
once carried by a failing fiber is spread equally among 
all the surviving fibers. This may often be a very good 
approximation. However, intuitively it is natural that 
fibers being closer to a failing fiber feel more of an effect 
than fibers further away — an effect reminiscent of stress 
enhancement around cracks. We will in this section dis- 
cuss three classes of models where there are local effects 
in how the forces carried by failed fibers are distributed. 
We start with the most extreme where the forces are to- 
tally absorbed by the nearest surviving fibers. We then 
move on to models where the stress is distributed accord- 
ing to a power law in the distance from the failing fiber, 
and lastly a model where we assume the clamps the fibers 
are attached to are soft and therefore deform due to the 



A. Stress alleviation by nearest neighbors 

The extreme form for local load redistribution 
is that all extra stresses caused by a fiber fail- 
ure are ta ken up by the n e arest-neighbor s urviv- 
ing fib ers fPuxburv and L eathl. Il99 4': 'Harlow", Il985 : 



Harlow and Phoenix, 1981, 1991: ,Kuo and Phoenix 



1987 : Phoenix and Smi th|, I1983D . The simplest geom- 



etry is one-dimensional, so that the N fibers are or- 
dered linearly, with or without periodic boundary con- 
ditions. In this case precisely two fibers, one on each 
side, take up, and divide equally, the extra stress, see 
Fig. [30l When the strength thresholds take only two 
values, the b undle streng th distrib ution has beeii found 
analytically (iDuxbury and Leathl . Il994l : iHarlowl Il985l: 
iHarlow and Phoenix . 1991). 

At a total force F on the bundle the force on a fiber 
surrounded by n; previously failed fibers on the left-hand 
side, and on the right-hand side, is then 

^{l + \[ni+nr)) = f{2 + ni+nr) . (229) 



Here 



27V 



(230) 



is one-half the force-per-fiber, is a convenient variable to 

u se as the driving f o rce pa rameter. 

Zhang and Ding ( 1994 ) and [Hansen and Hemmeil 



(|1994l ) studied numerically the burst distribution in 
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FIG. 31 Burst distribution in local model as found numeri- 
cally for 4 000 000 samples with = 20000 fibers (crosses), 
and calculated from Eq. (|233|l (boxes). The straight line 
shows the power law and the broken curve the function 
exp(— A/Ao) with Ao = 1.1. Note the small value of Aq. 



the local load-sharing model. In Fig. [3T] we exhibit 
simulation results similar t o those first appearing in 
[Hansen and Hemmeil ( 19941 ) using threshold strengths 
randomly distributed on the unit interval. Again a 
power-law distribution seems to appear. However, the 
burst exponent ^ seems much larger than in the global 
lad-sharing model, Eq. I140[ 



(231) 



Thus the relative frequency of long (non-fatal) bursts is 
considerably reduced. 

It was concluded from the numerics that systems with 
local load-sharing are not in the universality class of fiber 
b undles with gl o bal lo ad redistribution. 

iKloster et all ( 19971 ) set out to analytically calculate 
the burst distribution in the local load-sharing model, 
finding the surprising result that there is no power law 
di stribution of bursts at all; it is exponential. 

iKloster et all ( 19971 ) did not find the burst distribution 
for general threshold distribution. Rather, they limited 
their study to the uniform threshold distribution in the 
force parameter /, given by 



Pif) 



f for < / < 1 
1 for / > 1. 



(232) 



Bursts in the local and the global models have different 
characters. In the local model a burst develops with one 
failure acting as the seed. If many neighboring fibers have 
failed, the load on the fibers on each side is high, and if 
they burst the load on the new neighbors will be even 
higher, etc. In this way a weak region in the bundle may 
be responsible for the failure of the whole bundle. For a 
large number N of fibers the probability of a weak region 
somewhere is higher. This hints in a qualitative way that 
the maximum load the bundle are able to carry does not 
increase proportional to iV, but slower than linear. 



The result of a c alculation based on combinatorics 
(|Kloster et a/.l . ll997t) was the burst distribution 



l/{A+2) 





Pf{N 



l,L2;/)[i-(ii + i2 + 2)/]d/. 



where S{1; /) is the probability that a selected region of 
I consecutive fibers have all failed whereas the two fibers 
at both boundaries are still intact at force parameter y, 
p{l,a; f)df is the probability that a force increase from / 
to f + df leads to a burst of length I and magnitude a. 
Pf(n,L;f) is the probability at force parameter / that 
among the n first fibers there is no fatal burst, and that 
the last L fibers of these have all failed. 

The load distribution rule (Eq. I229P implies that burst 
of size A does necessarily lead to a complete breakdown 
of the whole bundle if the external force is too high, i.e., 
if X exceeds a critical value Xmax- Since here a fiber can 
at most take a load of unity, we have 



fn 



1 



(234) 



Let us now attempt to find an simple estimate for the 
maximal force per fiber that the fiber bundle can toler- 
ate. In order to do that we assume that the fatal burst 
occurs in a region where no fibers have previously failed 
so that the burst has the same magnitude and length. 
We know that a single burst of length A = f^^ — 2 is 
fatal, Eq. (|234|) . so our criterion is simply 



1 



(235) 



If we take into account that the two fibers adjacent to 
the burst should hold, and ignore the rest of the bundle, 
the gap distribution would be 

l/(A+2) 

N-^DiA)^ J [l-(2 + A)/]MA,A;/)d/ 



_ 2p(A,A) 

~ A(A-|-l)(A+2)^+i • 



(236) 



With the abbreviation 



Ra 



P(A,A) 
(A-1)! 



we have 



D{A)/N 



2(A+2)! n 
A2(A+Tp(A+2)^+^ 



^87r(A+2) 



-A-2 



Ra, 



(237) 



using Stirling's formula. 
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Taking logarithms we have 
ln£)(A) -ln7V=:= -(A + 2) 



In B a 
A+2 



-(A + 2) 



using that 



hm R}/" = 1 



(238) 



(239) 



for i?A when A is large. 

The failure criterion (Eq. I235P then takes the form 



IniV: 



1 

7 



(240) 



Since / — F/2N we have the following estimate for the 
maximum force F that the fiber bundle can tolerate be- 
fore complete failure: 



F 



2N 



(241) 



Due to the assumption that the fatal burst occurs in a re- 
gion with no previously failed fibers, the numerical pref- 
actor is an overestimate. The size dependence 



F oc 



N 
hTiV' 



(242) 



shows that the maximum load the fiber bundle can carry 
does not increase proportionally to the number of fibers, 
but slower. This is to be expected since the probability 
of finding somewhere a stretch of weak fibers that start a 
fatal burst increases when the number of fibers increases. 

The N / In dependenc e agre es with a previous esti- 
mate by IZhang and Dind ( 19951 ) for a uniform thresh- 
old distribution. The bimodal distribution used in 
(jPuxburv and Leathl . Il994 iHarlow and Phoenixl . Il99lf ) 
also shows this behavior. 

The burst distribution Eq. (|238p is exponential. The 
probability of a single burst zipping through the fiber 
bundle grows with the system size N. This con- 
trasts strongly with the global load sharing model whose 
strength (maximum force it can sustain) grows linearly 
with N, and whose burst distribution follows a universal 
power law. If the latter behavior is reminiscent of a sec- 
ond order transition with a critical point, the local load 
sharing model behaves more as if moving towards a first 
order phase transition. 

Let us discuss briefly the effect of a low cutoff (in fiber 
strength distribution) on the failure properties of LLS 
model. As in case of ELS model (see Section III. A), we 
consider uniform fiber threshold distribution having a low 
cutoff Cl (Eq. I126p . We present a probabihstic argument 
to determine the upper limit of Cl, beyond which the 
whole bundle fails at once. Following the weakest fiber 
breaking approach the first fiber fails at an applied stress 
Cl (for large N). As we are using periodic boundary 
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FIG. 32 Numerical estimate of the upper bound of Cl in LLS 
model: For Cl > 0.5 the average step values go below 1.5, 
i.e., the bundle fails at one step in most of the realizations. 



conditions, the Uc nearest neighbors (ric is the coordina- 
tion number) bear the terminal stress of the failing fiber 
and their stress value rises to Xf = Cl(1 -I- l/?ic)- Now, 
the number of nearest neighbors (intact) having strength 
threshold below Xf is {nn)fau = ncP{xf) (see Eq. I13ip . 
Putting the value of P{xf) and Xf we finally get 



{nn)faii = 



(Cl 



I-ClY 



(243) 



If (nn)faii > 1, then at least another fiber fails and this 
is likely to trigger a cascade of failure events resulting 
complete collapse of the bundle. Therefore, to avoid the 
'instant failure' situation we must have {nn)faii < 1, 
from which we get the upper bound of C^: Cl < 5- As 
the above condition does not depend on the coordination 
number n, at any dimension the whole bundle is likely 
to collapse at once for Cl > 1/2. It should be men- 
tioned that LLS model should behave almost like ELS 
model at the limit of infinite dimensions and therefore 
the identical bound (of Cl) in both the cases is not sur - 
prising. A numerical study ( Pradhan and Hansenl . [20051 ) 
confirms (Fig. [5^ the above analytic argument in one 
dimension. When average step value goes below 1.5, one 
step failure is the dominating mode then. One can find 
the extreme limit of Cl when all the nearest neighbors 
fail after the weakest fiber breaks. Then the LLS bundle 
collapses instantly for sure. Setting (nn) fau = Uc one 
gets the condition Cl > nc/{l + Uc), where stress level 
of all the nearest neighbors crosses the upper cutoff 1 
of the strength distribution. Clearly such failure is very 
rapid (like a chain reaction) and does not depend on the 
shape of the strength distributions, except for the upper 
cutoff. Also as Tie increases (ELS limit), Cl for instant 
failure assumes the trivial value 1. Similar sudden fail- 
ure in FBM has been discussed bv iMoreno et al. (2001]) 
in the context of a 'one sided load transfer' model. 

The local load sharing (LLS) scheme introduces stress 
enhancement around the failed fiber, which accelerate 
damage evolution. Therefore, a few isolated cracks 
can drive the system toward complete failure through 
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FIG. 33 The strength - Cl is plotted against 1/iV for dif- 
ferent C'l values: 0.3 (square), 0.35 (circle), 0.4 (up triangle), 
0.45 (down triangle), 0.5 (star). All the straight lines ap- 
proach value as N ~* oo. 



growth and coalescence. The LLS model shows zero 
strength (for fiber threshold distributions starting from 
zero value) at the limit N oo, foll owing a logar i thmic 
dependence on the system size (N) (Gomez ej (7illl993l : 
|Pradhan and Chakrabarti, 2003a; Smith, 198^). Now for 
threshold distributions having a low cutoff (Cl), the ulti- 
mate strength of the bundle cannot be less than Cl ■ For 
such a uniform distribution (Eq. I126|) . numerical simu- 
lations show (Fig. [55]) that as Cl increases the quantity 
(strength-Ci) approaches zero following straight lines 
with 1/A'', but the slope gradually decreases -which sug- 
gests that the system size dependence of the strength 
grad ually becomes weaker. 

In [Hansen and Hemmerl ( 19941 ) a model interpolating 
between the global load sharing fiber bundle and a vari- 
ant of a local load shari ng model was i ntrod uced and 
studied. iKimI (|2004D and iPradhan et ai\ l|2005l) followed 
up th is work. In the model studied by Pradhan et al\ 
(I2005D a fraction g of the load a failing fiber carries would 
be distributed among its surviving neighbors and a frac- 
tion 1 — g among all surviving fibers. Hence, for g — 1 
the model would be purely local load sharing, whereas for 
g = it would be purely global load sharing. We show in 
Fig. [33] space time diagrams of the one-dimensional ver- 
sion of the model for different values of g. For increasing 
values of g, t h ere is incre asing localizatio n. 



Both iKimI (|2004t ) and iPradhan et oil (|2005l) found a 
phase transition when interpolating between the global 
load sharing model and the local load sharing model dis- 
cussed earlier in this section. In the one-dimensional 
model shown in Fig. [Ml the critical value of g, is 
Qc = 0.79 lb 0.01 for a flat threshold distribution 
(IPradhan et"aLl . l2005l ). 



B. Intermediate load-sharing models 

A crucial mechanism in brittle fracture is the stress en- 
hancement that occurs at crack tips. The stress field has 
a 1/ singularity, where r is the distance to the crack 




FIG. 34 Space time diagram of bre aking sequence s of th e 
interpolating fiber bundle model of IPradhan et al\ (2005') . 
Black lines represent broken fibers and white regions arc un- 
broken parts of the bundle. As the interpolation parameter 
g is increased, there is increasing localization of the failing 
fibers. 



tip, in this region. It is the interplay between fracture 
growth due to this singularity and due to weak spots in 
the mat erial that drive the develop ment of the fracture 
process (" Herrmann and Rouxl . Il990[ ) . Clearly, there is a 
cutoff in the stress field as r ^ 0. This may be caused 
by non-linearities in the material constitutive relations or 
by microstructure in the material such as the presence of 
cr ystallites. 

iHidalgo et al\ (|2002[ ) has introduced a fiber bundle 
that contains a power-law dependence on the distance 
from a failing fiber on the force redistribution in order to 
model the stress singularity seen around crack tips. The 
fiber bundle is implemented as a regular two-dimensional 
grid of parallel fibers clamped between two stiff blocks. 
Assuming that fiber j has just failed, a force transfer 
function 



Z 



where 



-=y- 



(244) 



(245) 



and / is the set of intact fibers, redistributes the forces. 
7 is treated as a parameter on the unit interval. There 
are two limiting cases, 7^0 which recovers the global 
load-sharing fiber bundle and 7 — > oo which recovers the 
local load-sharing model with nearest-neighbor stress al- 
leviation (Section lIV.Ap . The load increase on fiber i is 
hence given by 

where B is the set of failed fibers up to that point. 
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FIG. 35 Strength at failu re, Gc as a func t ion o f 7 in the vari- 
able range fiber bundle of lHidalgo eloZI (|2002l ). 
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FIG. 36 Strength at failure, ctc as a function of the number of 
fibers A'^ for di fferent range exponen ts 7 in the variable range 
fiber bundle of lHidalgo etal] (|2002l '). 



This model is too complex for analytical treatment and 
numerical simulations must be invoked. Around 7 = 
2.0 there is a transition in behavior between essentially 
global load sharing and local load sharing as described 
in Section flV. A I For 7 < 2.0, the maximum sustainable 
force scales the number of fibers in the bundle at the 
outset, N, whereas for 7 > 2.0, a A^/ In A^-behavior is 
observed as in Eq. (|242p . This is seen in Figs. [55l and [551 

The burst distribution shows a power-law distribution 
with exponent ^ = 5/2 again signaling global load shar- 
ing behavior for smaller values of 7. As 7 is increased, 
deviations from this behavior is seen. This must be inter- 
preted as a crossover towards local load sharing behavior 
as described in the previous section, see Fig. [ST] 

Lastly, the structure of the clusters of failed fibers at 
breakdown is studied. In the global load sharing model 
implemented in two dimensions does not yield anything 
particular. There is a percolation transition in the clus- 
ter size distribution when the relative density of failed 
fibers reaches the percolation threshold, but this has 
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FIG. 37 Burst size distribution for diff erent range exp onents 
7 in the variable range fiber bundle of lHidalgo et aZ.I |2002'). 
The data points where 7 < 2.0 can be fitted to a power law 
with exponent ^ — 5/2. For larger 7, the slope becomes 
steeper and resembling the behavior seen in the local load- 
sharing model, see Fig. 1311 



no particular s ignific ance in the evolution of the model. 
iHidako et al\ (I2002D find the cluster distribution having 
two distinct behaviors, depending on whether 7 is smaller 
than or larger than 2. There is no clear power law be- 



h avior, see Fig 

iHidalgo et al. (|2008h have studied an anisotropic ver- 
sion of this model. The force transfer function p44p is 
in this work generalized to 



Z 



(l-a)Ay2^.))7/2 



(247) 



where a is an anisotropy parameter. The behavior of this 
model turns out to be quite similar to that found in the 
is otropic model. 



iRaischel et all ( 2006) introduced a low c utoff in the 
threshold distribution f Pradhan e t al, 20 Q3) in t he vari- 
able range fiber bundle model of Hidal go et all |2002). 
They studied the burst distribution as a function of 7 
the burst distribution as a function of 7 and cutoff in the 
failure thresholds in terms of deformation, e^. Fig. [351 
summarizes their findings: For the explored values of 7 
in the range 2.0 < 7 < 6.0, a crossover from burst expo- 
nent f = 5/2to^ = 3/2is seen for small el, whereas for 
larger 7 sa 6.0, the burst distribution may be fitted to a 

v alue ^ = 9/2 for small cl - 

Newman and Gabrielovl ( 199l[ ) and [Newman et all 



(|l994l ) introduced a fiber bundle model where the fibers 
are hierarchically organized. The fibers are paired two by 
two. Each of these pairs is seen as an "order one" fiber. 
This pairing is repeated for the order one fibers, creating 
order two fibers and so on. Within each sub-bundle, the 
fibers are subject to equal load sharing. If the cumula- 
tive threshold distribution for the fibers at order zero is 
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FIG. 38 Size distribution of clusters of bro ken bonds at 
collap se in the variable range fiber bundle of iHidalgo et al\ 
(|2002l ). There is no evident power law behavior. 
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FIG. 39 Burst size distributions in the variable range fiber 
bundle model for different 7 and tL values: a) 7 = 2.0, b) 
2.5, c) 7 = 3.0 and d) 7 = 6.0. From iRaischel et al\ 



Looel ) 



Po{x) = P{x), then the threshold distribution at level 



Piix) = Po{x)[2Po{2x) ^ Po{x)] , 



(248) 



a result which is readily generalized to any level. Even 
though the starting point here is global load sharing, 
the approximation introduced by treating the fibers at 
each level as fibers with given thresholds, leads to the 
introduction of spatial load dependence in the model. To 
our knowledge, bursts have not been studied within this 
framework. 



C. Elastic medium anchoring 

In this Section we generalize the fiber bundle prob- 
lem to include more realistically the elastic response of 
the surfaces to which the fibers are attached. So far, 
these have been assumed to be infinitely stiff for the equal 
load sharing model, or their response has been modeled 
as very soft, but in a fairly unrealistic way in the local 
load sharing models, see Section HV.AI We will end up 
with a description that is somewhat related to the mod- 
els of the previous Sectio n, in particular themodel of 
IHidalgo et cd\ (|2002t) . In iBatrouni et al. I (|2002f ). a re- 
alistic model for the elastic response of the clamps was 
studied. The model was presented in the context of the 
failure of weldings. In this language, the two clamps were 
seen as elastic media glued together at a common inter- 
face. 

Without loss of generality, one of the media may be 
assumed to be infinitely stiff whereas the other is soft. 
When a force is applied to a given fiber, the soft clamp 
responds by a deformation falling off inversely to the dis- 
tance from the loaded fiber. Hence, the problem becomes 
one of solving the response of the surface with respect 
to a given loading of the fibers. Fibers exceeding their 
maximally sustainable load fail, and the forces and de- 
formations must be recalculated. The two clamps can 
be pulled apart by controlling (fixing) either the applied 
force or the displacement. The displacement is defined 
as the change in the distance between two points, one in 
each clamp positioned far from the interface. The line 
connecting these points is perpendicular to the average 
position of the interface. In our case, the pulling is ac- 
complished by controlling the displacement. As the dis- 
placement is increased slowly, fibers will fail, eventually 
ripping the two surfaces apart. 

We now concretize these ideas in a model. It consists 
of a two-dimensional square L x L lattices with periodic 
boundary conditions. The lower one represents the hard, 
stiff surface and the upper one the elastic surface. The 
nodes of the two lattices are matched, (i.e., there is no 
relative lateral displacement). The fibers are modelled as 
in the previous sections: elastic up to a threshold value 
which has been individually chosen for each fiber from 
some threshold distribution. The spacing between the 
fibers is a in both the x and y directions. The force that 
each fiber is carrying is transferred over an area of size 

to the soft clamp: As the two clamps are separated by 
controlling the displacement of the hard clamp relative to 
the zero level, D, the forces carried by the fibers increases 
from zero. When the force carried by a fiber reaches its 
breaking threshold, it breaks irreversibly and the forces 
redistribute themselves through the deformation of the 
soft clamp. Hence, the fibers are broken one by one until 
the two clamps are no longer in mechanical contact. The 
force, fi, carried by the ith fiber is given by 



/, = ~kiu, - D) 



(249) 



where k is the spring constant and is the deformation 
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of the elastic clamp at site i. All unbroken fibers have 
k = 1 while a broken fiber has fc = 0. The quantity 
{ui — D) is, therefore, the length, and since k — 1, also 
the force carried by fiber i. The deformation of the soft 
clamp is described by the coupled system of equations. 



(250) 



where the elastic Green function, Gj.,- is given by 
(jJohnsonl . Il985t iLandau and Litshitj . Il958l ) 



ha/2 



-a/2 



dx' dy' 



-a/2 J -a/2 x',y~ y') 



(251) 



In this equation, s is the Poisson ratio, e the elastic con- 
stant, and |j — j| the distance between sites i and j. The 
indices i and j run over all sites. The integration over 
the area is done to average t he force from the fiber s 



over this area. As remarked bv IBatrouni et oL 
the Green function (Eq. I25ip applies for a medium occu- 
pying the infinite half space. However, with a judicious 
choice of elastic constants, it may be used for a finite 
medium if its range is small compared to L, the size of 
the system. 

By combining Eq. (1249^ and Eq. (|250p . one obtains 



(I-hKG)/ = K^ 



(252) 



where matrix-vector notation is used. I is the x 
identity matrix, and G is the Green function represented 
as an x dense matrix. The constant vector D is 
dimensional. The diagonal matrix K is also x L^. 
Its matrix elements are either 1, for unbroken fibers, or 
for broken ones. 

Once Eq. ((252)) is solved for the force, /, Eq. (|250l) 
yields the deformations of the elastic clamp. 

Eq. (|252p is of the familiar form Ax = b. Since the 
Green function connects all nodes to all other nodes, the 
X matrix A is dense which puts severe limits on 
the size of the system that may be studied. 

The simulation proceeds as follows: One starts with 
all springs present, each with its stochastic breakdown 
threshold. The two media are then pulled apart, the 
forces calculated using the Conjugate Gradient algorithm 



(CG) (IBatrouni and Hansenl . [l988t iPress et al. 1 . 1 19921 ) . 



and the fiber which is the nearest to its threshold is bro- 
ken, i.e., the matrix element corresponding to it in the 
matrix K is zeroed. Then the new forces are calculated, 
a new fiber broken and so on until all fibers have failed. 

However, there are two problems that render the simu- 
lation of large systems extremely difficult from a numer- 
ical point of view. The first is that since G is x 
dense matrix, the number of operations per CG itera- 
tion scales like L^. Even more serious is the fact that 
as the system evolves and springs are broken, the matrix 
(I-l- fcG) becomes ill-conditioned. To overcome the prob- 
lematic scaling of the algorithm, the matrix-vector 
multiplications are done in Fourier space since the Green 



function is diagonal in this space. Symbolically, these 
multiplications may be written as follows. 



(I + KF-^FG)F-^F/ = KD 



(253) 



where F is the FFT (fast Fourier transform) operator 
and F^-*- its inverse (F^-'-F = 1). Since I and K are 
diagonal, operations involving them are performed in 
real space. With this formulation, the number of op- 
erations/iteration in the CG algorithm now scales like 
L'^ \n{L) rather than L^. 

To overcome the ill-conditioning of th e matrix (I+fcG) 
we need to precondition the ma trix ( Batrouni et al. I , 
119861: IBatrouni and Hansenl . Il988l ). This means that in- 
stead of solving Eq. (|253p , one solves the equivalent prob- 
lem 



Q(I + KF^FG)F^F/ ^ QKD 



(254) 



where we simply have multiplied both sides by the 
arbitrary, positive definite preconditioning matrix Q. 
Clearly, the ideal choice is Qo = (I + KG)~^ which 
would always solve the problem in one iteration. Since 
this is not possible in general, we look for a form for Q 
which satisfies the following two conditions: (1) As close 
as possible to Qo, and (2) fast to calculate. The choice 
of a good Q is further complicated by the fact that as 
the system evolves and fibers are broken, corresponding 
matrix elements of K are set to zero. So, the matrix 
(I + KG) e volves from the i nitial form (I -I- G) to the 
final one I. IBatrouni et al. I ( 2002f) did not find a fixed 
Q that worked throughout the entire breakdown process. 
They therefore chose the form 

Q = I - (KG) + (KG)(KG) - (KG)(KG)(KG) + ... 

(255) 

which is the Taylor series expansion of Qo = (I+KG)~^. 
For best performance, the number of terms kept in the 
expansion is left as a parameter since it depends on the 
physical parameters of the system. It is important to em- 
phasize the following points, (a) As fibers arc broken, the 
preconditioning matrix evolves with the ill-conditioned 
matrix and, therefore, remains a good approximation of 
its inverse throughout the breaking process, (b) All ma- 
trix multiplications involving G are done using FFTs. 
(c) The calculation of Q can be easily organized so that 
it scales like nL^ ln(i) where n is the number of terms 
kept in the Taylor expansion, equation (j255p . The result 
is a stable accelerated algorithm which scales essentially 
as the volume of the system. 

Fig.HOlshows the force-displacement curve for a system 
of size 128 x 128 and elastic constant e = 10. Whether we 
control the applied force, F, or the displacement, D, the 
system will eventually suffer catastrophic collapse. How- 
ever, this is not so when e — 100 as shown in Fig. 311 
In this case, only controlling the force will lead to catas- 
trophic failure. In the limit when e — > oo, the model be- 
comes the equal load-sharing fiber bundle model, where 
F = {1 — D)D. In this limit there are no spatial cor- 
relations and the force instability is due to the decreas- 
ing total elastic constant of the system making the force 
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FIG. 41 Force-displacement curve, 128 x 128 systems with 
e = 100. 
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FIG. 43 Burst distribution for 128 x 128, e = 100. The slope 
of the straight line is —2.5. 



on each surviving bond increase faster than the typical 
spread of threshold values. No such effect exists when 
controlling displacement D. However, when the elastic 
constant, e, is small, spatial correlations in the form of 
localization, where fibers that are close in space have a 
tendency to fail consecutively, do develop, and these are 
responsible for the displacement instability which is seen 
in Fig. HOI 

We now turn to the study of the burst distribution. 
Figs. [32] and 23] show the burst distribution for e = 10 
and 100. In both cases we find that the burst distribu- 
tion follows a power law w ith an exponent £ = 2.6 ±0.1. 
It was argued in Ref. (Batrouni et al. 1 12002D that the 
value of £ in this case is indeed 5/2 as in the global 
load sharing model. These two figures should be com- 
pared with Fig. [37] showing the burst distribution in 
the variable range fi ber bundle model of Hidalgo et al., 
(|Hidalgo et aLl . l2002l ). where £ = 5/2 is recovered as long 
as the range exponent 7 is small, rendering the forces 
long range among the fibers. 

As the failure process proceeds, there is an increasing 



competition between local failure due to stress enhance- 
ment and local failure due to local weakness of material. 
When the displacement, D, is the control parameter and 
e is sufficiently small (for example e = 10), catastrophic 
failure eventually occurs due to localization. The on- 
set of this localization, i.e., the catastrophic regime, oc- 
curs when the two mechanisms are equall y important . 
This may be due to self-organized criticality ( Bak et a7~l . 
I1987D occurrin g at this point. I n ord er to test whether 
this is the case. iBatrouni et al. 1 (I2002D measured the size 
distribution of broken bond clusters at the point when D 
reaches its maximum point on the F — D characteristics, 
i.e. the onset of localization and catastrophic failure. 
The analys is was performed using a Ho shen-Kopelman 
algorithm fetauffe r and Aharonvl . ll99 4') . The result is 
shown the in Fig.mj for 56 disorder realizations, L — 128 
and e = 10. The result is consistent with a power law 
distribution with exponent —1.6, and consequently with 
self organization. If this process were in the universality 
class of percolation, the exponent would have been 2.05. 
Hence, we are dealing with a new universality class in 



38 



10' 



_g 10" 



10-' 



10"' 



■0 



0^ 



10 100 

Cluster size 

FIG. 44 Area distribution of zones where glue has failed for 
systems of size 128 x 128 and elastic constant e — 10. The 
straight line is a least square fit and indicates a power law 
with exponent —1.6. 



this system. This behavior should be contrasted to the 
on e seen in the va riable range fiber bundle model studied 
bv lHidalgo efdl ((2002: ) where no power law distribution 
was found (see Fig. [351 in Section [lEll) . 



V. FIBER BUNDLES IN MATERIAL SCIENCE AND 
OTHER APPLICATIONS 

The aim of this Section is to demonstrate how the fiber 
bundle model may be used as a tool for studying both 
important phenomena occuring in materials such as fa- 
tigue and for studying important classes of materials, no- 
tably fiber reinforced comp osites. A general review has 
recent ly been written by iMishnaevskv and BrOndstedl 
( 20091 ) on this subject. We also discuss appHcations of 
fiber bundle models in other contexts, ranging from traf- 
fic modeling to earthquake dynamics. 



A. Time dependent failure: Fatigue or creep phenomena 

Materials may undergo time dependent deformation 
under steady load. Sometimes when a load is applied, 
though surviving at first stage, the system fails after 
long time. This type of failure is referred as the fatigue - 
failure or creep rupture ( ChakrabartiL 1 199^ : iLawnl . 1 1993( ) . 
The fatigue- failure is basica l ly a t hermally activated pro- 
cess ( Phoenix and Tiernevl . Il983[ ) and originates at the 
atomic level of the fibers where the molecules are in ran- 
dom thermal vibrations. Eventually a molecule acquires 
sufficient thermal energy to overcome the local energy 
barrier and slips relative to other molecules. The fre- 
quency of such events is greatly enhanced by increases in 
temperature, stress and impurity level. After a molecular 
slip or rupture, neighboring molecules become overloaded 
and the failure rate increases. These molecular failures 
accumulate locally and produce microcracks within the 



material. Also micro-cracks can grow at the crack-tips 
with time du e to chemical diffusion in the atmosphere 
( Lawnl . [l993( ) which helps the growth of fractures. These 
failures nucleate around the defects in the solid and the 
failure behavior and its statistics therefore crucially de- 
pends on the disorder or impurity distribution within the 
sample. The system then fails under a stress less than 
its normal strength (cc) and the failure time (r) depends 
on both, the applied load and the impurity level. 

Fatig ue- failure in fiber bundle model was first studied 
by IColeman (,1956,. 1957a ) considering different classes 
of fibers and several breakdown rules. The probabilis- 
tic analysis gives the lifetime distribution under vari- 
ous loading conditions: Constant load, loads propor- 
tional to time and periodic loads. The 'time depen- 
dent fatigue' and 'cycle dependent fatigue' both are ad- 
dressed in this work introducing the concept of 'mem- 
ory' effect, i.e., the load history can affect the failure of 
fibers. Such time dependent failure in fiber bundles have 
been c onsiderably extended and generalized by Phoenix 
et al. (fN ewman and Phoenixl. 120011: iPhoenixL \l97§ . 1 19791: 
|Ph oenix and Tierney, 1983) for equal load sharing (ELS) 
and local load sharing (LLS) bundles. The approximate 
fatigue life time distributions have been achieved through 
probabilistic analysis introducing power law breakdown 
rule and exponential breakdown rule at the molecular 
level ( Phoenix and Tiernevl . [l983[ ) . 

When dry fibers are replced by viscoelastic elements 
having time dependent deformation properties, fiber bun- 
dle model exhibi ts creep behavior ( Hidalgo et all . [2003 
iKun et al\ . [2003h in terms of the macroscopic response 
under constant external load. There exists a critical 
load (or stress) below which the deformation attains a 
constant value (infinite lifetime) and above the criti- 
cal load, deformation increases monotonically resulting 
global failure (finite lifetime). Another extension of clas- 
sical fiber bundle model, the cq i itinuo us damage model 
( Hidalgo et aLl . [200H iKun et a/.l . l2003f ). captures similar 
creep behavior assuming that a fiber can fail more than 
once and at each failure its stiffness gets reduced by a 
constant factor. In both of these models the lifetime of 
the bundle diverges at the critical load following robust 
power law variation with the applied load. 

Als o, few ex perimen ts (iBaneriee and ChakrabartiL 
120011 : iKun et all l2007l : iPauchard and Meunieri . 11993! ) 
have been performed to observe the failure time of ma- 
terials and its statistics. The effect of thermal activa- 
tion and disordered noise on the failure have recently 
been measured for material breakdown and approxi- 



mate fatigue b ehavior has been obtained (iGuarino et al 
1999, 2002: Ku n et all. l2007t iPradhan and Chakrabarti . 



2003bi : iRoux. . ,2000l:IScorretti a/.IJ200l[) using fiber bun- 



dle models. 

In this section, we discuss several approaches of achiev- 
ing fatigue-failure behavior in equal load sharing fiber 
bundle models. The approaches differ basically in the 
way how time-dependence has been incorporated in the 
failure process. 
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1. Thermally induced failure in fiber bundles 

The influence of noise on macroscopic failure 
in fiber bundle model, has been studied numeri- 



callv (ICiliberto et al\ . 1200 ll : iGuarino et d[ Il999l 120021 : 
IScorretti et g/J . 12001 ). using both disorder noise and ther- 
mal noise. The strength of each fiber is characterized by a 

(c) 

critical stress xl which is a random variable that follows 
a normal distribution of mean x^"^^ and variance KTri: 



(256) 



where K is the Boltzmann constant and Nd is the dis- 
order noise. Again each fiber is subjected to an addi- 
tive time dependent random stress Axi{t) which follows 
a zero mean normal distribution of variance KT: 



Ax,{t)^ NTit,KT). 



(257) 



Here Nt is the thermal noise. 

Due to the equal load sharing scheme, if a number n{t) 
of fibers are broken at time t after force F is applied on 
the bundle, the local force on each of the remaining fibers 
will be: 



Xi{t) 



XqN 

N ~n{t) 



Ax^{t), 



(258) 



where N is the total number of fibers in the intact bundle 
and Xq = F/N is the initial force per fiber. If KT = 0, 
the model reduces to the static one. In that case the ap- 
plied force is increased linearly from zero to the critical 
value Fc above which the whole bundle breaks. There- 
fore, at a constant force F, the bundle breaks in a single 
avalanche only if _F > F^ otherwise it will never break. 
If KT then the system can break at an applied force 
F < Fc due to the thermal effect. Such thermal failure of 
the model has bee n studied numer i cally as a function of 
Xq, K T and KT^ (jCiliberto et all l200ll : iGuarino et aZI . 
I1999D . 

When KT ^ 0, the failure time Tf as a function of 
(1 — xq)'^ follows an exponential law for any fixed value 
of KTd-. 

Tf ^ exp [a(l - xof] , (259) 

where a is a fitting parameter [Fig. I45r a)]. which is a 
function of KT. At constant stress, the failure time de- 
pends on thermal noise KT as: 



Tf To exp 



(kt) 



(260) 



where A is a function of x [Fig. |45{b)]. Similar result 
has also been observed in case of a heterogeneous fiber 
bundle [Fig. gS]. 

One can com pare these results with Pomeau's theory 
( Pomeaul Il992i) for the failure time of solids: 



/ pdy(d-l) 

Tf = To exp ag ■- 

\ ' KTcff P^''''\ 



(261) 
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FIG. 45 Failure time r/ of an homogeneous bundle {KTd = 
0) in a creep test, (a) Tf as a function of the normalized 
force (1 — xo)'^ for several values of thermal noise variance 
KT: (cross) KT = 0.0045; (circle) KT = 0.006; (triangle) 
KT — 0.01. (b) Tf as a function of 1/KT for several values 
of So: (box) Xq — 0.45; (circle) xq = 0.54; (triangle) xq = 0.7. 
Continuous lines in (a) and (b) are the fits with Eq. (I259II and 
Eq. (|260p , respectively. Adapted from iCiliberto et a/.. (,2001< 'l. 



where Ps is the imposed stress, Tq is a constant, F^. the 
surface energy, Y the Young modulus, ag is a constant 
which depends on the geometry, T^f / is an effective tem- 
perature and d the dimensionality of the system. This 
theory is based on the phys ical argument that t h ermal 
activatio n of micro cracks ( Golubovic and Fend . Il99ll : 
iPomeaul [19 9 2) is responsible for the macroscopic fail- 
ure of the material. It can be noted that the functional 
dependence of Tf on stress is different for fiber bundle 
model and for solids. The main reason of which is the 
different geometry of the stress distribution in the fiber 
bundle and the sohds. 

These numerical studies suggest that disorder noise 
amplifies the effect of thermal noise and reduces the de- 
pendence of Tf on the temperature and this can explain 
the recent experimenta l obse rvation s on microcrystals 
(jPauchard and Meunieil . Il993l ). gels (iBonn et al\. fl998() 
and macroscopic composite materials (jGuarino et 
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FIG. 46 Failure time r/ of an heterogeneous bundle {KTd 7^ 
0) in creep test. The life time r/ is plotted as a function of 
1/KT at xq = 0.45. Different symbols correspond to different 
values of KTa: (Circle) KTa = 0; (triangle) KTg = 0.02; 
(box) KTd = 0.04. Adapted from lCiliberto ei al\ (120011 1. 



II999D . 

The numerical observations described above have been 
confirmed la ter through an analytic investigation by 
iRouxl (|2000l) . In a homogeneous (no disorder in fiber 
strengths) fiber bundle model, force or stress on each 
fiber is 



X ^ xo + ri, 



(262) 



where xq = F/N and rj is the random noise with a Gaus- 
sian distribution 



1 



y/2TTKT 



exp 



2Kt)' 



(263) 



with zero mean and variance KT. Now the probability 
that one fiber survives after t time step is 



(264) 



where P is the cumulative probability. Then the proba- 
bility that the entire bundle can survive after t time step 



PAr(t) = [1 - P(l - Xo)]^*. 

Therefore the average failure time is 

, > -1 



7Vln[l-P(l-a;o)]' 



(265) 



(266) 



After the first fiber breaks, the situation remains the 
same with a smaller bundle and larger stress. Thus the 
average failure time after i — 1 broken fiber is 



-1 



{N - i) In {1 - P[l - Nxo/iN - i)]} ' 



(267) 



Now the total failure time can be obtained by taking 
the sum over all i as 



N 



^ (A^ - i) In {1 - P[l - Nxo/{N 



(268) 



When A^ is large, one can replace the sum by a contin- 
uous integral: 



N 



-N 



iN-y)ln{l-P[l-Nxo/{N-y)]} 
— 1 dz 



dy 



ln[l-P(l-z)] z ■ 

To achieve a closed form equation it has been con- 
sidered that the above sum is dominated by the time re- 
quired for breaking the first fiber when xq is much smaller 
than the maximum load Xc and when KT << 1. Then 
P can be considered to be much smaller than 1. Now the 
derivative of the average time with respect to xq gives 



d{rf) 



1 



1 



(269) 



dxQ xo In [1 - P(l - Xq)] a:;oP(l-a;o)' 

When (1 — .-eq)^ ^ KT the error function can be ex- 
panded as 



P(l-a;o) 



/xrexp(-(ig^) 



[1 + OiKT)]. (270) 



V27r(l - xo) 

Finally taking into account the dominating terms, one 
gets 



, , V2ttKT 
{Tf) = z cxp 



Xo 



2KT 



(271) 



These analytic expressions a re identical to what 
have been observed ea r her (ICiliberto et 'dl I2OOII: 
iGuarino et all . Il999l 120021 : IScorretti et aLl . l200lh innu- 
merical simulations. A similar analysis (|Rouxl . I2OOOI ) 
shows that when fiber strengths are distributed (hetero- 
geneous case), the average first failure time can be ex- 
pressed as 



27r (l-a:o) 
VK{T + 6) 



N 



exp 



'2K{T + e) 



(272) 



where 9 is an effective temperature that is added to the 
temperature T due to the disorder. Hence, the disorder 
leads to an effective tem perature Tg// = T + Q, and this 
is what IScorretti et all (|2001.1 proposed: time indepen- 
dent heterogeneities of the system modify the effective 
temp erature. 

Pol iti et al ] (1200 Ih estimate the total time to failure 
for the thermally activated fiber bundle model and find 
the same behavior as in Eq. (|272[) : The disorder adds 
a constant to the temperature of the fiber bundle. In 



iGuarino et al. I (|2006[) . these results are generalized to the 
two-dimensional fuse model. 
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FIG. 47 The simulation results showing variation of average 
failure time r/ against (a) external stress a and (b) against 
noise T {K — 1 here), for a homogeneous bundle containing 
N = 10^ fibers. The insets show the variation of the frac- 
tion U of unbroken fibers with time t for different T values 
[1.2 (cross) and 1.0 (plus)] in (a) and different a values [0.15 
(cross) and 0.12 (plus)] in (b). The dotted and dashed lines 
represent the theoretical result Eq. (|277|l . 



2. Noise induced failure in fiber bundles 

Not only the temperature, but several other fac- 
tors can result in fatigue- failure in materials: Weather 
effects , che mical effects etc. ( Chakrabartil Il994l : 
iLawnl . 



19931). Recently, there h as been an attempt 



( Pradhan and Chakrabartl l2003bl ) to incorporate all the 
noise effects through a single parameter KT {K is Boltz- 
mann constant) which can directly influence the failure 
probability of the individual elements. Such a failure 
probability p{a, KT) at any applied stress a, induced by 
a non-zero noise KT, has been formulated as: 



p{a,KT) 



1, 



< CT < a; 
a > X 



(273) 

where x is the failure strength of an element. Clearly, 
the failure probability increases as a and KT increases. 
Without any noise {KT = 0) the model is trivial: The 
bundle remains intact for stress a < ac and it fails com- 
pletely for cr > (7c where Uc is the critical stress value. 

At KT and under any stress a (< CTc), some fibers 
(weaker fibers) fail and the total load has to be sup- 



ported by the surviving fibers, which in turn enhances 
their stress value, inducing further failure. The bundle 
therefore fails at a < ac after a finite time r/ . 

In case of homogeneous bundle (all fibers have the same 
strength x), the critical stress value for the bundle is 
CTc = X. Then t he time dynamics at an applied st ress a 
can be written ( Pradhan and Chakrabart i l2003bl) as 



Ut+i = Ut 



^-p\jy.KT 



(274) 



where Ut is the fraction of total fibers remains intact 
after time step t. In the continuum limit, we can write 
the above recursion relation in a differential form 



dU 



a 

— exp 

Oc 



1 ( Gc 



(275) 



The solution gives the failure time 



Tf = I at — — exp 
/o cr 



1 

KT 



I '^''P KT ( a 



U 



dU. 



Hence, for a < Uc 
Tf = KT exp 



1 

ICT 



exp 



f- 

\aKT 



(276) 
(277) 



Again, for cr > (Tc, one gets Ut+i — from Eq. 
(|274p . giving Tf — 0. For small KT and as ti — > ac, 
Tf ~ iCrexp[(crc/cr- l)/-ft:T]. These results agree 
qualitatively with the re cent exp erimental ob servations 
( Baneriee and Chakrabart i. 200ll: iGuarino et ai . 2002). 

In order to investigate the fatigue behavior 
in heterogeneous fiber bundles, uniform distri- 
bution of fiber strengt h s has been considered 
( Pradhan and Chakrabartil 2003bl ). The noise- 

induced failure probability has the similar form : 
p{(J,KT) = exp [--^ (f - 1)] for < (T < x and 
p((T, KT) — 1 for a > X, where x denotes the strength of 
the individual fibers in the bundle. Now it is difficult to 
tackle the pr oblem analytically. Howeve r , Mont e Carlo 
simulations ( Pradhan and Chakrabartil l2003b[ ) show 
(Fig. [48|) the variations of average failure time (t/) with 
noise [KT) and stress level (cr): 



Tf — KT exp 



1 

XT 



exp 



1 



aKT KT 



(278) 

where, Uc is the critical stress. This phenomenological 
form Eq. (I278|) is indeed very close to the analytic result 
E g. (I276p for the homog e neous fiber bundle. 

iNewman and Phoenixl (|200lh considered the breaking 
dynamics in a fiber bundle where each fiber has a failure 
probability p determined by the loading time t and the 
load cr on it as p{t,(j) = 1 — exp[—ta''], and analyzed 
the life time (tf ) distribution. For 1/2 < p < 1, ELS 
and LLS models have identical Gaussian distribution for 
tf. For p > 1, LLS shows extrem e statistics, while 
ELS gives Gaussian behavior, see e.g.. lCurtin and Sched 
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FIG. 48 Simulation results of fatigue behavior: (a) average 
failure time Tf vs. noise T {K = 1 here) for three different 
stress values and (b) r/ vs. a for three different noise val- 
ues. The bundle contains A*" = 10^ fibers with uniformly dis- 
tributed strength thresholds. The time variations of fraction 
of surviving fibers are shown in the insets. The dotted lines 
in (a) and (b) correspond to the fit wit h expression (Eg. I278p 
where ffe — 0.245 (exact value=l/4 ijBhattacharvva et all 
l2003l : lPradhan et ffld . [200^ 1. 



(|l997f l. lYoshioka et al. I ([200H) also studied similar tf 
distribution and their scaling property with load and 
temperature variations. 

3. Creep rupture in viscoelastic fiber bundles 



Creep behavior ha s been achieved ( Hidalgo et all . 
I2Q02I: Ik un et all I2OO3I) in a bundle of viscoelastic fibers, 
where a fiber is modeled by a Kelvin- Voigt element (see 
Fig. SHI and results in the constitutive stress-strain rela- 
tion 



Poi + Ye. 



(279) 



Here (Tq is the applied stress, e is the corresponding 
strain, /3o denotes the damping coefficient, and Y is the 
Young modulus of the fibers. 

In the equal load sharing mode, the time evolution 
of the system under a steady external stress CTq can be 



FIG. 49 The viscoelastic fiber bundle: intact fibers are mod- 
eled by Kelvin- Voigt eleme nts which consists of a spring and 
a dashpot in parallel. From lKun et all (|2003l ). 



described by the equation 



1 - Pie) 



(3oe + Ye. 



(280) 



As one can expect intuitively, there is a critical load 
CTc for the system and Eq. (|280p suggests two distinct 
regimes depending on the value of the external load cFo- 
When (To is below the critical value CTc Eq. (|280p has a 
fixed-point solution es , which can be obtained by setting 
£ = in Eq. (|280ll 



ao = Yes[l-P{es)]. 



(281) 



In this case the strain value converges to when t 00, 
and no macroscopic failure occurs. But, when Co > Cc, 
no fixed-point solution exists. Here, e remains always 
positive, that means in this case, the strain of the sys- 
tem e(t) monotonically increase s until the syst em fails 
globally at a finite time t f (Hidalgo et all |2002[ ) . 

The solution of the differential equation Eq. (|280p gives 
a complete description of the failure process. By separa- 
tion of variables, the integral becomes 



t = (3o / de 



1 - Pie) 



'ao-Ye[l-Pie)] 



+ C, 



(282) 



where C is integration constant. 

Below the critical point (Tq < the bundle slowly 
relaxes to the fixed-point value e^. The characteristic 
time scale of such relaxation process can be obtained by 
analyzing the behavior of e(t) in the vicinity of e^. After 
introducing a new variable Sq as 6oit) = — e(t), the 
differential equation can be written as 



dt 



Y 

To 



1 - 



£sPi£s) 
1 - P{es 



(283) 



Clearly, the solution of Eq. (|283p has the form Sq 
exp [— i/r], with 



Y 



1 - 



£sP(£s) 



l-^(£s) 



(284) 
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FIG. 50 The behavior of the time to failure tf for uniform 
and WeibuU distributions with two different WeibuU index 
for the ELS case. All the three curves are parallel to each 
other on a double logarithmic plot with an exponent close to 
0.5, in agreement w ith the general result Eq. (|288|l . From 
iHidalgo et all(|2002l ). 



where r is the characteristic time of the relaxation pro- 
cess. 

The variation of the relaxation time r with the exter- 
nal driving near the critical point ctc is crucial for any 
dynamical system. Since cro(es) has a maximum of the 
value (Tc at Ec, in the vicinity of one can use the ap- 
proximation: 



(285) 



where the multiplication factor A depends on the cumu- 
lative distribution P. Using the approximation Eq. (|285p . 
it can be shown from Eq. (|284p .that 



-1/2 



for (To < 



(286) 



Therefore, the relaxation time diverges following a uni- 
versal power law with an exponent —1/2. Note that a 
dry fiber bundle model, under constant load shows simi- 
lar power law divergence (see Section III. A). 

How does the system behave above the critical point? 
The behavior can be analyzed in the same way when (Tq 
is close to (Tc- Then one can write ctq = Cc + Actq, where 
AiTo << (Tc- It is obvious that the relaxation steps are 
too many when e{t) becomes close to Ec- Therefore, the 
integral in Eq. (|282p , is dominated by the region close to 
£c- Using Eq. (|285p . the integral in Eq. (|282p becomes 



l-P{e) 



(287) 



Evaluating the integration over a small e interval in the 
vicinity of £c, one gets 



tf « (cTo - ^Tc)"^/^ 



for 



(To > (Tc. 



(288) 



Thus, tf has a two-sided power law divergence at iTc with 
a universal exponent — ^ independent of the specific form 




FIG. 51 T he fibers are modeled by Maxwell elements. From 
iKun et al\ (|2003i ) 



of the disorder distribution P{£), similarly to r in case 
of dry fiber bundle (see Section III. A). 

To check the validity of the universal power law be- 
havio r of tf, simulations were performed I Hidalgo et all 
with various disorder distributions, i.e. uniform 
distribution and the WeibuU distribution of the form 
P(e) = 1 — exp [— (e/A)''], where A is the characteris- 
tic strain and p is the shape parameter. The simulation 
results (Fig. [50]) are in excellent agreement with the an- 
alytic results. 



4. Creep rupture in a bundle of slowly relaxing fibers 

A slow relaxati on following fiber fa i lure can also lea d 
to creep behavior ( Hidalgo et aLl . [200ll : lKun et al. |,[2003). 
In this case, the fibers are linearly elastic until they break, 
but after breaking they undergo a slow relaxation pro- 
cess. Therefore when a fiber breaks, its load does not 
drop to zero instantaneously. Instead it undergoes a slow 
relaxation process and thus introduces a time scale into 
the system. As the intact fibers are assumed to be lin- 
early elastic, the deformation rate is 



X 



(289) 



where x denotes the stress, Sf denotes the strain and Y 
is the Young modulus of the fibers. In addition, to cap- 
ture the slow relaxation effect, the broken fibers with the 
surrounding matrix material are modeled by Maxwell ele- 
ments (Fig. lSip . i.e. they are assumed as a serial coupling 
of a spring and a dashpot. Such arrangement results in 
the following non-linear response 



^+Bxl 



(290) 



where Xb is the time dependent stress and Eb is the time 
dependent deformation of a broken fiber. The relaxation 
of broken fibers is characterized by few parameters, Sb, B, 
and m, where is the effective stiffness of a broken 
fiber, the exponent m characterizes the strength of non- 
linearity and i? is a constant. 
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In equal load sharing mode, when external stress ctq 
is applied, the macroscopic elastic behavior of the com- 
posite can be represented by the coii stitutive equation 
(|Hidalgo et ail . 1200 ll : iKun et all l2003^ 



ao = x{t) [1 - P{x{t))] + Xb{t)P{x{t)). 



(291) 



Here Xb{t) is the amount of stress carried by the broken 
fibers and P{x{t)) and 1 — P{x{t)) denote the fraction of 
broken and intact fibers at time t, respectively. 

By construction (Fig. I5ip . the two time derivatives 
have to be always equal 



(292) 



Now the differential equation for the time evolution of 
the system can be obtained using Eq. (|29ip . Eq. (|290p 
and Eq. (p92)) as 



-4 



B 



1 



1 

(To - X[l 



1 



P{x) 
P{x)] 



P{x) 

p{xY 



(X - (To) 



P{x) 



(293) 



Similar to the viscoelastic model, two different regimes 
of x{t) can be distinguished depending on the value of 
(Jo'- If the external load is below the critical value ctc a 
fixed-point solution Xs exists which can be obtained by 
setting X = in Eq. ([^55)1 



CTo ^Xs[l- P{Xs 



(294) 



This means that the solution x{t) of Eq. (|293p converges 
asymptotically to Xs resulting in an infinite lifetime tjoi 
the system. When the external load is above the critical 
value, the deformation rate e = x/Y remains always pos- 
itive, resulting in a macroscopic failure in a finite time tf. 
Now we focus on the universal behavior of the model in 
the vicinity of the critical point. Below the critical point 
the relaxation of x{t) to the stationary solution Xs can 
be presented by a differential equation of the form 



dt 



-'0 ' 



(295) 



where 5o denotes the difference (So(0 = Xg — x{t). Hence, 
the characteristic time scale r of the relaxation process 
only emerges if m = 1. Also, in this case relaxation time 
goes as T ~ {gc — Oo)~^^'^ when approaching the critical 
point from below. However, for m > 1 the situation is 
different: relaxation process is characterized by Soit) — 
where a ^ with ao ac- 
Again, close to the critical point, it can also be shown 
that the lifetime t f shows power law divergence when the 
external load approaches the critical point from above 



tf (cTo - fc) 



-(m-l/2) 



for (To > (Tc 



(296) 



The exponent is universal in the sense that it does not 
depend on the disorder distribution. However it depends 




FIG. 52 Lifetime as a function of the distance from the 
critical point, ao — (Jc, for two different values of the parameter 
m. The number of fibers in the bundle is N = 10^. From 
iKun et al\ l|2003f) . 




FIG. 53 Asphalt samples set up for experimental testing of fa- 
tigue failure. Figure (a) shows how fatigue failure under these 
experimental conditions are modeled using a fi ber bundle 
model and (b) shows a post-failure sample. From lKun et al\ 
(|2007l ). 



on the exponent m, which characterizes the non-linearity 
of broken fibers. 

As a check, numerical simulations ha ve been performed 
(|Hidalgo et all l200ll: iKun et all l2003l) :br several differ- 
ent values of the exponent m (Fig. [52]). The slope of the 
fitted straight lines agrees well with the analytic predic- 
tions (Eq.iinni). 



5. Fatigue-failure experiment 

An interesting experimental and theoretic al study of 
fatigu e failure in asphalt was performed by iKun et al\ 
( 20071 ). The experimental set-up is shown in Fig. (53] 
The cylindrical sample was subjected to cyclic diametric 
compression at constant load amplitude (Tq, and the de- 
formation e as a function of the number of cycles Ncyde 
was recorded together with the number of cycles Nf at 
which catastrophic failure occurs. 

Fig. [54] shows deformation e as a function of num- 
ber of cycles Ncyde for two different load amplitudes 
o'o/'^c = 0.3 and 0.4. Here ac is the tensile strength 
of the asphalt. Fig. [55] shows the number of cycles to 
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FIG. 54 Deformation e as a function of number of cycles 
N cycle, showing both experimental data and theoret i cal cu rves 
based on the fiber bundle model. From iKun et al\ (|2007l '). 



catastrophic failure as a function of the load amplitude 
(To/cTc. This curve shows three regimes, the middle one 
being characterizable by a power law, 




FIG. 55 Number of cycles at catastrophic failure, Nf as a 
function of load amplitude ao/oc- The two curves show ex- 
perimental and fiber bundle model data respectively. From 
iKun et al\ IpOOTlV 



Or 



(297) 



This is the Basquin regime { Li and Metcalj . 1200! 



Si et al 



Kun et 



20021: ISornette et all 120021 : ISureshT Il99lh . 



(|2007l ) find a' = 2.2 ± 0.1 for the asphalt sys- 



tem. 

I n order to model the behavior found in Figs. [54] and 
[551 [Kun et ail (2003) introduce equal load-sharing fiber 
bundle model as illustrated in Fig. [55k . Each fiber 
1 < i < N is subjected to a time dependent load Xi{t). 
There are two failure mechanisms present. Fiber i fails 
instantaneously at time t when Xi{t) for the first time 
reaches its failure threshold ti. However, there is also 
a damage accumulation mechanism characterized by the 
parameter Ci{t). In the time interval dt, fiber i accumu- 
lates a damage 



dci{t) = axi{ty dt 



(298) 



where a > is a scale parameter and 7 > is an exponent 
to be determined. Hence, the accumulated damage is 



c,(t) = a / x^{t'y dt' . 



(299) 



When Ci (t) for the first time exceeds the damage accumu- 
lation threshold, Si, fiber i fails. The thresholds ti and Si 
are chosen from a joint probability distribution pt.sit, s). 
[Kun et all ( 2007( ) make the assumption that this distri- 
bution may be factorized, pt,s{t, s) = pt(t)ps{s). 

In addition to damage accumulation, there is yet an- 
other important mechanism that needs to be incorpo- 
rated in the model: damage healing ( Jo et a/l [20081 : 



[Si et all [200I . A time scale r is associated with this 
mechanism, and the ELS average force-load equation 



ao(t)[l ~ Pt{x{t))]x{t) 



(300) 



where Pt{t) = j^Ptit') dt' is the cumulative instanta- 
neous breaking threshold probability, is generalized to 



CTo(t) = 



I -P. 



vit'ydt' 



[1-Pt{x{t))]x{t) , 
(301) 



where Ps{s) = Ps{s') ds' . [Kun et al\ ([2007[ ) show that 
the Basquin law (|297p may be derived analytically from 
Eq. pOip leading to a' = 7. The solid curves in Fig. [54] 
are fits of the theoretical curves based on Eq. (|30ip to 
the experimental data for the e vs. N cycle- Likewise, Fig. 
[551 shows the fit of the theoretical Nf vs. ctq /^c to the 
experimental data. For this fit, 7 = 2.0 and r — 15000. 



B. Precursors of global failure 

A fundamental question in strength considerations of 
materials is when does it fail. Are there signals that can 
warn of imminent failure? This is of uttermost impor- 
tance in e.g. the diamond mining industry where sudden 
failure of the mine can be extremely costly in terms 
of lives. These mines are under continuous acoustic 
surveillance, but at present there are no tell-tale acoustic 
signature of imminent catastrophic failure. The same 
type of question is of course also central to earthquake 
prediction and initiates the search for precursors of global 
(catastrophic) failure events, see e.g.. [Sahimi and Arbabil 
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FIG. 56 Variation of and t'"^ with applied stress for 
a bundle having TV = 50000 fibers. Uniform distribution of 
fiber thresholds is considered and averages are taken over 1000 
sample. The dotted straight lines are the best linear fits near 
the critical point. 
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FIG. 57 Variation of a-nd with applied stress for 
a single bundle having A'^ = 10000000 fibers with uniform 
distribution of fiber thresholds. Straight lines represent the 
best linear fits near the critical point. 



(1X991 I1996D. iPradhan and Chakrabartil (HqoI), 
iPradhan and Chakrabartil (|2006t)" The precursor 



parameters essentially reflect the growing correlations 
within a dynamic system as it approaches the failure 
point. As we show, sometimes it is possible to predict 
the global failure points in advance. Needless to mention 
that the existence of any such precursors and detailed 
knowledge about their behavior for major catastrophic 
failures like earthquakes, landslides, mine/bridge col- 
lapses, would be of supreme value for our civilization. 
In this sub-section we discuss some precursors of global 
failure in ELS models. We also comment on how can 
one predict the critical point (global failure point) from 
the precursor parameters. 

(a) Divergence of susceptibility and relaxation time 

As we discussed earlier (Section III. A) in case of ELS 
fiber bundles the susceptibility (x) and the relaxation 
time (r) follow power laws (exponent = —1/2) with 
external stress and both diverge at the critical stress. 
Therefore if we plot ^-nd with external stress, 
we expect a linear fit near critical point and the straight 
lines should touch X axis at the critical stress. We in- 
deed found similar behavior (Fig. I56p in simulation ex- 
periments after taking averages over many sample. 

For application, it is always important that such pre- 



diction can be done in a single sample. For a single bun- 
dle having very large number of fibers, similar response 
of X and r have been observed . The estimation (through 
extrapolation) of the failure point is also quite satisfac- 
tory (Fig. [FT)) . 

(h) Pattern of breaking-rate 

When we apply load on a material body, it is impor- 
tant to know whether the body can support that load or 
not. The similar question can be asked in fiber bundle 
model. We found that if we record the breaking-rate, 
i.e, the amount of failure in each load redistribution - 
then the pattern of breaking-rate clearly shows whether 
the bundle is going to fail or not. For any stress be- 
low the critical state, the breaking-rate follows exponen- 
tial decay (Fig. [58]) with the step of load redistribution 
and for stress values above critical stress it is a power 
law followed by a gradual rise (Fig. [55)) . Clearly, at 
critical stress it follows a robust power law with expo- 
nent value —2 that can be obtained analytically from 
Eq. ()^0)) . As we can see from Fig. [55] that when the 
applied stress value is above the critical stress, breaking- 
rate initially goes down with step number, then at some 
point it starts going up and continues till the complete 
breakdown. That means if breaking-rate changes from 
downward trend to upward trend - the bundle will fail 
surely -but not immediately after the change occurs -it 
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FIG. 58 Log- log plot of breaking-rate with step of load redis- 
tribution for 7 different stress values. Circles are for stresses 
below the critical stress and triangles are for stresses above 
the critical stress. The simulation has been performed for a 
single bundle with A'^ = 10000000 fibers having uniform dis- 
tribution of fiber thresholds. The dotted straight line has a 
slope —2. 



takes few more steps and number of these steps decreases 
as we apply bigger external stress (above the critical 
value). Therefore, if we can locate this minimum in the 
breaking-rate pattern, we can save the system (bundle) 
from breaking down by withdrawing the applied load im- 
mediately. We have another important question here: Is 
there any relation between the breaking rate minimum 
and the failure time (time to collapse) of the bundle? 
There is indeed a universal relationship which has been 
explored recently ( Pradhan and Hemmed . l2009l ) through 
numerical and analytical studies: For slightly overloaded 
bundle we can rewrite Eq. (|8T|) as 



Ut = i- Vitan(A*t-S* 



(302) 



where 



A* = tan"^(2Ve) and B* = tan"i(l/2Ve 

From Eq. (|302p follows the breaking rate 



Rit) = = V^A* cos-2(A*t - B*). 

at 

R{t) has a minimum when 

0= — ocsin(2A*t-2S*), 
at 



(303) 
(304) 



which corresponds to 



tn = 



B* 

'a* 



(305) 



(306) 



When criticality is approached, i.e. when e — > 0, we have 
A* —> 0, and thus to ~^ oo, as expected. 
We see from Eq. (|302l) that Ut ^ for 




FIG. 59 The breaking rate R{t) vs. step t (upper plot) 
and vs. the rescaled step variable t/tf (lower plot) for the 
uniform threshold distribution for a bundle of A'^ = 10^ fibers. 
Different symbols are used for different excess stress levels 
a — ac'- 0.001 (circles), 0.003 (triangles), 0.005 (squares) and 
0.007 (crosses). 



This is an excellent approximation to the integer value at 
which the fiber bundle collapses completely. Thus with 
very good approximation we have the simple connection 
tf — 2to. When the breaking rate starts increasing we 
are halfway (see Fig. [55)1 to complete collapse! 



(c) Crossover signature in avalanche distribution 

The bursts or avalanches can be recorded from outside 
-without disturbing the ongoing failure process. There- 
fore, any signature in burst statistics that can warn of im- 
minent system failure would be very useful in the sense of 
wide scope of applicability. As discussed in Section III.B, 
when the avalanches are recorded close to the global fail- 
ure point, the distribution shows (Fig. I60p a different 
power law (^ ~ 3/2) than the one (^ = 5/2) characteriz- 
ing the size distribution of all avalanches. This crossover 
behavior has been analyzed analytically in case of ELS 
fiber bundle model and similar crossover behavior is also 



tt 



{B* + tan"i(l/2Ve) /A* = 2B*/A*. (307) seen (jPradhan et all 1200611 in the burst distribution and 
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FIG. 60 The avalanche size distributions for different values 
of xo in ELS model with uniform fiber strength distribution. 
Here bundle size — 50000 and averages are taken over 10000 
sample. Two power laws (dotted lines) have been drawn as 
reference lines to compare the numerical results. 



energy distribution of the fuse model which is an estab- 
lished model for studying fracture and breakdown phe- 
nomena in disordered systems. The crossover length be- 
comes bigger and bigger as the failure point is approached 
and it diverges at the failure point (Eq. I192p . In some 
sense, the magnitude of the crossover length tells us how 
far the system is from the global failure point. Most im- 
portant is that this crossover signal does not hinge on 
observing rare events and is seen also in a single sys- 
tem (see Fig. [T^ . Therefore, such crossover signature 
has a strong potential to be used as a useful detection 
tool. It should be mentioned that a recent observation 
(|Kawamural I2006D suggests a clear change in exponent 
values of the local magnitude distributions of earthquakes 
in Japan, before the onset of a mainshock (Fig. [6T|) . This 
observation has definitely strengthened the possibility of 
using crossover signals in burst statistics as a criterion 
for imminent failure. 



C. Fiber reinforced composites 

As we have seen, fiber bundle models provide a fertile 
ground for studying a wide range of breakdown phenom- 
ena. In some sense, they correspond to the Ising model in 
the study of magnetism. In this Section, we will review 
how the fiber bundle models are generalized to describe 
composites containing fibers. Such composites are of in- 
creasing practical importance, see e.g.. Fig. [S^ 

The status of modeling fiber reinforc ed c o mpos- 
ites has recently been revi e wed ([Mishnaevskvl . l2007t 
iMishnaevskv and BrOndstedl . l2009l ). These materials 
consist of fibers embedded in a matrix. During ten- 
sile loading the main part of the load is carried by the 
fibers and the strength of the composite is governed to 
a large extent by the strength of the fibers themselves. 
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FIG. 61 Crossover signature in the local magnitude distribu- 
tions of earthquakes in Japan. The exponent of the distribu- 
tion during 100 days before a mainshock i s about 0.60, muc h 
smaller than the average value 0.88. From|K awamur al lj2006D . 
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FIG. 62 This dental bridge is made from a fiber reinforced 
composite with braided fibers made from polyethylene. Due 
to the braided structure, this composite is four times tougher 
than a composit e made from the same materi als, but with- 
out the braiding (|Karbhari and Strasslerl . [2007h (Courtesy H. 
Strassler, Univ. of Maryland Dental School and V. Karbhari, 
Univ. of Alabama, Huntsville). 



The matrix material is chosen so that its yield threshold 
is lower than that of the fibers which are embedded in 
it. Common materials used for the fibers are aluminum, 
aluminum oxide, aluminum silica, asbestos, beryllium, 
beryllium carbide, beryllium oxide, carbon (graphite), 
glass (E-glass, S-glass, D-glass), molybdenum, polyamide 
(aromatic polyamide, aramid), Kevlar 29 and Kevlar 49, 
polyester, quartz (fused silica), steel, tantalum, titanium, 
tungsten or tungsten monocarbide. Most matrix materi- 
als are resins as a result of their wide variation in proper- 
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ties and relatively low cost. Common resin materials are 
epoxy, phenolic, polyester, polyurethane and vinyl ester. 
When the composite is to be used under adverse condi- 
tions such as high temperature, metallic matrix materials 
such as aluminum, copper, lead, magnesium, nickel, sil- 
ver or titanium, or non-metallic matrix materials such 
as ceramics may be used. When the matrix material is 
brittle, cracks open up in the matrix perpendicular to 
the fiber direction at roughly equal spacing. In metal- 
lic matrix materials, plasticity sets in at sufficient load. 
Lastly, in polymer matrix composites, the matrix typi- 
cally responds linearly, but still the fibers carry most of 
the load due to the large compliance of the matrix. When 
a fiber fails, the forces it carried are redistributed among 
the surviving fibers and the matrix. If the matrix-fiber 
interface is weak compared to the strength of the fibers 
and the matrix themselves, fractures develop along the 
fibers. When the matrix is brittle, the fibers bridging 
the developing crack in the matrix will, besides binding 
the crack together, lead to stress alleviation at the ma- 
trix crack front. This leads to the energy necessary to 
propagate a crack further increases with the l ength of 
the crack, ( Sorensen and Jacobsen . 19981 . |200(J) i.e., so- 
called R-curve behavior (jLawnl . 119931 ). When the bridg- 
ing fibers fail, they typically do so through debonding at 
the fiber-matrix interface. This is followed by pull-out, 
see Fig. m 

The Cox shear lag model forms the basis for the stan- 
dard tools used for analyzing break down in fiber rein- 
forced composites ( Choul Il997l : ICox . 1952). It consid- 
ers the elastic response of a single fiber in a homegeous 
matrix only capable of transmitting shear stresses. By 
treating the properties of the matrix as effective and 
due to the self-consistent response by the matrix ma- 
terial and the rest of the fibers, the Cox mod el be- 
comes a mean- field model ( Raisanen et al. 1. ll997l ). Ex- 
tensions of the Cox single-fiber model to debonding 
and sli p at the fiber-matrix interface have been pub- 
lished (I Aveston and Keiivl . Il973t iBudianskv et all . 1 19861: 
lHsuehl . ll990lll992l ). In 1961 the single-fiber calculation 
of Cox was extended to two-dimensional unidirectional 
fibers in a compliant matr ix, i.e., a m a trix i ncapable of 
carrying tensile stress, by iHedgepet h* ('1961^ ). In 1967, 
this ca lculation was followed up bvlHcdgcpct h and Dvkg 
(1963) for three-dimensional unidirectional fibers placed 
in a square or hexagonal pattern. They found the aver- 
age stress intensity factor (i.e., the ratio between local 
stress in an intact fiber and the applied stress) to be 



2i 



2i + l 



(308) 



after k fibers failing. In his Ph.D. thesis, iFichteil 
(I1969D extended these calculations to aligned arrays 
of broken fibers mixed with intact fibers. This ap- 
proach was subsequently generalized to systems where 
the matrix has a non-z ero stiffness and hence is 
able to transmit stress ( Beverlein and Landid . Il999t 
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FIG. 63 Post mortem micrograph showing a fiber reinforced 
composite where the matrix has undergone brittle failure fol- 
lowed by_JajluTe_of_tlM_^^ through debonding. 
From iKarbhari and Strassled (|2007l '). 



iLandis and McMeekingj . Il999|). Viscoelasti c ity o f the 
matrix has been included bv iLagoudas et all ( 1989( ) and 
iBeverlein and Phoenix! (|l998fK 

ICurtid (119911) demonstrated that when the fibers re- 
spond under global load-sharing conditions, a mean- 
field theory may be constructed where the breakdown 
of the composite is reduced to that of the fai l ure of 
a single fiber i n an effective matrix (ICurtin j_ 1993 ; 
Hild et a/.l.[T99llHild and Feillardl. ll997HRoux and Hild . 



2002| ). IWagner and Eita^ (1993) studied the redistribu 



tion of forces onto the neighbors of a single failing fiber 
within a two-dimensional uni-directional composite using 
the shear-lag model, finding that within this scheme, the 
stress-enhancement is less pr onounced than earlier cal- 
culations had shown. IZhou and Wagner (1999 , ,2000.) in- 
troduced a multi-fiber failure model including debonding 
and frictional effects at the fiber-matrix interface, finding 
that the stress intensity factor would decrease with in- 
creasing interfiber distance. An important calculational 
principle, the Bre ak Influence Superpo s ition Technique 
was introduced bv Sastrv and Phoenixl ( 19931 ) based on 
the method of Kac hanovl (|l985l ) in order to handle mod- 
els with multiple fiber failures. The technique consists 
in determining the transmission factors, which give the 
load at a given position along a given fiber due to a unit 
negative load at the single break point in the fiber bun- 
dle. The multiple failure case is then constructed through 
superposition of these single-failure transmission factors. 
This method has proven very efficient from a numerical 
point of view, and has been generalized through a s e 



ries of later pape r s, see ( Beverlein and Phoenbl 1997allb 



Beverlein et al. I . Il996t Landis et al. . 2000; iLi et al. 



20061) 



llbnabdelialil and Cm^ (Il997a'.'b'), 
I5cia and Curtin (2001) and [Xia e t al. 



'Curtinl (I1998D . 

(2002.) analyzed 
the interaction between multiple breaks in uni-directional 
fibers embedde d in a matrix using a la ttice Green func- 
tion technique ( Zhou and Curtinl , l200lh to calculate the 
load transfer from broken to unbroken fibers including 
fiber-matrix sliding with a constant interfacial shear re- 
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FIG. 64 The discretization of a three-dim ensional unidirec- 
tional fiber reinforced composite used by IZhou and Curtiiil 
(|200ll ). 



sistance r, given by either a debonded sliding interface or 
by matrix shear yielding. The differential load carrying 
capacity of the matrix is assumed to be negligible. In 
the following we describe the Zhou and Curtin approach 
in some detail. The load-bearing fibers have a strength 
distribution given by the cumulative probability 



P(cr,L) = 1 



(309) 



of failure over a length of fiber L experiencing a stress ct, 
where 



^0 Vfo, 



(310) 



where p is the WeibuU index. When a fiber breaks, the 
load is transferred to the unbroken fibers. We will return 
to the details henceforth. The newly broken fiber slides 
relatively to the matrix. The shear resistance r provides 
an average axial fiber stress along the single broken fiber 



(t(z) = min I , cr°(z) 

r 



p{z) 



(311) 



at a distance z from the break, where r is radius of the 
fibers and cr"(z) is the axial fiber stress prior to the failure 
at point z. This defines a length scale 



2t 



(312) 



The total stress change within a distance ±Zs of the break 
is distributed to the other fibers. A key assumption in 
what now follows is that the total stress in each plane 
z is conserved: The stress difference a°{z) — a{z) is dis- 
tributed among the other intact fibers at the same z-level. 

In order to set up the lattice Green function approach, 
the system must be discretized. Each fiber, oriented in 



the z direction, of length Lz is divided into elements 
of length S = L^/Nz. The fibers are arranged on the 
nodes of a square lattice in the xy plane so that there is 
a total of Nf = x Ny fibers. The lattice constants in 
the X and y directions are Ux and Uy, respectively. Each 
fiber is labeled by n where 1 < n < Nf. This is shown 
in Fig. 1641 The stress on fiber n in layer m along the 
z direction is given by an,m- Fiber n at layer m may 
be intact. It then acts as a Hookean spring with spring 
constant kt responding to the stress (Jn,m- If fiber n has 
broken at layer m, it carries a stress equal to zero. The 
third possibility is that fiber n has broken elsewhere at 
m', and layer m is within the slip zone. It then carries a 
stress 



2t5. 



"^'l:Cr°.m 



Pn 



(313) 



which is the discretization of Eq. (j31ip with zero spring 
constant. 

Each element m of fiber n has two end nodes associ- 
ated with it. At all such nodes, springs parallel to the xy 
plane are placed linking fiber n with its nearest neighbors. 
These springs have spring constant kg . The displacement 
of the nodes is assumed confined to the z direction only. 
Zhou and Curtin denote the displacement of node con- 
necting element m with element m -I- 1 of fiber n, 
and the displacement of node linking element m with el- 
ement m — 1 of fiber n, u~ ^. The force on element m of 
fiber n from element m of fiber n -I- 1 is 

fm{n; n + 1) = fcs(w^+i.m — "^t.m) + ^s(u„+i.„i ~ Un,m) ■ 

(314) 

The reader should compare the following discussion 
with that which was presented in Section HV. CI We now 
assume that it is only layer m = that carries any dam- 
aged or slipped elements, the rest of the layers m ^ are 
perfect. Let u = {u^^}- If a force f = {f^.m} apphed 
to the nodes, the response is 



u = G f , 



(315) 



where G is the lattice Green function. Given the dis- 
placements from solving this equation combined with Eq. 
(|314p , the force carried by each broken element is found. 
The inverse of the lattice Green function is D = G~^. 
The elements of D are either zero, kg or fcj, reflecting 
the status of the springs; undamaged, slipping or bro- 
ken. When there are no breaks in layer to = 0, we define 
D" = (G")-i, and 5T> = D^-D. Hence, 5T> plays a role 
somewhat similar to the matrix K defined in Eq. (|252p . 
By combining these definitions, Zhou and Curtin find 



G= (l--G°<5D)-i G" 



(316) 



The matrices G and D have dimension N x N where 
N = iV^ X Ny X Nz- By appropriately labeling the rows 
and columns, the matrices D and G may be written 



G 




(317) 
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where the {2NxNy) x {2NxNy) matrix Gdd couples ele- 
ments within the layer m = 0, where all the damage is 
located. The matrix Gpp couples elements within the rest 
of the layers. These are undamaged — "perfect" . The 
two matrices Gdp and Gpd provide the cross couplings. 
The matrix ST) becomes in this representation 



6T> 



6Ddd 




(318) 



Combining this equation with Eq. (|316p gives 



Gdd = {I -SUddyG'^dd, (319) 

where the intact Green function G^^ may be found an- 
alytically by solving Eq. (|315p for the intact lattice in 
Fourier space {q}, 



sin^ (a^) + ks sin^ (^) -I- kt sin^ 



(320) 



By using that = — / n, 0, Zhou and Curtin find 
that 

(321) 

where and n refer to the upper and lower node at- 
tached to element n is layer m = 0. Before completing 
the model, the WeibuU strength distribution, Eqs. (|309p 
and (|310p . must be discretized. Each element {n,m) is 
given a maximum sustainable load Sn,m from the cumu- 
lative probability 



P/(s) = 1 - , 

where = {Lo/Sy/faQ. 

The breakdown algorithm proceeds as follows: 



(322) 



1. A force per fiber set equal to the smallest break- 
ing threshold, /o = minn.m Sn,m, is applied to the 
system. 

2. The weakest fiber or fibers are broken by setting 
their spring constants to zero. 

3. Decrease the stresses in the element below and 
above the just broken fibers according to Eq. (|313p . 



4. Solve Eq. (|319p for the layers in which the breaks 
occurred. 

5. Calculate the spring displacements in the layers 
where the breaks occurred using Eq. (|32ip and 
an effective applied force /+„j = /o — Pn,m- The 
force on each intact spring in such a layer is then 



6. With the new stresses, search for other springs that 
carry a force beyond their thresholds Sn,m- If such 
springs are found, break these and return to (2). 
Otherwise proceed. 

7. Search for the spring which is closest to it break- 
ing threshold. This spring is the one with A = 
min„ , m(sn,m//n,m)- Incrcasc the load by a factor 
A?7, where rj is equal to or somewhat larger than 
unity. This factor is present to take into account 
the non-linearities introduced in the system due to 
the slip of the fibers. 

8. Proceed until the system no longer can sustain a 
load. 

By changing the ratio kt / kg between the moduli of the 
springs in the discretized lattice, it is possible to go from 
fiber bundle behavior essentially evolving according to 
equal load sharing (ELS) to local load sharing (LLS). 

Whereas the computational cost of finite-element cal- 
culations on fiber reinforced composites scales with the 
volume of the composite, the Break Influence Superposi- 
tion Technique and the Lattice Green function technique 
scale with the number of fiber breaks in the sample. This 
translates into systems studied by the latter two tech- 
niques can be orders of magnitud e larger than the former 
(|lbnabdelialil and Curti"^ . ll997aD . 

After this rather sketchy tour through the use of fiber 
bundle models as tools for describing the increasing im- 
portant fiber reinforced composite materials, we now turn 
to the use of fiber bundle models in non-mcchanical set- 
tings. 



D. Failure phenomena in networks, traffic and earthquake. 

The typical failure dynamics of the fiber bundle model 
captures quite faithfully the failure behavior of several 
multicomponent systems like the communication or traf- 
fic networks. Similar to the elastic networks considered 
here, as the local stress or load (transmission rate or 
traffic currents) at any part of the network goes beyond 
the sustainable limit, that part of the system or the net- 
work fails or gets jammed, and the excess load gets redis- 
tributed over the other intact parts. This, in turn, may 
induce further failure or jamming in the system. Because 
of the tectonic motions stresses develop at the crust - 
tectonic plate resting (contact) regions and the failure 
at any of these supports induces additional stresses else- 
where. Apart from the healing phenomena in geological 
faults, the fiber bundle models have built-in features to 
capture the earthquake dynamics. Naturally, the sta- 
tistically established laws for earthquake dynamics can 
be easily recast into the forms derived here for the fiber 
bundle models. 

We will consider here in some more details these three 
applications specifically. 
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1. Modelling network failures 



Th e fiber bundle model has been applied ()Kim et all . 
[2005) to study the cascading failures of network struc- 
tures , like Erdos-Renyi networks ( Erd5s and Renyil 
I1959D . known as ER netwo r ks an d Watts-Strogatz net- 
works ( Watts and Strogat3 . 1 199 4^. known as WS net- 
works, to model the overloading failures in power grids, 
etc. Here, the nodes or the individual power stations are 
modelled as fibers and the transmission links between 
these nodes are utilized to transfer the excess load (from 
one broken fiber or station to another). 

The load transfer of broken fibers or nodes through the 
edges or link s of the underlying network is governed by 
the LLS rule j Harlow and Phoeni^ . ll978HPlioenixl . ll979t 
ISmith and Phoenix . 19"8lh . Under a non-zero external 
load Na, the actual stress ai of the intact fiber i is given 
by the sum of a and the transferred load from neighbor- 
ing broken fibers. The local load transfer, from broken 
fibers to intact fibers, depends on the load concentration 
factor Ki = Gi/a with Ki = 1 + "^/rnj/kj, where the 
primed summation is over the cluster of broken fibers di- 
rectly connected to i, mj is the number of broken fibers 
in the cluster j, and kj is the number of intact fibers 
directly connected to j. 

Let the external stress a be increased by an infinites- 
imal amount Sa starting from a = 0. Fibers for which 
strength < Kia break iteratively until no more fibers 
break. For each increment of a, the size s(cr) of the 
avalanche is defined as the number of broken fibers trig- 
gered by the increment. The surviving fraction C/(cr) of 
fibers can be written as 



(323) 



One can also measure directly the response function x, 
or the generalized susceptibility, denoted as 



dU 
da 



(324) 



The critical value Uc of the external load, can be defined 
from the condition of the global breakdown U{ac) — . 

The critical value CTc and the susceptibility x have been 
calculated numerically for the model under the LLS rule 
on various network structures, such as the local regu- 
lar network, the WS network, the ER network, and the 
scale- free networks (Fig. [SS]). The results suggest that 
the critical behavior of the model on complex networks 
is completely different from that on a regular lattice. 
More specifically, while CTc for the FBM on a local regu- 
lar network vanishes in the thermodynamic limit and is 
described by CTc ~ 1/ ln( A^) for finite-sized systems (see 
the curve for p = in Fig. (HSl corresponding to the WS 
network with the rewiring probability p = 0), ac for all 
networks except for the local regular one does not dimin- 
ish but converges to a nonzero value as N is increased. 
Moreover, the susceptibility diverges at the critical point 
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FIG. 65 (a) The system size (TV) dependence of critical points 
((Tc) for various networks with A'' = 2*, 2^, . . . , 2^^ vertices, 
(b) The susceptibility for the networks with A'^ — 2^*. p and 
7 are the rewiring probability in the WS networks and the 
exponent of degree distribution P(fc) ~ fc ' respectively. The 
data points are obtained fr om the averages over 10* (10^ for 
N = 2^^) ensembles. From lKim eTall (|2005l 1. 



as X ~ (cc — cr)~^/^, regardless of the networks, which is 
again in a sharp contrast to the local regular network [see 
Fig. [SHb)]. The critical exponent 1/2 clearly indicates 
that the FBM under the LLS rule on complex networks 
belongs to the same un i versal ity class as that of the ELS 
regime (|Pradhan et al although the load-sharing 

rule is strictly local. The observed variation for CTc is only 
natural for LLS model. In an LLS model if n successive 
fibers fail (each with a finite probability pf), then the 
total probability of such an event is NpJ{l — pf Y as the 
probability is proportional to the bundle size N . If this 
probability is finite, then rt ~ IniV (for any finite pj, 
the failure probability of any fiber in the bundle). For 
a failure of n successive fibers, the neighboring intact 
fibers get the transferred load ~ na which if becomes 
greater than or equal to their strength, it surely fails giv- 
ing (Tc ^ \/n ^ \/\nN (see Section IV. A). 

The evidence that the LLS model on complex networks 
belongs to the universality class of the ELS model is 
also found (|Kim et al ], [looi) in the avalanche size dis- 
tribution D(A): Unanimously observed power-law be- 
havior i3(A) ~ A~^/^ (a) for all networks except for 
the local regular one (the WS network with p = 0) 
is in perfect agreement with the behavior for the ELS 
case ( Hemmer and Hansenl . Il992[ ). On the other hand, 
the LLS model for a regular lattice has been shown 
to ex hibit completely diffe r ent avalanche size d i stribu - 
tion ( Hansen and Hemmerj . 11994 iKloster et al. 
Also one can observe a clear difference in terms of the fail- 
ure probability F{a) defined as the probability of failure 
of the whole system at an external stress a. While F{a) 
values for LLS on complex networks fall on a common 
line, L LS on regular net work show s a distinctly different 
trend ( Kim et see also iDivakaran and Duttal 

(|2007cD . 
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2. Modelling trafPic jams 

One can apply the equal load sharing fiber bundle 
model to study the traffic failure in a system of paral- 
lel road network in a city. For some special distributions, 
like the uniform distribution, of traffic handling capaci- 
ties (thresholds) of the roads, the critical behavior of the 
jamming transition can be studied analytically. This, in 
fact, is exactly comparable with that for the asymmet- 
ric simple exclusion process in a single channel or road 
(IChakrabartil2006D . 

Traffic jams or congestions occur essentially due to the 
excluded volume effects (among the vehicles) in a sin- 
gle road and due to the cooperative (traffic) load sharing 
by the (free) lanes or roads in multiply connected road 
netwo rks (see e.g., IChowdhurv et aZ.f (|2000l ): Ide Silveiral 
( 19991 )) . Using FBM for th e traffic network, it has been 
shown (|Chakrabartill2006^ that the generic equation for 
the approach of the jamming transition in FBM corre- 
sponds to that for the Asymmetric Simple Exclusion Pro- 
cesses (ASEP) leading to the trans port failure transitio n 
in a single channel or lane (see e . g . . I S t inchcomb j ( 20051 )). 

Let the suburban highway traffic, while entering the 
city, get fragmented equally through the various narrower 
streets within the city and get combined again outside the 
city (see Fig. IM]). If Iq denotes the input traffic current 
and It is the total output traffic current, then at the 
steady state, without any global traffic jam. It = Iq- In 
case It falls below Iq, the global jam starts and soon 
It drops to zero. This occurs if Iq > Ic, the critical 
traffic current of the network, beyond which global traf- 
fic jam occurs. Let the parallel roads within the city 
have different thresholds for traffic handling capacity: 

, , • ■ • , foi' the N different roads (the n-th road 
gets jammed if the traffic current i per road exceeds jc„). 
Initially i — Iq /N and increases as some of the roads get 
jammed and the same traffic load Iq has to be shared 
equally by a lower number of unjammed roads. Next, we 
assume that the distribution p{ic) of these thresholds is 
uniform up to a maximum threshold current (correspond- 
ing to the widest road traffic current capacity), which is 
normalized to unity (sets the scale for /c). 

The jamming dynamics in this model starts from the 
n-th road (say in the morning) when the traffic load i per 
city roads exceeds the threshold ic„ of that road. Due to 
this jam, the total number of unconjested roads decreases 
and the rest of these roads have to bear the entire traffic 
load in the system. Hence the effective traffic load or 
stress on the uncongested roads increases and this com- 
pels some more roads to get jammed. These two sequen- 
tial operations, namely the stress or traffic load redis- 
tribution and further failure in service of roads continue 
till an equilibrium is reached, where either the surviv- 
ing roads are strong (wide) enough to share equally and 
carry the entire traffic load on the system (for Iq < Ic) 
or all the roads fail (for Iq > Ic) and a (global) traffic 
jam occurs in the entire road network system. 

This jamming dynamics can be represented by recur- 



Suburban Highway 



Suburban Highway 



City roads/lanes 

FIG. 66 Tlie highway traffic current Iq gets fragmented into 
uniform currents i in each of the narrower roads and the roads 
having threshold current ic„ < i get congested or blocked. 
This results in extra load for the uncongested roads. We 
assume that this extra load per unconjested roads gets equally 
redistributed and gets added to the existing load, causing 
further blocking of some more roads. 



sion relations in discrete time steps. Let Ut{i) be the 
fraction of unconjested roads in the network that survive 
after (discrete) time step t, counted from the time t — 
when the load (at the level Iq = iN) is put in the system 
(time step indicates the number of stress redistributions) . 
As such, Ut{i = 0) = 1 for all t and Ut{i) 1 for t = 
for any i; Ut{i) = U*{i) 7^ for i ^ oo if /q < /c, and 
Ut{i) = for < ^ 00 if /o > I^. 

Here Ut{i) follows a simple recursion relation 



t+i 



Iq 
UfN 



or, Ut+i = 1 - — 



(325) 



The critical behavior of this model remains the same as 
discussed in Section III. A in terms of i and the exponent 
values remain unchanged: a~l/2^(3 = 6, rj~l for 
all these equal (traffic) load sharing models. 

In the simplest version of the asymmetric simple ex- 
clusion process transport in a chain (see Fig. [57]), the 
transport corresponds to movement of vehicles, which is 
possible only when a vehicle at site I, say, moves to the 
vac ant site / -I- 1. The transport current / is then given 
by ([Stinchcombi,[200i) 



Jp'i{l-p'i+i), 



(326) 



where pj denotes the site occupation density at site I and 
J denotes the inter-site hopping probability. The above 
equation can be easily recast in the form 



Pi+i 



(327) 
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J 



/ i+j 



FIG. 67 The transport current / in the one dimensional lane 
or road is possible if, say, the l-th site is occupied and the 
(l + l)-th site is vacant. The inter-site hopping probability is 
indicated by J. 



where a = I/J. Formally it is the same as the re- 
cursion relation for the density of uncongested roads in 
the FBM model discussed above; the site index here in 
ASEP plays the role of time index in FBM. Such ex- 
act correspondence indicates identical critical behavior 
in both the cases. The same universality for different 
cases (different threshold distributions) in FBM suggests 
similar behavior for ot her equivalent ASEP cases as well 
( Bhattachariee I . l2007t) . For e xtension of the mo del to 
scale free Traffic Networks, see IZheng et al\ (|2008f) . 



As mentioned already, because of the tectonic motions, 
stresses develop at the crust-tectonic plate contact re- 
gions and the entire load is supported by such regions. 
Failure at any of these supports necessitates load redis- 
tributions and induces additional stresses elsewhere. In 
fact the avalanche statistics discussed in Section III.B 
can easily explain the Gutenberg-Richter law (Eq. I328P 
with the identification ^' = ^ — 1. Similarly, the decay 
of the number of failed fibers N{t) — N{1 ~ Ut) at the 
critical point, given by N{t) ^ t^^ (see e.g., Eq. BU]) . 
crudely speaking, gives in turn the Omori law behavior 
(Eq. I329P for the fiber bundle model, with the identifi- 
cation rj — 1 + rj. 

For some recent d iscussions on further s tudies along 
these lines, see e.g., iTurcotte and Glassed ( 20041) . The 
stick-slip moti o n in th e Burridge-Knopoff model (see e.g., 
ICarlson et al\ ( 1994I )). where the blocks (representing 
portion of the solid crust) connected with springs (repre- 
senting the elastic strain developed due to tectonic mo- 
tion) are pulled uniformly on a rough surface, has the 
same feature of stress redistribution as one or more blocks 
slip and the d ynamics was ma pped onto a ELS fiber bun- 



sup ana tne d ynamics was ma ppea oi 
die model bv lSornettd (|1992D . The power law distribu- 
tions of the fluctuation driven bursts around the critical 
points have been interpreted as the above two statistical 
laws for earthquakes. 



3. Modelling earthquake dynamics 

The Earth's outer crust, several tens of kilometers 
in thickness, rests on tectonic shells. Due to the 
high temperature-pressure phase changes and conse- 
quent ionizations of the metallic ores, powerful magneto- 
hydrodynamic convective flows occur in the earth's man- 
tle at several hundreds of kilometers in depth. The tec- 
tonic shell, divided into about ten mobile plates, have got 
relative veloc ities of the or der of a few centimeters per 
year, see e.g.. IScholj (|2002[) . 

The stresses developed at the interfaces between the 
crust and the tectonic shells during the (long) sticking 
periods get released during the (very short) slips, causing 
the releases of the stored elastic energies (at the fault 
asperities) and consequent earthquakes. 

Two well known phenomenologically established 
laws governing the earthquake statistics are (a) the 
Gutenberg-Richter law 



(328) 



relating the number density {N) of earthquakes with the 
released energy greater than or equal to E; and (b) the 
Omori law 



d{Af{t))/dt = l/i" , 



(329) 



where M denotes the number of aftershocks having mag- 
nitude or released energy larger than a preassigned small 
but otherwise arbitrary threshold value. 



VI. SUMMARY AND CONCLUDING REMARKS 

The fiber bundle model enjoys a rare double position in 
that it is both useful in a practical setting for describing 
a class of real materials under real working conditions, 
and at the same time being abstract enough to function 
as a model for exploring fundamental breakdown mecha- 
nisms from a general point of view. Very few models are 
capable of such a double life. This means that a review of 
the fiber bundle model may take on very different char- 
acter depending on the point of view. We have in this 
review emphasized the fiber bundle model as a model for 
exploring fundamental breakdown mechanisms. 

Failure in loaded disordered materials is a collective 
phenomenon. It proceeds through a competition between 
disorder and stress distribution. The disorder implies a 
distribution of local strength. If the stress distribution 
were uniform in the material, it would be the weakest 
spot that would fail first. Suppose now that there has 
been a local failure at a given spot in the material. The 
further away we go from this failed region, the weaker 
the weakest region within this distance will be. Hence, 
the disorder makes local failures repel each other: They 
will occur as far as possible from each other. However, 
as the material fails locally, the stresses are redistributed. 
This redistribution creates hot spots where local failure 
is likely due to high stresses. Since these hot spots oc- 
cur at the boundaries of the failed regions, the effect of 
the stress field is an attra ction between the local failures 
( Roux and Hansenl . fl990l) . Hence, disorder and stress has 
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opposite effect on the breakdown process; repulsion vs. 
attraction, and this leads to competition between them. 
Since the disorder in the strength of the material, leading 
to repulsion, dimishes throughout the breakdown pro- 
cess, whereas the stresses create increasingly important 
hot spots, it is the stress distribution that ends up dom- 
inating towards the end of the process. 

The fiber bundle model catches this essential aspect of 
the failure process. Depending on the load redistribution 
mechanism, the quantitative aspects change. However, 
qualitatively it remains the same. In the ELS case, there 
are no localized hot spots, all surviving fibers are loaded 
the same way. Geometry does not enter into the redis- 
tribution of forces, and we may say that the "hot spots" 
include all surviving fibers. This aspect gives the ELS 
fiber bundle model its mean field character, even though 
all other fluctuations are present, such as those giving 
rise to bursts. 

As shown in Section III, the lack of geometrical aspects 
in the redistribution of forces in the ELS model enables 
us to construct the recursion relations (e.g., Eq. [5]) which 
capture well the failure dynamics. We find that the even- 
tual statistics, governed by the fixed points for a < ac, 
the average strength of the bundle, essentially shows a 
normal critical behavior: order parameter O ~ (ctc — <j)", 
breakdown susceptibility x ^ (cc ^ cr) ^nd relaxation 
time r = k|(Tc — "'I"* with a — 1/2 = f3 = 6 and 
k^/k+ = lniV/27r for a bundle of N fibers, where sub- 
scripts -|- and — refers to post and pre critical cases re- 
spectively. The statistics of fluctuations over these av- 
erage behavior, given by the avalanche size distributions 
D{A) ^ with ^ = 3 for discrete load increment and 
— 5/2 for quasi-static load increment in such EIS cases. 
The critical stress CTc of the bundle is of course nonuni- 
versal and its magnitude depends on the flber strength 
distribution. 

For the LLS model, we essentially find (see Section 

IV. A), the critical strength of the fiber bundle <Jc ^ 
1 / In TV which vanishes in the macroscopic system size 
limit. The avalanche size distribution is exponential for 
such cases. For range-dependent redistribution of load 
(see Section IV. B) one recovers the finite value of CTc and 
the ELS-like mean field behavior for its failure statistics. 

Extensions of the model to capture creep and fatigue 
behavior of composite materials are discussed in Section 

V. A. Precursors of global failure are discussed in Section 
V.B. It appears, a detailed knowledge of the critical be- 
havior of the model can help very precise determination 
of the global failure point from the well defined precur- 
sors. Section V.C provides a rather cursory review of 
models of fiber reinforced composites. These models go 
far beyond the simple fiber bundle model in complex- 
ity and represent the state of the art of theoretical ap- 
proaches to this important class of materials. However, 
as complicated as these models are, the philosophy of 
the fiber bundle model is still very much present. Finally 
we discussed a few extensions of the model to failures in 
communication networks, traffic jams and earthquakes in 



Section V.D. 

As discussed here in details, the fiber bundle model 
not only gives an elegant and profound solution of the 
dynamic critical phenomena of failures in disordered sys- 
tems, with the associated universality classes etc, but also 
offers the first solution to the entire linear and non-linear 
stress-strain behavior for any material up to its fracture 
or rupture point. Although the model had been intro- 
duced at about the same time (1926) as the Ising model 
for static critical phenomena, it is only now that the full 
(mean-field) critical dynamics in the fiber bundle model 
is solved. Apart from these, as already discussed, several 
aspects of the fiuctuations in this model are now well 
understood. Even from this specific point of view, the 
model is not only intuitively very attractive, its behavior 
is extremely rich and intriguing. It would be surprising 
if it did not offer new profound insights into failure phe- 
nomena also in the future. 
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